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Preface 


The treatment of quantum field theory given here is based on the 
Hamiltonian formalism and the theory is developed as far as possible 
in a logical manner, with results following from equations of motion. 
This involves a considerable departure from the usual treatment of 
quantum field theory in which one rather abandons the hope of logical 
deduction and is content with setting up working rules and showing 
their consistency, the equations of motion getting lost in the process. 

In physics one should aim at a comprehensive scheme for the de- 
scription of the whole of nature. A vast domain in physics can be 
successfully described in terms of equations of motion. It is necessary 
that quantum field theory be based on concepts and methods that can 
be unified with those used in the rest of physics. This necessity forces 
one to think of quantum field theory in terms of equations of motion. 
The usual treatment should thus be considered a stopgap, without any 
. * ng future. 

’. 1>se lectures were given at Yeshiva University during the academic 
year 1963-1964. The tape recording of the lectures was transcribed by 
Daniel Wisnivesky, to whom J am greatly indebted for giving so much 
time and care to the work. I have made a number of alterations to 
present the material in a form suitable for a book. 


P. A. M. Dirac 
Cambridge, England 
August 1966 
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§1. Relationship of the Heisenberg and Schrodinger Pictures 


In atomic theory we have fields and we have particles. The fields 
and the particles are not two different things. They are two ways of 
describing the same thing—two different points of view. We use one 
or the other according to convenience. In general the particle point of 
view is the more convenient one for problems involving only a few 
particles. It is the one usually used for elementary problems. When 
we deal with many particles of the same kind, the field point of view 
is the more convenient. For a general theory, which is to be formulated 
accurately, we have many particles coming in, maybe an infinite number, 
and the field concepts are useful. The present course of lectures will be 
concerned with general theory, so it will be based mainly on the field 
point of view. I shall assume that you have some knowledge of the 
basic ideas of quantum theory and also that you are familiar with 
relativity and Maxwell’s electrodynamics. 

Our object is to get a single comprehensive theory that will describe 
the whole of physics. This theory should consist of a scheme of 
equations, together with the rules for applying and interpreting the 
equations. The equations by themselves do not form a physical theory. 
Only when they are joined onto the rules for using them do we have a 
physical theory. A primary requirement is that the equations and the 
rules for using them should be consistent. 
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We need a single consistent comprehensive theory. Any special 
theory that one sets up for dealing with a particular problem should be 
consistent with this general theory. I do not need to tell you that such 
a theory has not yet been attained. It is the ultimate goal, toward 
which all physicists are working. These lectures are intended to show 
you how far one can go toward this goal at the present time. 

Quantum field theory is a subject that has been very much developed 
during the last fifteen years. The recent development started with the 
discovery of the Lamb shift and its theoretical explanation, a very 
brilliant discovery because of the close agreement between the theory 
and the experiment. Many papers have been written on the subject, 
and now there are quite a number of books about it, and people who 
want to learn about this recent work on quantum field theory will be 
better advised to follow some of the books instead of going through the 
original papers. 

This work, though, has serious difficulties in it. There are many 
quantities appearing in the theory which turn out to be infinite although 
they ought to be finite; and physicists have followed all sorts of tricks 
to avoid these difficulties, but the result is that the theory is in rather 
amess. The departures from logic are very serious and one really gives 
up all pretense of logical development in places. In this course of 
lectures I propose to follow different lines, which avoid some of the 
difficulties. I won’t be able to avoid all of them, but I will avoid some 
of the worst ones and I will, at any rate, be able to preserve some 
semblance of logic. By that I mean that I will neglect only quantities 
that one can have some physical feeling to count as small, instead of 
doing what people often do in this kind of work: neglect something 
that is really infinitely great. It will not be possible for me to set things 
out rigorously. I shall neglect things I shan’t be able to prove are really 
small, but still I hope you'll have some feeling that these things are 
small. 

The system of approximations I shall use will be somewhat similar 
to the approximations that engineers use in their calculations. Engineers 
have to get results and there are so many factors occurring in their 
problems that they have to neglect an awful lot of them; they don’t 
have time to study everything seriously and they develop a sort of 
feeling as to what can be neglected and what can’t. I believe that 
physicists working in quantum field theory will have to develop a similar 
sort of feeling as to what can be neglected and what can’t. The final 
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test is whether the resulting theory is coherent and in reasonable 
agreement with experiment. 

Well, those are the general features of the theory that I shall be 
presenting to you and I would like to say at the outset in what way 
my form of the theory differs from the usual form. In quantum 
mechanics as it is usually taught there are two pictures that one may 
use. There is Heisenberg’s picture, in which one works with Heisen- 
berg’s dynamical variables, and in which one has Heisenberg’s equation 
of motion, 


ih = uH — Hu. (1.1) 


H, the Hamiltonian, is the quantity that represents the total energy 
of the system. In the Schrédinger theory we work with state vectors 
that represent the states. We have a state vector representing a par- 
ticular state, denoted by a symbol such as |A), and it varies with the 
time according to the equation 


_d 
in|) = HA). (1.2) 


H is the same in both pictures. The essential difference between the 
two pictures is that in Heisenberg’s the dynamical variables vary with 
the time while the state vectors are fixed; in Schrédinger’s the state 
vectors vary with the time while the dynamical variables are fixed. 

The two pictures are connected in this way: any Schrédinger dy- 
namical variable is connected with the corresponding Heisenberg 
dynamical variable by the transformation 


us = BE ey se (3) 
and state vectors are correspondingly connected by 
|4s) = e*7"/"| Ay). (1.4) 


One easily checks that, with ug constant, (1.3) makes wy vary with ¢ 
according to (1.1), and with |A4;;) constant, (1.4) makes |Ag) vary with 
t according to (1.2). 

People usually say that these equations establish the equivalence of 
the two pictures—that one can use whichever one likes indiscriminately. 
However, the argument is valid only if e*”“/” exists and can be applied 
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to one state vector to give another, and for the Hamiltonians one meets 
with in quantum field theory there is good reason to believe that this 
is not so, because of convergence difficulties, and so the two pictures 
are not equivalent. 

The reasons for believing that the two pictures are not equivalent 
arise from examples for which one can find the Hamiltonian and be 
confident that it is correct. There is the example of quantum electro- 
dynamics, in which one considers électrons and positrons interacting 
with the electromagnetic field. One can set up the Hamiltonian for 
that problem and be pretty confident it is correct, because when one 
uses it in the Heisenberg equations of motion one gets sensible field 
equations, which are in close analogy to the classical field equations, 
are relativistic and look satisfactory from every point of view. On the 
other hand, when one uses that Hamiltonian in the Schrédinger picture, 
one gets a Schrédinger equation that has no solutions (at any rate no 
solutions of physical value). People have been trying very hard for 
very many years, decades even, to get solutions and they have failed, 
and now it seems pretty certain that solutions do not exist. At any rate, 
if they do exist they are of no use to us if we cannot find them. 

One would expect in any case that there should be a trivial solution 
corresponding to the physical state for which there are no particles 
and no fields present—the perfect vacuum. The vacuum is quite a 
trivial thing physically and we should expect it to correspond to 
a trivial solution of the Schrédinger equation, but there isn’t even 
this trivial solution. That makes one think that the Heisenberg picture 
is a good picture and the Schrodinger picture is a bad picture and that 
the two pictures are not equivalent. 

One can consider e~*”"/* to exist in a formal way and multiply it 
into a constant state vector |4). Then the product e7*#*/*|4) will 
satisfy the Schrédinger equation (1.2). The Schrédinger wave function 
consists of the coordinates of e~*¥/*|4). The nonexistence of the 
Schrodinger wave function may be attributed to e~*#*/*|4) not having 
coordinates, rather than to its nonexistence. It lies in a space of too 
many dimensions to have coordinates. 

The state vectors are usually supposed to be vectors in Hilbert space, 
in which case they would have coordinates. I should say here that lam 
using the term Hilbert space in the sense in which Hilbert originally 
used it—as a space with an infinite number of dimensions but still a 
separable space, that is, a space that can be spanned by an enumerable 
infinity of basic vectors. Mathematicians today use the term Ajilbert 
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space in a more general sense—to include also spaces that are not 
separable—but whenever I mention Hilbert space it is to be understood 
as a separable Hilbert space. The vectors in such a separable Hilbert 
space can be represented by coordinates. 

We are led to the conclusion that in quantum field theory the state 
vectors are not in Hilbert space. We may start off a particular state 
vector in Hilbert space at a particular time. Suppose we then make it 
vary with time in accordance with the Schrédinger equation: What 
would happen to it? Roughly, what happens to it is that it gets knocked 
right out of Hilbert space in the shortest time interval possible. The 
interactions that are physically important in quantum field theory are 
so violent that they will knock any Schrédinger state vector out of 
Hilbert space in the shortest possible time interval. It is thus not 
possible to get a solution of the Schrédinger equation for which the 
state vector stays in Hilbert space. So the Schrédinger picture is 
unworkable. 

The Heisenberg picture is a good picture, because the Heisenberg 
equations of motion are reasonable. They are not entirely free of 
difficulty; I am not asserting that we don’t get into any difficulties at 
all using the Heisenberg equations in quantum electrodynamics, for 
instance, but we do not get into such profound difficulties as we would 
if we tried to base everything on an equation that doesn’t have any 
solutions at all. 

We have the problem therefore of setting up quantum field theory 
working entirely in the Heisenberg picture. That is what I shall try 
to do in this course of lectures. We may take the existing theory given 
in the textbooks, for instance, and then cut out all references to the 
Schrodinger picture, and see how well we can get on without them. 
That is what I would ask you to do if you are going to use the text- 
books at all. Wherever you see references to the Schrédinger picture, 
cut them out. You will find that very often you can get on fairly 
well without these references, and that by cutting them out the whole 
theory becomes more logical and more intelligible. I would say that 
we are cutting out a lot of deadwood from the usual presentation of 
the quantum field theory. What is left, referring only to the Heisenberg 
picture, is what is essential and what we should concentrate our attention 
on. 

We have now to survey the usual quantum theory and see how much 
of it survives with this present point of view. If we just consider 
relationships at one particular time, you see that just about everything 
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survives. At one particular time tg we have dynamical variables a(t). 
These Heisenberg dynamical variables will operate to the right on some 
vectors we shall call ket vectors. The situation that we have at one 
particular time is essentially the same whether we are thinking in terms 
of the Heisenberg picture or the Schrédinger picture. The difference 
between the Heisenberg picture and the Schrédinger picture comes 
into play only when we vary fp. So all the usual quantum theory, insofar 
as it concerns relationships between dynamical variables and kets at 
one particular time, can be taken over and preserved in the new 
formulation. 

We have the kets at one particular time and we may picture them as 
corresponding to physical states at that time, so we retain the concept 
of a physical state at a certain time. We do not have the concept of a 
physical state throughout all time. The latter concept, which lies at 
the basis of the Schrédinger picture, will not occur in the present 
formulation. We have only the much weaker concept of a physical 
state at a particular time. It is a nonrelativistic concept, a very special 
thing which refers to one particular Lorentz frame. The Schrédinger 
concept of the state throughout all time would be a relativistic concept 
if it existed and would therefore be something one could think of as 
being rather fundamental in nature. But it does not exist, and all 
we have about states in the present formulation is a physical state at 
a certain time, a nonrelativistic concept—a very special thing, but 
useful in the calculations. 

We can preserve the usual quantum relationships between dynamical 
variables and kets when we keep to the one time fy). For example, 
take the equation a(fo)|A) = alA), where a is a number. If we had 
that equation, we could say that |A) represents the state at time fo 
for which the dynamical variable a at time fo certainly has the value a. 
The usual theory of eigenvalues and eigenvectors and their physical 
interpretation remains valid provided we stick to one time. 

There is one respect, though,.in which the present formulation differs 
from the usual one even when we think of only one time, and that is 
that the totality of kets which we have for our dynamical variables at 
one particular time to operate on do not form a Hilbert space. 

We cannot get solutions in the Schrédinger picture with the kets 
remaining in Hilbert space, and corresponding to that we cannot get 
solutions in the Heisenberg picture in which our dynamical variables 
operate only on the kets of a Hilbert space. We have to imagine that 
the dynamical variables of the Heisenberg picture operate on the 
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vectors of some space that is a bigger space than a Hilbert space. | 
don’t know the mathematical nature of this more general space—it 
is better to leave it unspecified for the present than make a guess which 
might later turn out to be unsuitable for physical purposes. I expect 
the mathematicians will make guesses at what the space is, but they 
might very well guess wrong, and they already have made a wrong 
guess when they thought that Hilbert spaces were adequate for physical 
purposes. It is better to concentrate on the physics and try to develop 
a physical theory that will be successful in examples. 

Question: Does the difficulty correspond entirely to the fact that 
there are an infinite number of degrees of freedom? Is that the sole 
cause? 

Answer: I won't say that it is the sole cause. You might very well 
have an infinite number of degrees of freedom and an interaction 
between them which is sufficiently soft that there exists a solution of 
the Schrédinger equation in Hilbert space. That would be possible, 
but that kind of interaction is not the kind of interaction that occurs in 
nature. For the kind of interaction that occurs in nature, when we try 
to get the solution in the usual way by a perturbation method, we get 
integrals that do not converge. We might set up a nonrelativistic theory 
with a different kind of interaction for which the integrals do converge, 
by putting in some sort of cutoff, for instance, or just an interaction 
that is not so big for the high frequencies. That would be mathematically 
acceptable if you don’t insist on relativity, and would lead to a 
Schrédinger equation you could solve, even though there is an infinite 
number of degrees of freedom with interaction between them. So the 
mere fact that you have an infinite number of degrees of freedom with 
interaction between them is not sufficient to throw out the Schrédinger 
picture. It is thrown out by the interactions which physicists are 
interested in being so violent in the high frequencies, and it doesn’t 
seem to be possible to get interactions satisfying relativity which do 
not have this violence in the high frequencies. 

I would agree that maybe physicists in the future will be able to think 
of a different Hamiltonian for which the interaction is not so violent 
at the high frequencies, and then the Schrodinger picture will come 
back into its own. That is quite a possibility. The theory I am presenting 
to you now I won’t say is the final theory. I would only say that for 
the kind of Hamiltonians we work with at present, which are successful 
in giving results such as the Lamb shift and the anomalous magnetic 
moment, one cannot use the Schr6édinger picture. It 1s better to 
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abandon all attempts at using the Schrédinger picture with these 
Hamiltonians. The best quantum field theory one can formulate at 
the present time is one that works entirely in the Heisenberg picture; 
but one cannot be sure that at some future time people might not think 
of a new Hamiltonian for which the Schrédinger picture wi// work. 

One has to preserve an open mind on these subjects if one is a 
research worker, and one musn’t believe in anything too strongly; one 
must always be prepared that various beliefs one has had for a long 
time may be overthrown. So please understand that what I am saying 
is, I believe, the best that can be formulated at the present time, but I 
certainly don’t want to assert that it will survive a long time. I don’t 
want to assert that the Schrddinger picture will not come back; in 
fact, there are so many beautiful things about it that I have the feeling 
in the back of my mind that it ought to come back. I am really very 
loath to give it up. That is why it has taken so many years to think 
that one may have to give it up. Still one has to face the realities of 
the situation at the present time, and they force one to give it up for 
the present. 


§2. Basic Quantum Concepts: Hilbert-Space Case 


In the quantum theory our dynamical variables are quantities that 
satisfy noncommutative algebra. We make the fundamental assump- 
tion that the quantum variables that correspond to the classical 
Hamiltonian variables g, and p, satisfy at any time the commutation 
relations 


q4s — 4sGr = 0, — PrPs — PsPr = O, 
(2.1) 
GrPs — PsGr = HE Orgs 


We define the quantum Poisson bracket of any two quantities u 
and v by 


ih[u, v]) = uv — vu. (2.2) 


With this definition we see from (2.1) that in the cases when u and v 
are any two of the quantities q,, ps, the quantum Poisson bracket has 
the same value as the classical Poisson bracket. 

With this definition of the quantum Poisson bracket we see that the 
classical Poisson-bracket laws hold also for the quantum Poisson 
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bracket; we have in the quantum theory, as in the classical theory, 
the laws 


{u,v] = —[, uv], [uy + ue,v) = [u,,v] + [ue, v], 


C3) 
[uiuvo,0) = uy[u2,v]) + [uy, vJute. 


In this equation one has to be careful to keep the order of uw, and us; 
that is, all through the equation, uw, should occur to the left of wo. 
Finally, we have the Jacobi identity 


[w, [v, w)] + [v, [w, v]]} + [w, [v, v]] = 0. (2.4) 


Now the question arises: What is the mathematical nature of the 
dynamical variables of quantum theory? In simple examples, one can 
consider them to be matrices. In fact, in Heisenberg’s original quantum 
mechanics they were simply matrices. They are matrices usually with 
an infinite number of rows and columns. That is, each quantity w 
can be represented as a two-dimensional array of numbers, where the 
number of rows and columns is infinite and there is a one-to-one 
correspondence between the rows and columns. Such matrices can be 
added and multiplied according to the standard rules of matrix algebra; 
and matrices are noncommutative, so they can be fitted in with the 
scheme of equations (2.1). 

Now if we look upon our dynamical variables as matrices, we may 
allow them to be multiplied into a column matrix, a matrix con- 
sisting of a single column, by putting this column matrix to the right 
of the square matrix. The column matrix can, of course, be con- 
sidered as a vector: It is just a set of numbers and we may look upon 
the numbers as the coordinates of a vector in some vector spacc. The 
column matrices that are important physically are the ones for which 
the sum of the squares of the moduli of all the numbers is finite; and 
then the column matrix corresponds to a vector in Hilbert space. The 
dynamical variables are then operators in Hilbert space. Each of them 
can be multiplied into one of these vectors and gives another vector 
according to a linear law. 

I should like to remind you that I shall always be using the term 
‘Hilbert space” to denote the space in which the vectors can be 
represented by coordinates as above. Using the term “Hilbert space” 
in this sense we see that in simple examples the quantum dynamical 
variables can be looked upon as operators in Hilbert space. 

There is another way in which we may look upon them as 
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operators in a Hilbert space. We may consider a row matrix and we 
may multiply it on the left into a square matrix by matrix multiplication 
and get another row matrix. The elements of the row matrix may 
also be considered as the coordinates of a vector in a Hilbert space. 
Then we have our dynamical variables operating to the left on the 
vectors of a Hilbert space, producing other vectors in that Hilbert 
space. We have two Hilbert spaces, one Hilbert space with the co- 
ordinates expressed as a row matrix, the other Hilbert space with the 
coordinates expressed as a column matrix, and our dynamical variables 
can operate to the left on the first kind of Hilbert vector and to the 
right on the second kind of Hilbert vector. 

A convenient notation is the following: The vectors on the right are 
called kets, and I write them |A), with a label such as A when necessary ; 
the vectors on the left are called bras, and I write them as a sort of 
mirror image of the ket vectors (B|. The reason this notation is con- 
venient is that we may multiply a row vector by a column vector, 
and then we get just a number (B|A); and this, you see, is a com- 
plete bracket expression. The bra vector is the first half of this 
complete bracket expression; the ket vector is the second half of 
this complete bracket expression; and the two of them multiplied 
together form a number. 

We may form quantities like (Blu|A) if we have associative multi- 
plication between our kets, bras, and dynamical variables. We shall 
assume associative multiplication for the quantities occurring in quan- 
tum theory. Mathematically it isn’t always valid. One can construct 
examples of matrices for which associative multiplication doesn’t 
hold, but these examples have not so far turned out to be important 
physically, and one can get on very well physically if one assumes 
associative multiplication. 

To introduce coordinates in a general way in a Hilbert space we 
may introduce a set of basic kets |r), r = 1,2,..., 0, with con- 
jugate bras (r|, to satisfy the following conditions: 


(1) They are complete; that is, any ket |P) can be written |P) = 
DEE ln. 

(2) They are orthonormal, that is, (rls) = 6,5. 

(3) They are independent; that is, there is no linear relationship 
among them. 

(4) Closure relation }°,|r)(r| = 1. 


One can assume (1) and (2) and deduce (3) and (4), or, alternatively, 
assume (3) and (4) and deduce (1) and (2). 
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In the Hilbert space the basic kets satisfying conditions (1), (2), 
(3), and (4) provide a system of coordinates. Any ket |P) has co- 
ordinates (r|P); any bra (Q| has coordinates (Q|r); and any dynamical 
variable a has coordinates (r|a|s), which form a square matrix. From 
condition (4) we deduce the law of matrix multiplication: 


(rloB|s) = >> (rlalt)(t|6|s). 


There are some dynamical variables which are represented by diagonal 
matrices: 


(rlols) = f(r) brs - 


There might be several such v’s, and they all commute. 

In a particular physical problem we look for a complete set of com- 
muting observables which we want to be represented by diagonal 
matrices. We then take their simultaneous eigenkets as basic kets: 


|v’) = lor’, Oo 


where |v,’) is an eigenvector of v, belonging to the eigenvalue v,’. 

We may introduce a second complete set w,, Wo,... of commuting 
observables w and a corresponding set of basic kets |w’). A ket |P) 
will then have coordinates (u’|P) and (w’|P) in the two systems of 
coordinates. These are related by 


(IP) = D7 &'|w')(w'[P). 


The above discussion in Hilbert space applies to simple cases and 
the question arises: What are the simple cases? Well, it seems that any 
dynamical system involving only a finite number of degrees of freedom, 
or any problem in which only a finite number of degrees of freedom 
comes into play, can be treated in this way. The quantum dynamical 
variables can then be looked upon as matrices or as operators on the 
vectors of a Hilbert space. 

But for quantum field theory, that’s not good enough; in quantum 
field theory we have an infinite number of degrees of freedom coming 
in and then in general it doesn’t seem to be possible to look upon our 
dynamical variables as matrices or as operators in Hilbert space. 
People have done a great deal of work in quantum field theory, on the 
assumption that the dynamical variables are, even then, operators in 
Hilbert space; this work leads to a lot of difficulties, many of which 
can be evaded by not making the assumption. 
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§3. Basic Quantum Concepts: General Case 


In the general case the question arises: What is the precise mathe- 
matical nature of the dynamical variables? It could be that they are 
just operators in the more general Hilbert spaces mathematicians now 
work with, but I’m not sure that is the case. I think it’s better not to 
have any preconceived ideas about it or we may be led astray like 
people were led astray for a long time by assuming that the dynamical 
variables are always operators in Hilbert space. It’s better, for the 
present, to keep an open mind about these dynamical variables and 
just call them g-numbers. 

g-numbers, of course, is just a name. Ordinary numbers, when one 
wants to make a distinction between them and these dynamical 
variables, may be called c-nuwmbers. The reason for the letters g and c 
are that g might make you think of quantum or queer, and c of classical 
or commuting. The essential difference between c-numbers and q- 
numbers is that c-numbers commute with everything, and g-numbers 
do not in general commute with other g-numbers. 

g-numbers obey the following algebraic axioms: the ordinary axioms 
of addition; the distributive axiom of multiplication u(v, -++ v2) = 
uv, + uve ; and the associative axiom of multiplication w(uw) = (uv)w. 
With these axioms we can establish the quantum Poisson-bracket 
laws (2.3) and (2.4), in analogy to the classical Poisson-bracket laws. 
But we do not have commutative multiplication, and if a product uv 
is equal to zero, it does not follow that either of the factors equals zero. 
In fact, this situation occurs also with matrices. You can very well 
have the product of two matrices equal to zero without either of the 
factors being equal to zero. Well, g-numbers are just used in con- 
nection with these algebraic axioms; calculations in quantum mechanics 
are just calculations working out the consequences of the basic equations 
(2.1) with the use of noncommutative algebra. 

When it’s a question of forming limits, we cannot proceed with 
mathematical rigor. Of course, we cannot have mathematical rigor 
without knowing the precise mathematical nature of the quantities we 
are dealing with. We shall continually be using g-numbers in integrals 
and differentiations, and we shall be making various assumptions which 
seem reasonable physically about the limits we use, but we’re not able 
to establish the theory with mathematical rigor. 

We proceed to build up the theory of g-numbers by making assump- 
tions that look reasonable as we need them. The first thing we need 
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is some standard of reality. In the classical theory we have a clear- 
cut distinction between real dynamical variables and complex dynamical 
variables involving the root of —1, and we need some corresponding 
clear-cut distinction between real g-numbers and complex ones in 
the quantum theory. The qg’s and the p’s we start off with in our 
quantum theory we assume to be real qg-numbers. In the classical 
theory we have the general rule that the Poisson bracket of two real 
things is real, and we make the corresponding assumption in the 
quantum theory: If w and v are real, their Poisson bracket is real. 
It follows that wo — vu equals a pure imaginary g-number except 
when it vanishes. We can infer that wo and vu are not both real, except 
possibly when wu and v commute. 

To preserve symmetry between the two products wv and vou, ’'m 
going to make the assumption that wv + vw is real, if uw and v are real. 
One can check that it doesn’t lead to inconsistencies; in fact, it is well 
known that it doesn’t lead to inconsistencies for the case in which the 
g-numbers can be represented by matrices. 

The conjugate of wv (let’s write it wv) is 


uo = dun + ou + Suv — vu = 4(uv + on) — B(uv — vu) = on. 
We can now work with complex g-numbers. Suppose 
E= 1, + ito, 1 = 01 + We with w,, wo, 01, V2 real. 
Then the conjugate complex of the product &y is 


En = (uy + ittg)(0y + ie) = Uy, — Valtg — Wolty — Wile 


(vo, — ive)(uy — ive) = 7 E. 


So to take the conjugate complex of a product, you take the conjugate 
complex of each factor and reverse their order. One finds that the rule 
extends to the product of more than two factors: Ent = €7% & and so 
on. We now have a standard of reality for our g-numbers, and we are 
able to determine whether any of them is real, or what its real and 
imaginary parts are. 

We assume that a g-number u can be multiplied into vectors on the 
right, called ket vectors and written |A), to give further ket vectors 
u|A), and into vectors on the left, called bra vectors and written (Bl, 
to give further bra vectors (Blu. A ket |A) and a bra (B| may be multi- 
plied, with the bra on the left, to give a number (B|A). 
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It is no longer possible to introduce coordinates for the kets and 
bras. If we tried to introduce basic kets and bras in accordance with 
the scheme of the preceding section, we should find that we could not 
satisfy conditions (1) and (4), because the number of dimensions is 
too large. We shall develop the theory by making suitable assumptions 
about the kets and bras which agree as closely as possible with the 
situation in the Hilbert-space case. 

In the Hilbert-space case each ket |A) with coordinates A, has a 
conjugate bra for which the coordinates are the conjugate complex 
numbers A,, which bra is labeled by the same letter A, so it 1s written 
(A|. We have the general relationships 


(BIA) = (AJB), (AIA) > 0. (3.1) 


We are going to assume that in the case when the vectors are not in 
Hilbert space and cannot be represented by coordinates, each ket |) 
still has a conjugate bra (A| and the relationships (3.1) still hold. 

We also assume that we have associative and distributive multiplica- 
tion for g-numbers and kets and bras and that for any product of q- 
numbers and kets and bras we get the conjugate product by reversing 
the order of all the factors and taking the conjugate of each factor. 

Suppose we take any g-number U with a reciprocal V. V is the 
reciprocal of U if and only if UV = 1 and VU = 1. Both conditions 
are necessary in a noncommutative algebra. If we have only UV = 1, 
we can infer UVU = U or U( VU) — 0. soreither } 2 — lon 
there exists € # 0 such that UE = 0. It is quite possible for the 
product of two g-numbers to be equal to zero without either of them 
being zero. 

Let us assume that V is the reciprocal of U; then we write VV = U—!. 
For any g-number @ we can then set up the similarity transformation 
oP == UR. 

It is easy to see that this transformation preserves the algebraic 
relationships between g-numbers. The transformation does not neces- 
sarily preserve reality relations. We should like to impose further 
conditions on U such that the similarity transformation preserves 
reality relations, that is, such that when a is real so is a*. 

Suppose that a = @ implies a* = a*. Then UaU7! = at = 
a* = UaU-1 = U-igU. Now U-1 = U-),s0 


Ua) = Uo 
or 


Be Basic Quantum Concepts: General Case 15 


We now have that UU commutes with every q-number, because any 
q-number a can be split into 


a= a, + lag, 


with both a; and ag real, and UU has to commute with both of those 
parts. A g-number that commutes with everything may be taken to bea 
c-number, so UU = k, a c-number. It is a real c-number because 
UU = UU. Furthermore, it is positive, because for any ket |A) we 
may form a new ket U|A) and then form (A|UU|A), which can be 
looked upon as the scalar product of the ket U|A) with its conjugate 
bra, so 
(A|[UU|A) = k(A|A) > 0, 


and as (A|A) > 0, we must have k > 0. 

Take now U’ = k—'/*U; then U’U’ = 1. Such a U’ is called a 
unitary q-number, and the transformation a* = U’aU'—! a unitary 
transformation of q-numbers. 

Let us consider now an infinitesimal unitary transformation defined 
by U’ = 1 + €A, where ¢ is an infinitesimal real c-number. Then 
U'~' = | — eéA, and for any g-number a we have 


a* = (1 + €A)a(l — €A) = a + €(Aa — aA) = a — ihefa, A]. 
We have the unitary condition 
l= UU’ = (14+ eA\(1 + €A) = 14+ €(A + A), 


so A + A = 0 and we can put A = iB, with B real. Then we can 
write the effect of the infinitesimal unitary transformation on a@ as 


a®* — a = hefa, B]. 


Now if we refer back to the Heisenberg equations of motion we see 
that the connection between a dynamical variable at time ¢ and that 
variable at time ¢ + 6f is just of this nature. For 6¢ small we have 


i452 — 8: = Stlg, fEi}h 


In classical mechanics the connection is an infinitesimal contact 
transformation. The contact transformations of classical mechanics 
correspond to unitary transformations in the quantum theory. 
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There is another way of writing 211.2, 


etffotlh 


ere 


8i+6t 
if we assume that 


et Hbtih — | 4 jHét/h for small 41. 


We can now integrate to get the connection between the dynamical 
variables at two times which differ appreciably: 


eth (ta—tih —iH (tg—t}){h% 


Say = §t1€ 


It can be proved formally that this formula is correct. However, 
it is not in general a useful equation, because of the difficulty of giving 
a meaning to the exponential. One could try to define it by a power 
series but one does not even know whether the series converges or not, 
beeause there is no general definition of convergence of g-numbers. 
In practice, to get the correction between g,, and g;,, one has to in- 
tegrate the Heisenberg equations of motion. 

In the present general case we cannot define a ket by coordinates. 
We may define a ket |S) by imposing suitable linear conditions on it, 
say 

us) = Waly) = he. eS) = ©, 


where the é’s are q-numbers, not necessarily real. These conditions 
need not be independent. They imply 


(es — eeaviie) = 0 or Gig ar Bees) = 0, 


which are the consistency conditions. 


84. The Fock Representation 
Let us suppose that we have a harmonic oscillator with energy 
UE Se a) 
Introduce the complex g-numbers 


a= (2h) "(p+ ig), B= (2A) Se a) (4.1) 


4. The Fock Representation AP 
Using the usual relation gp — pq = ih, we find 


so 7 Is like the operator 0/0n. 
We also have 


Bin = ig)p — 1g) = p + g —h = 2H — kh, 
(2 - 5) = a = Ga) = n(Z 4 5): 
or Hn = 7(H + nh). We can infer 
HEIGE -o= CAl ST EE te esa 
Let us now assume there is a ket |S) such that 


7|S) = 0. (4.2) 
Then 
H\S) = (nj + $)h|S) = 4h |S) 


and |S) is an eigenket of H belonging to the eigenvalue 44. Similarly, 
Hn|S) = n(H + h)|S) = $hn|S), 
Hy"|S) = (n + $)hn"|S). 


Thus the kets 7”|S) are eigenkets of H belonging to the eigenvalues 
(7 )n torn = 0,1, 2,..... 

One can easily see that if g and p are the only dynamical variables, 
there cannot be more than one independent ket satisfying (4.2). The 
way one would prove that would be to assume that there are two 
such kets, |.$,), |S2), with #|S,) = 0, %|S2) = O and assume that they 
are independent; and then introduce a projection operator w to satisfy 


wn”|S1) = ”|S1), wn” |S) = 0. 


One easily checks that it is quite consistent to assume the existence of 
such a w in the case when |S,) and |S2) are independent, One can also 
check that w commutes with » and 7, and hence with q and p, so it 
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is another g-number, another dynamical variable, independent of qg 
and p. That would contradict our hypothesis that g and p are the only 
independent dynamical variables. So ifg and pare the only independent 
dynamical variables, there is only one independent ket |S) satisfying 
(4.2). 

We may write a general ket |A) in the form 


where ¥(m) is a power series. We then have the function y(n) to rep- 
resent the ket |A). It is different from the usual kind of representation 
in that the variable y in the wave function is complex. 

The conjugate bra to |S) is (S|, satisfying (Sly = 0. We assume 
(SiS) 1 

If we have two kets 


|4) = Do ann"|S) and |B) = > dan"S), 
we have for the scalar product 


(BIA) = >> (SlbnW”"Gn 0” |S). 


nm 


We have 


(Sla'n™|S) = (S/R "(n" + ny” S) 

(Sly 'g"—"|S) = nla — 1)(S|y"—? 4”? |S) = ++ - 
_ (° fOr 77 sn, 

(al (OL) ae 


I 


and hence (BJA) = >,n!b,an. If we normalize (A|A) = 1, then 
Lenit!|an|? = 1 and n!la,|? is the probability of finding the system in 
the nth excited state for the state | A). 

One can make an interesting development, as was done by Fock, by 
introducing »~* to satisfy »—', = 1. We cannot have fg = I, 
because (S]y = 0, leading to (S|ny~! = 0. We then find S| 
(S|%”, because, as is easily checked, both sides give the same result 
when multiplied on the left into n"|S), or into 7”|S) with m > n, or 


a ne with m <n, provided we make the additional assumption 
a |S) = ©. 
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We now have 
(B| = (S| >) bam” = (S| Do alban” = (Sl6(n), 


say, and |B) is thus represented by ¢(n), a descending power series in 
n. The scalar product is now 


(BIA) = Qni) $ BCom) 2 


taken around the origin. 

If there exist other dynamical variables besides p and q, then we may 
construct kets of the type (4.3), in which the coefficients a, are not 
c-numbers but functions of the other dynamical variables. In that 
way we can get a more general set of kets built up from the one standard 
ket |S). For example, in the case I gave above, where there was just 
the projection operator w, we can introduce dynamical variables which 
don’t commute with w and which enable one to pass from |S,) to 
|S); and such dynamical variables could occur in the coefficients a, 
and would lead to more general kets obtainable from |S,) than the 
ones we should have if we took the a’s tc be c-numbers. 

Those are the basic ideas of the Fock representation applied to one 
degree of freedom. The Fock representation is useful for any problem 
in which $(p? + q’*) plays an important part. The Fock representation 
is always permissible, but it would only be useful if 4(p? + q) plays 
an important part in the calculations we are doing. A representation 
in quantum theory is like a system of coordinates in geometry; it is 
always permissible to take any system of coordinates that we like, but 
a system of coordinates is useful and convenient only when it is related 
in some special way to the things that are important for the problem 
under discussion. And here the Fock representation is useful when 
4(p? + q°) is important in the problem under discussion. Such 
problems are problems where we have an oscillator with energy 
3(p? + q”) and perhaps perturbations on this oscillator energy. 

This theory can be extended to several oscillators, several q’s and 
p’s, where for each q and p we have some quantity $(p? + gq”), perhaps 
with numerical coefficients coming in, playing an important part in our 
calculations. Then we introduce several 7’s, one referring to each p 
and q, so that we have a theory with the dynamical variables na, ja, 
the suffix a taking on different values to correspond to the different 
degrees of freedom. They will satisfy the conditions that the y’s all 
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commute with one another, 


faith — fea — 9; 
while 
Halb — Nd7a = Caer (4.4) 


If we form yea, we can call that , the degree of excitation of the 
nth oscillator. We define the standard ket |S) of the representation by 


the conditions 
Fas = 0 for all 7. (4.5) 


Those conditions are consistent with one another because the 7’s 
commute with one another. As before, we can show that if the y’s and 
Ws are the only dynamical variables in our problem, there is only one 
independent |S) satisfying the conditions (4.5). 

We can set up kets now of the form 


[P) = aa 1 2? |S 


where 11;’, M9’, ... are c-numbers taking the values 0, 1, 2,3,.... It 
is easily checked that na7a|P) = 1a’\P), from the commutations relations 
between the 7’s and 7's. That shows that |P) is an eigenket of each ng 
belonging to the eigenvalue ,’. We therefore get a ket that is an eigenket 
of the degrees of excitation of all the oscillators. 


We may form a more general ket |Q) by superimposing all these: 
o>, = Dae sae ny none oe -8 |S) 


summed over all the values of the n’, where the coefficients are c- 
numbers. 


§5. Second Quantization: Bosons 


The Fock representation leads to the subject of second quantization, 
which forms a cornerstone in quantum field theory. I expect most of 
you know about it already, and to those I shall just say that you can 
keep all your previous knowledge on second quantization provided you 
apply it only to the relationships at one particular time 79. Iam going 
to restrict the discussion here to one particular time and to develop 
the theory of second quantization at one particular time. 


5. Second Quantization: Bosons Bil 


Let us suppose that we have a particle that can exist in various states, 
at a certain time of course. The various states are described by eigen- 
values for a complete set of commuting observables. Now, in talking 
about a particle we have only a finite number of degrees of freedom, 
and all the usual work about the equivalence of the Heisenberg and 
Schrédinger representations is then valid. In particular, all the in- 
dependent states for one particle form a Hilbert space, so that in 
dealing with just one particle we have various states at a certain time 
corresponding to the vectors of a Hilbert space, and we can introduce 
coordinates for these vectors. We may use the usual notation in terms 
of the basic vectors |a'), |a”), |a?), ete. We may take the alphas here 
to be a complete set of commuting observables for the one particle, 


and the suffixes 1, 2, 3,... refer to eigenvalues of those observables. 
For an assembly of u such particles, we may have the state 
Jara”) loro”) ++ + Jew) 


obtained by just multiplying the kets for the individual particles. 
That is a possible state for the assembly. Such a state does not occur 
in nature; in nature we have all our particles satisfying either Bose 
statistics or Fermi statistics. If they satisfy Bose statistics, only sym- 
metrical states occur in nature. If they satisfy Fermi statistics, only 
antisymmetrical states occur in nature. 

Let’s take the Bose case first. For bosons a possible state is the state 
that we get by symmetrizing the above ket between all the particles: 


So) |”) oo [ler Slay%a2” - EDC) (5.1) 


where S is the symmetrizing operator, S = (u!)~'/?OP, P being a 
permutation operator. The numerical coefficient (u!)—'/? is inserted 
for convenience of normalization. This is a possible state for an 
assembly of bosons. If the number of bosons is finite, the states for 
them form a Hilbert space. 

One needs to check up on the normalization of the ket (5.1). One 
may assume that the original ket vectors are all orthonormal: 
(a*\a’) = 6°, so they form the base of a representation in Hilbert 
space. Under these conditions one finds that the squared length of 
Sla%ao +++ ay’) is 1 when the a,b,...,g states are all different, 
because all the terms in the sum (5.1) are then orthogonal to one 
another. If there are repetitions, some of them are not orthogonal. 
For example, if a and b are the same, the permutation that interchanges 
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aand b will give us a ket that is not orthogonal to the one we had before, 
but the same ket as before. The result is that we get a factor 1;’! m9’! --- 
in the squared length, where the n’’s are the numbers of states that are 
equal. n,’ is the number of a’s equal to ay, etc. 

Now let us suppose that the number of bosons wu is variable and that 
u takes various eigenvalues uw’. The totality of all the independent kets 
then consists of the following; first, there will be a ket corresponding to 
no bosons present, which we shall denote | ); then there are kets 
with one boson present, |a,), kets with two bosons present, Slate’) ; 
kets with three bosons present, S|a%a’a*); and so on. We can now set 
up a whole series of kets describing various states for the assembly of 
bosons where the number of bosons present is allowed to vary. 

A general ket is a sum of all these kets: 


a\) + agle*) + AapS|a%a”) fs, 


where the coefficients are c-numbers and the coefficients involving more 
than one suffix are symmetrical. 

Well, this is the way in which one describes the states for an assembly 
of bosons at a certain time. There is a connection between these kets 
and the kets we have for the set of oscillators described by the Fock 
representation. The connection is 


Slata® +++ a%) = 41"10"2 + +-(S) = game’ ** nglS), 5.2) 


where the n’’s are defined as before. (Don’t be disturbed by having the 
same letter S for the symmetrizing operator and the standard ket.) 
We may put these kets equal because those on the left have the same 
properties as those on the right. The kets on the left are orthogonal 
to one another, just as the kets on the right. The squared length of 
each of the kets on the left is the same as the squared length of the 
corresponding ket on the right. 

The fact that we have this connection means that the two dynamical 
systems are mathematically equivalent—one dynamical system being 
the set of bosons, the number of bosons being variable, and the other 
dynamical system being the set of oscillators. Those two dynamical 
systems are equivalent; they are just two ways of describing the same 
physical reality. It is this connection which establishes the equivalence 
between the wave and corpuscular theories of light. You may describe 
light either as waves, using the » variables, or as photons, using the 
kets |a) for the various states of the photons. 


5. Second Quantization: Bosons US: 


I should like you to note that each oscillator corresponds not to a 
boson but to a state for a boson. There may be any number of bosons 
in a boson state, and their number nv is the same as the degree of 
excitation of the corresponding oscillator. That is what we want 
physically. There may be any number of photons in a certain photon 
state; and the photon state corresponds to a Fourier component of the 
field, which is an oscillator. We have a boson state corresponding to 
an oscillator. 

Now a question that suggests itself is the following: Suppose we 
started off the discussion of bosons referring to a different set of basic 
states for each boson. Above we had the set of states fa,),|a2) °°: 
and so on. We might use a different set of states |6,),|62)--- and so on, 
where the 6’s are another complete set of commuting observables for 
one particle. For example, the a’s might be the position variables of 
the particles, and the 6’s might be the momentum variables. We can 
just as well formulate our discussion for the states of the assembly of 
bosons in terms of the 6 ket vectors instead of the a ket vectors. We 
have the a ket vectors satisfying the orthornormal conditions 
(a? |) = 8%, and we impose similar orthonormal conditions on these 
new ket vectors. I shall write a general one |64), because there is 
usually no one-to-one correspondence between the eigenvalues of the 
a’s and the §’s, and if we used the same index a in |8), it would imply 
such a one-to-one correspondence. Thus 


(IE) =o 


The kets for one particle are in Hilbert space; one particle is a system 
with just three degrees of freedom, or maybe an extra degree for the 
spin, and for such a simple particle the usual theory in terms of Hilbert 
vectors is quite correct. Both the a and the 8 vectors are in Hilbert 
space. The connection between them is given by the usual transforma- 
tion of coordinates in Hilbert space; one of the basic ket vectors 8) 
is connected with the basic ket vectors |w) by the equation 


6“) = |a*)(a*|6*), (5.3) 


with summation implied over all values of a. 

We may work with these new kets to discuss an assembly of bosons, 
and this implies setting up the ket S|8,4B82”---6u%) to describe u 
bosons. We may again introduce a set of oscillators described by a 
Fock representation, so that the states of these oscillators are con- 
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nected with the states of the assembly of bosons by the same formula, 
S\8 1482" - +> Bu) = nanz-- + ngl5). 


The question arises: What is the connection between the new set of 
oscillators and the previous ones? There are these two sets of oscil- 
lators, both of them offering a description of the same assembly of 
bosons. Well, the ket |.S') is the same in both cases; it refers to the state 
with no bosons present at all, and that, of course, doesn’t make any 
reference to what kind of ket vectors we use to describe a boson. It is 
evident that the connections among the 7’s are linear ones, so that 
there is a homogeneous linear connection between the new n’s and the 
old ones. Exactly what that connection is we can easily work out by 
considering the states of one boson. If we just have one boson, it may 
be described by the states |v”) = n,4|S) and also by |84) = n4|S). 
Now we have the connection (5.3) between the |a”)| and |84)| and hence 
nalS) = nala|84)|S). Here (a|8*) is just a numerical factor, so we 
get na = na(a|84). This formula shows us how the oscillator variables 
transform when we make a transformation of the basic kets to which we 
refer our bosons. One can sum up this result by saying that the y’s 
transform in the same way as the basic kets for one boson. 

When we describe one particle in the Schrodinger way we have a 
ket |P) to describe the state of one particle, which can be represented 
by a Schrédinger wave function, consisting of (a7|P), or (6“|P). Now 
the wave function transforms like a basic bra; the transformation law 
is 


(e4|P) = (64 


The transformation law for the wave function is therefore the same as 
the transformation law for the 7's. 

The »’s have the physical meaning of being emission operators of a 
boson, or operators that increase the excitation of an oscillator by 
one stage ; the #’s are correspondingly absorption operators, or operators 
which reduce the degree of excitation of an oscillator by one stage. 
So we see that the wave function for one particle transforms like the 
absorption operators in the Fock representation. 

This connection is known as second quantization. The reason for 
that name, of course, is that when you introduce the Schrédinger wave 
function you have performed a first quantization, if you start from the 
classical theory. Then second quantization consists in introducing 
Operators 7 to correspond to the wave function (a@"|P), and na, to 


at*){a"|P). 


6. Second Quantization: Fermions ae) 


correspond to the conjugate complex wave function (Pla). Second 
quantization consists of passing from the c-numbers, which give us a 
wave function and its conjugate complex, describing one particle, to 
g-numbers that describe an assembly of bosons. The g-numbers 
differ from the c-numbers only in the fact that the c-numbers all 
commute with one another while the g-numbers don’t, but satisfy 
Fanb — nba = Sap instead. So second quantization consists of intro- 
ducing a lack of commutation into the quantities that describe one 
particle in the Schrédinger theory, and this lack of commutation then 
provides us with a set of oscillator variables suitable for describing an 
assembly of bosons. 

I have written the formulas for the case when the basic kets for one 
particle are all discrete, but one can very well apply the same ideas to 
the case in which we have a continuous range of basic kets. The 6’s 
might very well be observables having a continuous range of eigen- 
values. That, in fact, is usually the case; if the 6’s refer to momenta of 
a particle or position coordinates of a particle, they certainly have a 
continuous range of eigenvalues. 

We have to make only a formal change in the equations when we 
have a continuous range of basic vectors. We have to replace the 
Kronecker 6 by the 6-function; thus for continuous kets we have 
(e4|83) = 6(64 — B®). When we pass over to the corresponding 
oscillators referring to the Fock representation, we shall get » variables 
having suffixes that take on a continuous range of values, and those 7 
variables will then satisfy the conditions 


Aang — neta = 6(64 — B”) 
Obi on Jobs abe Jobs = Gs )- >. 


86. Second Quantization: Fermions 


The above theory of second quantization can be applied to an 
assembly of bosons. There is a corresponding theory for an assembly 
of fermions. We may set it up in a way similar to the one we used for 
bosons. We start now with fermion particles, which can again be in 
various states; and we introduce, instead of the symmetrizing operator, 
an antisymmetrizing operator A defined by A = (u!)~'/?) + P, 
where the sum is over all permutations, with a plus or minus sign in 
front of each, according to whether the permutation is even or odd. 
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We can form Ala ;%a2°---a,%), which is antisymmetric between the 
various particles and which describes an assembly of fermions. It is 
now necessary that all the fermions be in different states; otherwise we 
get zero when we apply the antisymmetrizing operator. There is no 
longer any possibility of repetitions among the e’s; the number of 
fermions in any fermion state must be 0 or 1. 

In the case of bosons, the number of particles in any state is 
0, 1, 2, 3,..., any integer up to infinity, and that suggests a connection 
with oscillators, because the energy values for one oscillator are an 
arithmetical progression. That is why the assembly of bosons is con- 
nected with a set of oscillators. 

We want to do something corresponding to that for fermions, but 
it is certainly no good introducing oscillators, because the oscillators 
will have the wrong set of eigenvalues for their degrees of excitation. 
The number of particles in the various states are now variables, each 
of which has only two eigenvalues. That suggests variables like those 
which occur in the Pauli theory of spin one-half. We can develop the 
theory of fermions analogously to the case of bosons, by introducing 
operators 7 and 3%, which are not connected with harmonic oscillators 
but are connected with variables similar to the spin variables that 
describe a half-quantum of spin. 

Let’s introduce some complex qg-numbers 7, to satisfy 


nand + nona = O, Hand + Tota = O, 


naib + ona = Sad. ee 

For b = a these relations give 
mo = © (6.2) 
Naa + Fata = 1. (6.3) 


Well, the first question that arises is: Are these relations consistent? 
Is it possible to introduce such 7’s? One can easily check their con- 
sistency in the following way. We can consider matrices whose squares 
are equal to unity and which anticommute with each other. We have 
three such in the Pauli spin theory for an electron, and we have a 
larger number when we go to the relativistic theory of the electron. By 
extending the number of rows and columns in our matrices we can 
obtain any number of these quantities, which all have their squares 
equal to unity and anticommute with one another. Let us suppose 
we have a large number of such quantities & which all have their 
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Squares equal to | and anticommute; and let us label these £’s, so that 
there are two of them corresponding to each value of a. Call them Eq 
and &,’. Now take 


Va = 5(Eq ar 1), ja = glee Fi iE’). 


Then the y’s and 7's constructed in that way satisfy the required 
conditions, as is easily verified. 

We can treat the n’s and 7’s as we did in the Fock representation. 
We introduce the standard ket |S), satisfying 


7s) = 0 for all Ha, (6.4) 


noting that these equations are consistent with one another, because 
we don’t get into an inconsistency by forming 


HatolS) = 0, 
Weal) = 0, 
(Ha%» + Te%a)|S) = 0. 
We can now set up more general kets yansn-- + |S) and so on. All 


the y’s here have to be different. If two of them were the same, for 
example, a and c, we could bring the n, to the left of the factor 7») at 
the expense of a change of sign; and then we have the square of one », 
which is zero from (6.2). 

These kets can be ccnnected with the kets for an assembly of fer- 
mions. We may put 


Aley*ag"-++ oy") = nanp-* + 19/5). (6.5) 


The reason we can do that is similar to the reason it was permissible 
in the case of bosons. Two of the left-hand-side kets referring to 
different selections of states are orthogonal to each other. They are 
orthogonal to each other if they refer to different values of uw; and if 
they refer to the same value of u they are also orthogonal to each 
other if they refer to different sets of states. In the same way, the 
kets on the right side are orthogonal to each other when they refer to 
different » factors. One can easily check that. I am going to check that 
nanbne’** NglS) is orthogonal to nary +++ nq/|S), except when they 
have the same factors in front of |S). We form the scalar product of 
the conjugate of one of them with the other: 


(S| Figs "7° Ho Na’Nanb’ °° gl)» 
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We can change the order of the factors in accordance with the anti- 
commutation relations (6.1). If we take any one of the %’s, for instance, 
7), we can bring it over to the right and apply it directly to |S), and 
then it gives zero. The only condition under which we may have a 
nonzero result is when one of the 7’s is equal to y,-. Similarly, for each 
of the other 7, factors, there must be a corresponding 7; factor; other- 
wise the whole product is zero. 

We have therefore on the right side of (6.5) a set of kets that are all 
orthogonal to each other when they refer to different sets of suffixes. 
If we take two of them referring to the same set of suffixes in a dif- 
ferent order, then those two kets differ in sign or are equal to each 
other, according to whether the permutation that goes from one 
order to the other is odd or even. That is just the same situation that 
we have on the left side. 

The kets on the left side of (6.5) are all normalized. In the boson 
case they were not all normalized, because when we have more than 
one boson in a certain state, that brings a factor not equal to unity 
into the square of the modulus of that ket. But in the case of the 
fermions we cannot have more than one particle in the same state, and 
that situation cannot arise. So the kets on the left of (6.5) are all 
normalized. The kets on the right of (6.5) are also normalized, provided, 
of course, we make the standard ket |S) normalized. That is easily 
checked in this way. We form (S|%,°° + He%anans*** ng/S). Now let 
us take these two middle factors jana. From (6.3), we can replace 
them by | — xia. The factor %, here can be shifted over to the 
extreme right, as there cannot be any other 7,, and 7, then gives zero 
by acting on |S). So the middle factors can be replaced by unity. 
Then, of course, we have 7,7, as the two middle factors, and we treat 
them in the same way. They also give the factor unity, and so on. 
So the normalization is all right. 

Now we have proved all that we need to establish that the two sets 
of kets in (6.5) may be identified with each other. We now have a way 
of describing a set of fermions in terms of these operators 7 and 3. 
We have the ket |S) corresponding to no fermions present; 7/5) 
corresponds to one fermion in state a; yano|S) corresponds to two 
fermions in states a and b; and so on. 

We have in this scheme that the number of fermions in any state 
is given by the operator neg = Na. The proof consists in applying 
Na to each of the kets in the series above and seeing that we get the 
right result. 74|S) = 0, so |S) is an eigenket of N, belonging to the 
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eigenvalue 0. Nana|S) is equal to 
naNana|S) = na(1 a natia)|S) = TES): 


SO nal|S) is an eigenket of N, belonging to the eigenvalue 1. We check 
in the same way that Nanon-:-- |S) equals O if b,c,... are all dif- 
ferent from a; and Nanany:--|S) = nano: +: |S). The factors npn- 
don’t affect things here at all. That shows that for all the kets 
none ++ * |S), Na has the eigenvalue zero; and for all the kets nany -- - |S), 
N, has the eigenvalue 1. Thus N, does represent the number of fermions 
in state a. 

nq can now be looked upon as the operator of creation of a fermion 
into the state a. If we apply 7, to any ket representing a state with no 
fermions present in state a, the resulting ket refers to a state with one 
fermion present in state a. The fact that we cannot put two fermions 
into the state is connected with y,” = 0; if we apply the operator 7a 
twice we always get zero, no matter what we start with. Similarly, 
7a is the operator of destruction of a fermion from state a: 7q gives 
zero unless there is a fermion previously present; and then it gives a 
new ket for which this fermion has been removed. 

We can prove algebraically that the number of fermions in state a 
is QO or 1. For that we form 


= = Se = 22 => — 
Nae = Nalanata — na(l = Taha ie = Ucila — Wa te > Wale = Nae 


Thus NV,” = N,; 80 the eigenvalues of N, are 0 or 1. 

We may make transformations to different states from our original 
fermion states. We can do that just as we did for bosons. That will 
give us a new set of kets for the assembly of fermions, which can be 
connected with new 7 and 7 variables. The connection between the 
old and new variables is similar to that for the case of bosons. The 
calculation goes exactly the same for fermions and bosons. Thus we 
have a theory of second quantization for fermions. 

The final results for fermions and bosons can be connected in a 
single scheme of equations, 


nan —E nba = 0, Tih 2= Ta = bab OF é(a a b), 


where we have plus for fermions, minus for bosons. We have jaja = Nae 
which holds for both fermions and bosons. Reversing the order of the 
factors, Jana = 1 Na, where the minus sign is for fermions and the 
plus sign for bosons. 
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Those are the important equations connecting the emission and 
absorption operators for fermions and bosons. It is really remarkable 
that there is such a similarity between the cases of fermions and 
bosons. Looked at from the particle point of view, they are very dif- 
ferent; but expressed in terms of the operators 7 the equations become 
very similar in the two cases; there is just the difference in sign. 

We can make one development of the theory of fermions to which 
there is nothing corresponding in the theory of bosons. Let us go back 
to the equations (6.1) for fermions. We note that there is symmetry 
between the 7’s and the 7’s. Put all the n’s equal to 7s, and all %’s equal 
to n’s, and the equations remain unchanged. In the case of the bosons 
that is not true, because the sign goes wrong. 

Now let us see what this symmetry means physically. If we have a 
theory of an assembly of fermions, we can get another theory, in which 
we interchange the operators y and 7. Now 


Nata = ies Nana = t= Nie 


So if we interchange y’s and 7s, N, gets altered to 1 — N,. If the 
original NV, is equal to I, the new JN, is zero and vice versa. It means 
that if the original state is occupied, the state in the new description 
is unoccupied, and vice versa. We get, by this transformation, a new 
description in which the concepts of full and empty states are inter- 
changed. 

The original description would be useful if most of the N’s were 0 
and there were just a few particles corresponding to a few of the N’s 
being equal to I. But if we have a situation in which most of the N’s 
are equal to I and just a few of them equal to 0, we can make this 
transformation and get a new description in which most of the states 
are empty and just the few holes of the original description become the 
occupied states. 

Of course, we don’t have to apply this transformation to all the 
states. We might just apply it to some of them. For some values of a 
we might interchange 7, and },, and not for others; and then the rela- 
tions (6.1) are still all right. Thus just for some of the states we change 
the concepts of full and empty. Well, you know that that is what one 
does in the relativistic theory of electrons. Considering the electron 
as one particle by itself, studying its relativistic quantum theory, one 
finds that it has states of positive energy and states of negative energy. 
It becomes reasonable to interchange the concepts of full and empty 
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states for the negative-energy states, so that we concentrate our 
attention on the originally unoccupied negative-energy states. 

The ket |S) has to be changed when we make a transformation in- 
volving the interchange of » and %. The original |S) satisfied the con- 
ditions %,|S) = 0. When we interchange the n’s and %’s we are led to 
work with a new |S), satisfying 7a/Snew) = 0. |Snew) is a ket for which, 
with the original picture of particles, we have all the states occupied. 
Multiplying it by 72 we get a state with one hole. 

When we interchange y’s and 7s only for some degrees of freedom 
and not all, we need a ket |S) for which some 7’s applied to |S) are 
equal to 0 and some 7’s applied to |S) are equal to 0. The general 
|S) satisfies 


HalSgenn — © for some a-values, 


nb|Sgen) = O for b-values differing from the above a’s. 


When we are dealing with electrons and positrons it is convenient to 
set up such a ket for which the a-values include all the positive-energy 
states and the b-values include all the negative-energy states. 

Quantum field theory is formulated entirely in terms of these opera- 
tors 7 and %, referring to either bosons or fermions. It is true quite 
generally that all quantum field theories are built up in terms of such 
n’s and 7's; they are the basic bricks one uses. 
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The n’s and 7s are the g-numbers that come into the description of 
nature at one particular time. We shall now discuss how they vary 
with time. We shall assume that there is a Hamiltonian and that we 
have Heisenberg’s equations of motion giving us their variation with 
time. 

As far as the conditions at one time are concerned, we have the kets 


nano’ * |S), which lie in a Hilbert space if you impose suitable con- 

vergence conditions on the coefficients when you form sums 
Coen. ci |S). Ga) 
Ue wer 


One can in that way get an adequate description of the different states 
at one time, while sticking to one particular Hilbert space. 
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When we go over to conditions at another time, we would again 
have n’s and 7s referring to the new time, and we could again introduce 
|S) referring to the new time, and we would again get a set of Hilbert 
vectors. But these Hilbert vectors referring to the new time would 
not, in general, lie in the same Hilbert space as the set of ket vectors 
we started with. This is so for the Hamiltonians of interest in physics. 
The Hilbert space provided by the ket vectors changes, and if we want 
to consider the totality of ket vectors for all times we have to have a 
space larger than a Hilbert space. 

' This situation can best be discussed in terms of a model Hamil- 
tonian. I propose to take a model Hamiltonian that is not one which 
occurs in practice at all, but is a good model because it is sufficiently 
simple for us to be able to calculate explicitly everything that we want 
to calculate and yet is sufficiently complicated to show up the need 
for passing from one Hilbert space to another and the impossibility 
of representing everything in terms of one Hilbert space and using the 
Schrédinger picture. 

We shall work with fermions, so we have a set of dynamical variables 
"m » Tm Satisfying (6.1); we shalltakem = 1,2,..., «. It is necessary 
to have an infinite number of fermion states, as the difficulties do not 
arise when there is only a finite number. 

We take as the Hamiltonian, 


A= 4 lee ivenibs ae Dap dean (72) 


where the coefficients a are c-numbers. This Hamiltonian is real. It 
just involves quadratic terms in the emission and absorption operators. 

Now, because qman = — nm, any symmetrical matrix a will give 
zero; only the skew part of this matrix contributes anything, so we 
take Qn» = —Gym. We then have a Hamiltonian that we may also 
write 


eel => Cer es 
H= Bintan are Gretna): 


Now let us introduce a ket vector |S) to satisfy 7|S) = 0 for all m. 
This ket vector corresponds to no fermions present. 

Let us now form H|S). It is, of course, 44mn1m7n/S), the second 
term in H contributing nothing. This H|.S) just consists of terms each 
referring to two fermions being present. 

Now let us form H*|S). We get terms referring to the existence of 
four fermions and terms with no fermions at all. The terms con- 
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taining four fermions present are pretty straightforward and there is 
no problem with them. I want to concentrate on those terms which 
contain no fermions present at all. They come from 


1 ae. 
—44mn Thntin@ogtinta\ )- 


Using the anticommutation relations (6.1), we find that this is equal to 


= One oat On ae Wain) iy) 


iad | Gean@notata as Ginn ApghmNp( Sng 7 nan) || S) 
= SO nCuy tata S) = Son dnal Oma NaTim)|S) 
ae ly eels) = —4(Ga)|S), (3) 


where (_ ) indicates the diagonal sum of the enclosed matrix. The 
matrix @ is defined by (@)mn = Gmn. It is not the conjugate Hermitian. 

Now we may very well have the diagonal sum infinite even with a 
perfectly respectable matrix a, a bounded matrix, for example. We 
may have 


fe = =I) (7.4) 


for example. This could be achieved by taking as the nonzero element 
a= 0 1 : C2) 


This is a matrix satisfying the antisymmetry condition and is a very 
simple and well-behaved matrix, and yet we get an infinity in (7.3) 
and we have to say that H*|S) does not exist. 

Suppose we try to get a solution of the Schrédinger equation; i.e., 
we try to solve 


ee 
in =|) = HIA), (7.6) 


starting off with |A) = |S) att = 0. A formal solution of this equation 
is |A) = e *##/715) because, looking at things from the abstract point 
of view and assuming that one can differentiate the exponential ac- 
cording to the ordinary rule for differentiating an exponential that 


34 Lectures on Quantum Field Theory 


contains a numerical parameter, this would just satisfy (7.6), and it 
also satisfies the right initial condition. But, you see, if you try to give 
a meaning to the exponential as a power series, the terms of the second 
degree in ¢ already fail, as they involve the second degree in H and 
H?|S) does not exist. 

You might wonder whether e~*“*/*|S) does have a meaning, in spite 
of its being impossible to expand the exponential. Maybe the ex- 
ponential applied to the ket does exist, even though H® applied to the 
ket does not exist. Before giving a definite answer to this question, 
I should like to discuss the problem from the Heisenberg point of 
view. 

Take this same Hamiltonian and apply it in the Heisenberg picture. 
The Heisenberg equations of motion give 


Ss = 
ih ae = npll _ lelaes = = ol ete = Teepe 


I 


jl = = = — 
sn oon aa Tat) We =, Gen Cae ica Nplin)} 
= = 5 agi ma Oot) > — Gyn 7Fin : 
Then, of course, we have the conjugate equation 


Phase 
th Tip = ~—Agntn- 


Well, those Heisenberg equations of motion are very simple. The 
7s and 7's all vary linearly in terms of themselves. The equations can 
very easily be integrated. 

Writing the equations in matrix form, we have 


d a a7 
i ins oe = —1an, 
and then 
Ae 
i? We = aan. 


We see that the general solution is of the type 
nt) = An(O) + B7(0), 


where A and B are matrices whose elements are c-numbers and functions 
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of t. Then one gets 


Ch es dB _— 2 Be 
h 7 Uae: 7 = iaA, N= = GaAs 


Putting it like this, we get equations in terms of matrices instead 
of q-numbers. q-numbers are rather mysterious things, and we like 
to get them out of the calculation when we can and replace them by 
matrices involving c-numbers as their elements, which are quite well 
understood. 

Now let’s take the example (7.4). Then the elements of A and B 
just vary according to the law of cos (¢/h) and sin (t/h). We can 
easily fit them in to have the right initial conditions, at t = 0: 


dA_ 4 4B 


A=], B= 70, opie a ia/h, 
and the solution is 
nr(t) = cos (t/h)n-(0) + isin (t/h)G,s75(0). (7.7) 


There you see we have integrated the Heisenberg equations of 
motion, and we ran into no difficulty. The kind of infinity we had 
earlier is a difficulty peculiar to the Schrédinger picture, which one 
avoids by working with the Heisenberg picture. That kind of difficulty 
is typical of quantum field theory. In quantum field theory we deal, 
of course, with substantially more complicated Hamiltonians than 
this model I’ve dealt with, but we have the same general features that 
the Heisenberg equations of motion are easier to work with than the 
Schr6dinger ones. We may very well be concerned in practice with 
an example for which there are infinities even with the Heisenberg 
picture. J am not saying that the Heisenberg picture eliminates all the 
infinities, but it eliminates some that occur in the Schrédinger picture 
and is superior to the Schrédinger picture in that way. 

Now let us go back to the question I left unanswered: Does the 
ket e*#*/*|S) exist? Looking at things from the Heisenberg point of 
view, we have our 7’s varying with the time in accordance with (7.7). 
In the Heisenberg picture we may introduce a ket |So) to satisfy 
7n(0)|So) = 0. It is reasonable to assume that there is such a ket. 
If we can introduce kets in the Heisenberg picture, that is presumably 
one of them, about the simplest one that one can think of. It is the 
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same as our previous Schrodinger ket |S). But we may equally well 
introduce |S;) to satisfy 7(¢)|S;) = 0. One time is as good as another 
in the Heisenberg picture and if |So) exists, |S,) should exist also. 
What is the relationship between those two kets? Well, that is 
easily worked out. The formal solution of the Heisenberg equation is 


Tin) = erin. (Ojegie as 


Let us substitute this expression for 7, in the equation 7,(1)|S;) = 0. 
We get 


eu (Oe TSO) Ort O lemmas )S7 aa 


Let us compare this equation with the equation 7,,(0)|So) = 0. We 
see that 
Ce) = |So), 


apart from a numerical factor, a c-number factor. 

So |S,) = e*”#!5,), and it would seem that e’”“/"|S9) exists. The 
fact that we cannot expand the exponential means that the ket is not 
in the Hilbert space of kets (7.1) with |So) for |S). Here we have a 
ket that does not lie in that Hilbert space; apparently it does not lie 
in any Hilbert space at all. Thus we cannot get a solution of the 
Schrédinger equation if we keep to Hilbert space. We can get a solution 
if we are willing to accept kets in a more general space than a Hilbert 
space. But these more general kets cannot be described by coordinates. 
As far as is known at present, there is no way of introducing co- 
ordinates to describe the vectors in this more general ket space. The 
Schrédinger wave function consists of the coordinates of the ket and 
therefore the Schrédinger wave function does not exist. So we get the 
situation in which there is an abstract solution of the Schrédinger 
equation, but there is no Schrédinger wave function. 

Suppose we started off with a different ket, not this |S9). Could 
we perhaps have better luck with the Schrodinger picture? We might 
try to start off with a ket for which the Hamiltonian is diagonal. Well 
let’s see what happens then. I shall take the example (7.4) for the 
matrix a. It is good enough to illustrate the points we want to discuss. 

Let us go back to the solution (7.7) of the Heisenberg equations. 
The conjugate solution is 


A(t) = cos (t/h)7(0) — isin (1/h)an(0). 
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Suppose we form 
W(t) + a(t) = e*/"[n(0) + 2%(0)]. 
We may write it out explicitly with the matrix a given by (7.5), 


mit) + Het) = e/*[n1(0) + 7200)]. 
na(t) — F(t) = e*/*[2(0) — 71 (0)], 
na(t) + Ha) = e*/*[n3(0) + 74(0)], 
na(t) — a(t) = e”"[n4(0) — 73(0)], 


and so on. This way of writing things shows how the variables connect 
up to give quantities that vary with time according to the law e™/*, 
that is, quantities that just get multiplied by a numerical factor in the 
Heisenberg picture as times goes on. 

We might now set up a ket vector |S*) in the Heisenberg picture 
satisfying 


(n1 + 7e)|S*) = 0, (n2 — 71)|S*) = 0, 


7.8 
(nz + Fa)|S*) = 0, (na — %3)|S*) = 0, a 


and so on. These equations remain valid for all time with the same 
|S*), because the operators on |.S*) just get multiplied by numerical 
factors as time goes on. We should check, of course, that the equations 
(7.8) are consistent. It won’t do setting up equations for a ket which 
are not consistent with one another. Let us take the first two and let 
us form the anticommutatort 


(jee fees — Tike — —In, ily + (he, 72], = 0. 


The consistency of the others is checked in the same way. All these 
equations are therefore consistent and (7.8) is just as good a way of 
defining a ket as any of the previous ones. 

So there we have defined a ket |S*). Now how would that ket 
behave in the Schrédinger picture. Well, we must express the Hamil- 
tonian in terms of the variables occurring in (7.8) and their conjugates. 
We can do that in this way: Using the given value (7.5) for the matrix a, 


tWe use the notation, 
[u,v], = ue + cu, 
[u, v]_ uv — va, 


[u,v] = (av — vn)/ih. 


MI 


| 
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H will contain as the first two terms 

ning — Hit2- 
We can express them in the following way: 


nine — Fite = 31m — Fe)(ne — 71) — (ne + Fd F 72) 
+ 91%1 + Fone + noe + Fin] 
= U(n1 — H2)(m2 — F1) — (no + AVC + F2)] + I. 


Let us now replace H by another Hamiltonian, 


Ay, = 4Qn1 — Fe)(n2 — 71) — (no + FDO + 72)] 
+ 4(same with 4394) + °°. 


H, differs from H by a c-number. The c-number has a contribution 1 
from each pair of fermion states, so the total difference is an infinite 
c-number. 

From the point of view of the Heisenberg picture, H is equivalent 
to H; they both give just the same Heisenberg equations of motion. 
From the point of view of the Schrédinger equation, H, and H are 
different. Now, H,|S*) = 0, so |S*) is a solution of the Schrédinger 
equation for H, . If we try to work with H we should have |S*) changing 
by a phase factor and changing infinitely rapidly. If you want to 
interpret the Hamiltonian as the total energy, you would have to say 
that H, and A correspond to two different physical systems in which 
the total energy differs merely by an infinite constant. From the 
physical point of view, H, is at least as reasonable to work with as H, 
and if we work with H, then we do have a solution of the Schrodinger 
equation. Not merely a formal solution but a solution that stays in 
Hilbert space; in fact, it doesn’t move at all. 

Well, then the statement, “One can’t get a solution to the Schrodinger 
equation in Hilbert space” really needs modification. One should 
put it this way. One cannot get a solution approximating to |So). 
One can get a solution differing very substantially from |S) if one 
suitably defines the zero-point energy. 

That is the mathematical situation. The question is: What are the 
important features that arise when we want to apply this theory to the 
Hamiltonians that occur in practice? We have our Hamiltonian 
provided for us and we want to get solutions that are physically useful. 
Solutions that are physically useful will be those close to the vacuum 
state, with just a few particles present, if any. Solutions close to the 
vacuum state correspond to kets close to |So), for which one fails to 
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get a solution to the Schrédinger equation in Hilbert space. If you take 
quantum electrodynamics with the Hamiltonian I shall obtain later, 
you could get solutions corresponding to |S*), but that kind of solution 
would not be of any help to the physicist, because it would differ too 
much from the physical vacuum. 

I am really anticipating future lectures when I talk about the Hamil- 
tonian for quantum electrodynamics, but I would like you to have a 
clear picture of what the situation is. You know that in applying the 
second quantization theory to electrons you have to suppose that the 
physical vacuum is the state for which all the positive-energy electron 
states are unoccupied and all the negative-energy states are occupied. 
Such a state is like the |S») that occurs above. You would get a state 
like |S*) if you supposed that all the electron states are unoccupied, 
both the positive-energy and negative-energy ones. You might start 
off with the mathematical situation in which all the electron states are 
unoccupied to begin with, and you could treat the ket corresponding to 
that state as we did |.S*). This |.S*) didn’t change at all in the Schrédinger 
picture. But going over to the correct Hamiltonian of electrodynamics, 
we get an |.S*) that changes a little but not seriously. However, that 
ket is not one a physicist can use, because it differs too drastically from 
the kets that describe physical situations. An example of such a ket 
for electrodynamics is given in my paper, Comm. Dublin Inst. Adv. 
Studies, Ser. A (No. 3) (1946), Sec. 3. 

To say that the Schrodinger equation has no solutions is putting 
things too drastically. In the first place, it has symbolic solutions that 
don’t stay in Hilbert space; second, it does have some solutions that 
stay in Hilbert space but differ too much from what the physicist wants 
to be of any use. Under those conditions I would maintain that the 
Schrodinger picture is unworkable and the Heisenberg picture is the 
one we have to restrict ourselves to. 

I’ve given this discussion in terms of a very primitive Hamiltonian. 
One could develop the formalism to work with rather more complicated 
Hamiltonians which illustrate the same features. There is a more 
powerful formalism given in a paper of mine (The Vacuum in Quantum 
Electrodynamics) in Nuovo Cimento, Suppl. 6, p. 322 (1957). In Section 4 
of this paper I prove a number of theorems that are useful for getting 
kets which don’t lie in a Hilbert space. I shall just enumerate the 
theorems. I shall use y’s instead of y’s to agree with the notation of 
the paper, so we have 


Wave + vava = 9, vate + Weta = Sap. 
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In this work w will be considered as a column matrix and ¥ as a row 
matrix. Let us introduce a square matrix \ whose matrix elements 
are c-numbers. 

One of the theorems is 


ly, PrAyJ_ = ry. (7.9) 


Here WAy is a single element, because is a square matrix, and when it 
is multiplied by a column matrix on the right and a row matrix on the 
left, the result is a single element. The right side of (7.9) is a column 
matrix, because we have here a column matrix y multiplied on the 
right into the square matrix \. One checks the theorem by using the 
anticommutation relations, quite a straightforward piece of work. 
There is the conjugate theorem to (7.9), which is 


ly, Pry]. = —p. (7.10) 
When you have proved those two you can easily go on to prove 
uv, Py]_ = Yu, 1_¥, (7.11) 


where w is another square matrix with c-number elements, like \. 
One can go on from that to prove that 


oY ye YW = ey. (7.12) 


Here Yay is a single element, and we can form its exponential in the 
usual way. There is a conjugate theorem 


oY We = ye™ (a3) 


Those are the theorems that are useful for the present kind of work. 
Suppose a ket |X) satisfies 


(ay + Pb)|X) = 0, (7.14) 


where a is a row matrix of c-numbers and } is a column matrix of c- 
numbers. Define 7 


[Y) = ex). (7.15) 
Then from (7.12) and (7.13), 
(Gene ve*b)|Y) = 0. (7.16) 


It takes a whole lot of equations (7.14) to define the ket |X). These 
equations might be, for example, that all the absorption operators for 
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a positive-energy electron applied to |X) give zero and all the emission 
operators for a negative-energy electron give zero. Then we are able 
to construct another ket | Y) from it, satisfying a different set of linear 
conditions (7.16). These conditions (7.16) would have to be con- 
sistent if the conditions (7.14) satisfied by |X) are consistent. That iS, 
if the operators ay + bd all anticommute with each other, the operators 
ae *y + Jeb must necessarily all anticommute with each other also. 
Thus we get a set of consistent conditions to specify the new ket | ¥). 

It might very well be that if we try to expand the exponential in 
(7.15) we get infinities coming in, just as with the more primitive 
example we were considering earlier. Under those conditions the ket 
|Y) would not be in the same Hilbert space as the ket |X). That is a 
fairly general way of constructing kets that lie outside Hilbert space. 

When I was working on this some years ago, I was hoping to be able 
to get a ket lying outside the usual Hilbert space but suitable to define 
the physical vacuum at all times, and so satisfy the Schrédinger equa- 
tion. That proved to be impossible. One can get the ket to satisfy the 
conditions up to a certain degree of accuracy, but going to greater 
accuracy the conditions fail. It has led me to believe that it is impossible 
to get a ket to describe the physical vacuum at all times. 


$8. Value of the Classical Theory 


We must now face the problem of finding the right Hamiltonian to 
be used in our quantum theory to get a description of nature. For this 
we work from a classical analog. This is simply because it would be too 
difficult to think of a quantum theory straight away. We want to have 
something to help us. The only thing that we can have to help us is a 
classical model. 

We must put the classical theory in a form suitable for passing to 
the quantum theory; the form that we need is the Hamiltonian form. 
Now you might say: Can we not immediately start off with the Hamil- 
tonian form of classical theory and work from that to a quantum 
theory? That would not be convenient, because we want our theory 
to be relativistic, and for the Hamiltonian form of classical mechanics 
it is not very easy to see whether the theory is relativistic or not. 

It is best to start from a classical action principle. It is very easy to 
see whether an action principle is relativistic or not: If the action 
quantity is Lorentz-invariant, everything has to be relativistic. So we 
see that it is best to start off from an action principle. 
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We then make an assumption of localization: that each bit of action 
is localized in time. This is sufficient to provide a Lagrangian, which is 
just the action per unit time. For a relativistic theory, we should also 
require each bit of action to be localized in space, so that it becomes 
localized in space-time. 

There are standard methods for passing from the Lagrangian to the 
Hamiltonian. These methods are discussed in my book Lectures on 
Quantum Mechanics, Belfer Graduate School of Science Monograph 
Series No. 2. 

Let us consider now the connection between classical theory and 
quantum theory; and let us take the example of a classical theory which, 
when expressed in Hamiltonian form, does not involve any con- 
straints. I want to leave the question of how to handle constraints 
until later and for the present just consider a classical Hamiltonian 
theory without constraints. To quantize this theory we take the dy- 
namical variables g and p occurring in the classical Hamiltonian theory 
and make them into g-numbers. We have a Hamiltonian function in the 
classical theory, and in simple examples we can assume that we have 
the same Hamiltonian in quantum theory. The Hamiltonian of the 
classical theory is a certain function of the classical dynamical variables 
and we set up the quantum Hamiltonian which is the same function 
of the quantum dynamical variables; but I would like you to note 
that this procedure of quantization is not uniquely determined. Given 
a classical Hamiltonian, we cannot in general say with certainty what 
the corresponding quantum Hamiltonian is. The classical Hamiltonian 
may involve a product of two factors that don’t commute; now, when 
we go over to the quantum theory, what are we to do about this 
product? The classical theory doesn’t tell us the order in which to put 
the two factors; we have to decide on an order when we go over to the 
quantum theory. 

It is, of course, necessary that the quantum Hamiltonian be real, 
as it is the total energy of the system; that’s something physicists can 
measure, and has to be a real thing. The quantum Hamiltonian has 
to be a real g-number, but this condition of reality is not in general 
sufficient to fix which quantum Hamiltonian corresponds to a par- 
ticular classical Hamiltonian. Given one real quantum Hamiltonian, 
we can add to it # times any real quantity, and get another quantum 
Hamiltonian which corresponds to the same classical one. 

The result of this is that when we’ve got a given classical theory, 
in general there is not a unique quantum theory corresponding to it. 
There is no well-defined unique process for passing from classical 
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theory to quantum theory. That means that when we set up a quantum 
theory we have to set it up to stand on its own feet, independent of the 
classical theory. The only value of the classical theory is to provide 
us with hints for getting a quantum theory; the quantum theory is then 
something that has to stand in its own right. If we were sufficiently 
clever to be able to think of a good quantum theory straight away, 
we could manage without classical theory at all. But we’re not that 
clever, and we have to get all the hints that we can to help us in setting 
up a good quantum theory. 

The classical theory can help us quite a lot. It provides us with 
suitable Hamiltonians with which we can start to work in the quantum 
theory. We can study them and see if they are good or not. When 
you’ve got one of the Hamiltonians and you’ve studied it, you might 
find that you have to modify it; but still you’ve got a start. Without 
the classical theory you wouldn’t have this start at all. 

I think we may ultimately reach the stage when it is possible to set up 
quantum theory without any reference to classical theory, just as we 
already have reached the stage where we can set up the Einstein gravi- 
tational theory without any reference to the Newtonian theory. But 
from the point of view of teaching students, I think one would always 
have to proceed by stages—not expect too much from them, teach 
them first the elementary theories and gradually develop their minds; 
and that will always involve working from the classical theory first. 

Well, as far as constructing the Hamiltonian goes, it is just a question 
of hints. But there is a transformation theory of classical mechanics 
which is closely connected with the transformation theory of quantum 
mechanics; that connection arises from the similarity of the Poisson 
brackets in the two cases. The similarity of the Poisson brackets is 
something stronger than hints; it’s a very close connection between 
classical theory and quantum theory. You can’t modify the definition 
of the quantum Poisson bracket; that is uniquely determined by 
analogy to the classical Poisson bracket. 


§9. The Scalar Field 


As a simple example of a field that may be useful for the description 
of nature, we consider a real scalar field V, associated with particles 
of zero rest mass. 

The coordinates of space-time will be denoted x,(u = 0, 1, 2, 3). 
We raise and lower suffixes according to the rules: For any contra- 
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variant vector A,, we take 49 = A°, A, = —A(r = 1, 2,3). The 
derivative of a scalar function of position such as V will be written 
aV/ax, = V'". People who work in general relativity usually write 
the coordinates of space-time as x“, which would make all the suffixes 
opposite to the present notation. 

Beginning with the classical theory, we take for the action density 


Ge erm Cay 
It leads to the field equation 
t= 0 or [|V = 0, CLe 
where [_] = 07/dx, Ox". 
The Lagrangian is 


L= [ec CPs = en | vv ax 


(81) if CI WIE GR a, 


This equals kinetic energy — potential energy, as is usual in particle 
dynamics. 

The dynamical coordinates are the values of V at all points of three- 
dimensional space at a certain time. We may write them V,, where 
the suffix x denotes x; , X2, x3. The derivatives in space-like directions 
V’" are functions of the coordinates. The velocities are V, = V®. 

The momentum U, is the functional derivative 6L/6V,. To define 
this functional derivative we vary the velocities in L and pick out the 
coefficient of 6V,. We have 6L = (4r)~!{VéVd°x, and hence 
U = V/4xr. It is more convenient to redefine U equal to V, so that 
it is 4x times the momentum conjugate to V. 

With the new definition of U, the classical Hamiltonian for the 
scalar field is 


H= (dry | ov ax —L= en | w + VV") dx, (9.3) 


Note that H is positive-definite. 
We have the Poisson-bracket relations 


[V2 Uz] = 4a 6X, — x1’) &Xo — Xo") 6(x3 — X3/) 
= 4r 63(x — x’), 
o Loe U,"] = 0. 


[V 9 Vz") 
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They are the natural generalization of the Poisson-bracket relations 


for particles, applied to the infinite number of degrees of freedom of 
the flelae 


The Hamiltonian equations of motion are 


2 = [V., H) = [Vz, (81)! fu. d x] = U,, 


Us 


[oat Ce: | Vier Vat xa 


which reproduce the field equations (9.2). 

We shall now make a Fourier resolution of the field. Let V(x) = 
V(x0, X1, X2, X3), to be distinguished from V, = V(x1x2x3), in which 
a particular x9 1s understood. The field equation (9.2) shows that 
we can put 


V(x) = | (vee aL Pie Ve, 


where 
fe len I oxo = igXa — Korg — (Ks x) 
and 
ko = (ky? + ko? + ky?)'? = |k| > 0, 
and 


d’k = dk, dks dks. 


Here V;,° means V;.,x,<, and is a constant Fourier coefficient. 
For a given xo, 


V,= | (Vie 2? + Vue) Pk = i (Vi + Viaje" ak, 
(9.4) 


where V; = V;‘e**o%o is a time-dependent Fourier coefficient. Also 


U, = (~) > i fixlcr. = Vac |) d°k, (9.5) 
XZ 1L_xr3 


0X9 
We can invert these equations and get 
Voie = On)” | ee oe 
. (9.6) 
ilk|(Vi ~ V_«) 


(20) ° | Ure ee dx 
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V, and U, = dV,/dxo at a given Xo fix V, for this x9. The important 


thing is that the equations connecting V,, U, with V;, V_, do not 
involve x9. So we can consider Viz, V_; as new dynamical variables 


and they will satisfy V;, = [V,, H] and not 
(OV;,/dt) + [ee Hi, 


The relationship of the V, and V;,° should be clearly understood. 
The V; are the important quantities, of the nature of Hamiltonian 
dynamical variables. They and their conjugate complexes V, at any 
time xq are defined in terms of the field quantities at x9 by equations 
that do not involve x9. If the field V is interacting with other things, 
one may still use the V;,, V, as Hamiltonian dynamical variables. 
The V;°, on the other hand, are integrals of the equations of motion, 
not working dynamical variables. They are integrals only for the 
case when the field V is not interacting with other things. They cease 
to be well defined if the field is interacting with other things. 

Next we get the Hamiltonian in terms of the new variables V;, V;,. 


We have 
[ors = — fff - roe 
X Var — Viele ? Bk dk’ ax 


= || [k| |k’'|\Vie — Vie) Var — Vez) 
X 89(k + k’) d’k d’k! 


Vi 


I 


l 


I 


82° [Keo = Vl) Sa a 


/ ver’ Bx 


- | {fam + Vie" ky Vir + Var) 
K eT BE Bk! Gx 


= 87° / k?Ve + Vii\QV_, + Vy) d®k, 


and 
H = (ny [cw PE Soe on? [AP + V_4V_,)a®k 
or (9.7) 


H = 4x | V7, d°k, (9.8) 
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For the Poisson-bracket relations, we have from (9.6), 


ee ee Fe — 0 (9.9) 
ee V2, — 0 (9.10) 
We (eV Vee = Vg — ony ff tre, Upe™ =) 
Mode d «! 
= (162°)? / i| CRT eee Nise) 
id? x dx’ 
= CGF" / Cee on 
= Om) Ge Se 1), Ol) 
From (9.10) and (9.11) we find 
[Vis Vir — Viel = ie ran d3(k + k’). 
From (9.9) and (9.11) with & and k’ interchanged we have 
[Vis Ver + Vel = i ea 53(k + k’). 
Thus finally 
[Vi, Ver] = 0 (9.12) 
ee a om 4 aay 53(k + k') 
or 
[Vis Vier] = Te aa OAS == 1a (9.13) 


We may check the equation of motion: 


Vi, = (Vi, H] = 40°k? = Vn. 


i 
Ve 42k 
which agrees with the definition 


V, = Vice ele, 
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We may quantize the theory according to Bose statistics. This 
involves replacing (9.12) and (9.13) by 


ViVi — ViVi = 0 ae 


and 
—h 


ViV ar = ViVi = 4r2\k| 63(k = k’). (9.15) 


We compare these expressions with the ones coming from the Fock 
representation for a set of oscillators, 


Gets — tee — Ost Oe 


We see that V; corresponds to an emission operator 


af #Nt? 
= OS (i) nk 


and V; to an absorption operator 


af bs h 1/2 
Vs = (Oey (3) 7h » 


We cannot have V; corresponding to the absorption operator and 
V,, to the emission operator, because the right side of (9.15) is a negative 
quantity when not zero. 

We have in the quantum theory, as in the classical theory, 


H = ea fw HE OVERS = an? [eM i 
from (9.7). In the quantum theory this equals 
an? [evs a°k + infinite c-number. (9.16) 


In the Heisenberg equations of motion we may use for the Hamil- 
tonian either one of the two expressions 


H = Gay [cw +") dx or an? [eV a°k, 


but whereas the first equation gives an infinite contribution to the 
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vacuum energy, the second applied to the vacuum state |0) satisfying 
Viioy=o gives zero and is thus the physical energy. 
If we quantize according to Fermi statistics we put 


ViVi + VeeVe = 0, ViVi + ViaeVa = rae — k’). 
The expression H = (8r)—! f (U2 + V7") d°x for the Hamiltonian is 
no longer any good, because it does not give the correct variation of 
V;, with time. It is then not possible to localize the energy. For this 
reason the quantization of the scalar field according to Fermi statistics 
is not of interest. 


$10. The Scalar Field: Relativistic Invariance 


The question arises whether the quantum theory of the scalar field 
with Bose statistics is really relativistic. The classical theory is certainly 
relativistic. 

There are new assumptions made when one goes over from classical 
theory to quantum theory. One cannot be sure that the new assumptions 
are relativistic, because they refer to Hamiltonian variables, and all 
the Hamiltonian variables are defined in a nonrelativistic way. They’re 
defined as referring to one particular time. We have to check whether 
our quantum theory is really relativistic. We use our classical theory 
only to suggest the quantum theory, and then we have to check whether 
that quantum theory is satisfactory by making calculations entirely 
in the quantum theory. 

In the present case we have to check whether our commutation 
relations are independent of the Lorentz frame of reference with which 
we are working. We can do that in a direct manner by working out 
the commutator of V at two different points in space-time and seeing 
whether the result is relativistic. For the present let us take the case 
of the scalar field by itself. If there were other fields interacting with 
it, the calculation would have to be on different lines. 

We have 

[View Ver] = ian? |kl)~! d9(k — K’). 


For our present calculation it is convenient to work with the constant 
Fourier coefficients and not the time-dependent ones; they are related 
by Vio = VyeW*!*!*0, and, of course, Vz° = Vze'!*!%0, Let us work out 
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the Poisson bracket of two of these constant Fourier coefficients: 
Vie, Vu] = (Ve, Purlemi*ltocil* 20 = i(4e?|k|)~* 33(k — kK’). 


The result is that we get here just the same formula for the constant 
Fourier coefficients as for the time-dependent ones. And of course, 
we have (Mee, V,.°] = 0. Noncommutation appears only when we 
have one V and one V. 

Let us now take 


(Vix), V(x’)] = i) (ae ee 2 = Hisar Vine ae ae Vico 2) 
Maka 
Here k - x is the four-dimensional scalar product. We get 
V(x), Vx’)) = cua 
A Nasr |) a 
x (ete ee se a aa 


= taf vee em ee 
2 


Here kg means |k|. Now the convenient way of working this out is 
to rewrite it 


etn __ Cae) d(Ko _ \k|) d'‘k, 


where kg is now an independent variable. Now let us change the sign 
of all the k’s in the second term; this implies giving negative values 
to ky, but there’s no harm in that, because kg is now an independent 
variable. So we get 


i —! 3h-(2—" 
ia} Ie ae om [olka 


That’s the result and we want to see whether it 1s Lorentz-invariant. 

Suppose we have any four-vector with components ao, @,, @2, 43. 
Then 6(a,a") is Lorentz-invariant. The function 6(a*a,) splits off 
into two parts, one referring to positive a) and one referring to negative 
dg; and we may change the sign of that part for which ag is negative 
and still get a Lorentz-invariant function, (a9/|ao|)6(a"a,). It doesn’t 


(10.1) 
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really look Lorentz-invariant because you’ve got ag appearing in it, 
which refers to one particular Lorentz frame, but the Lorentz frame 
we take would not affect things, because a9/|do| is just a plus or minus 
sign, the same for any Lorentz frame of reference, because we don’t 
consider time reversal. 

We define the function A(a) as equal to (2a9/|do|) 6(a,a"). This 
A-function is quite an important function in the field theory of particles 
with zero rest mass and it is a Lorentz-invariant function. Let us study 
it a bit. We can write aa, = ag”? — |al|? = (ao — |al)(@o + |al). 
Now whenever we have a quantity like 6(xy), it is equal to 


|x\-* 8) + Lyl~ 8@d. 


That’s a general formula for the theory of distributions and we very 
easily check it by inserting both members as a factor in an integral 
and integrating things out. 

Let us use that result in our definition for A. We get 


ai 


M6) = ell ae Jal|* (ao — Jal) + Jao — [al |~* 5(a0 + Ial)]. 


Now in the first term we can put do equal to |a| and get 2|a| in the 


denominator; so the first term becomes |a|—! 5(a9 — |al). Similarly 
the second term is —|a|~* 5(a9 + |a|) and we have 
A(a) = (\a\)~* [8(ao — al) — dao + |al)}. (10.2) 


Expression (10.1) now becomes 


a PC) NOR ie. (10.3) 
4x? 
That proves that it is a Lorentz-invariant function. 
The function A is interesting and is worth further study. Id like 
to mention some of its properties. If we make a Fourier resolution, 
what do we get? I need not go through the analysis. The result is 


| A(k)e*? d*k = 437i A(a). (10.4) 


We get the same function again. We may interchange a and x here 
and get 


[aoe d‘a = 4n7i A(X). (10.5) 
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So A(k) is a function that transforms into itself under Fourier resolution, 
apart from a numerical coefficient. 
We can use that result in expression (10.3) and we get 


[Vix), Vx] = —A@ — x’). (10.6) 


That is the final result. It means that the field function V at two dif- 
ferent points in space-time commute except when the two points lie 
on the light cones of each other. When that is the case, we have a sort 
of 6-function singularity. 

I’d like to mention some other properties of the A-function. We 
have A(—a) = — A(a): It’s an odd function; it changes sign when 
you change the sign of all four a’s. If ag = 0, A(a) = O. It’s not clear 
from the definition just what A(a) is at the origin, but when you take 
into account its antisymmetric property you see that it has to vanish 
at the origin and everywhere on the three-dimensional surface ay = 0. 

Another property of the A-function that I want to deduce is the 
value of its derivative at ag = 0. Let us evaluate (0/daq) A(a) and 
then put a9 = 0. We can get the result most directly by using the 
Fourier resolution (10.4), 


dA(a) _ 
dag _ 


fan | koe A(k) dtk 


| k ik-a 
= Gab] ej OlKo — Nel) — 8(ko + Iki)e™* dk 


1 : 
ee i [3(Ko — |k|) + 6(Ko + |k|)e** d*k. 
Now let us put ag = 0; we get 


dA(a) 
0ag 


! . 
ao=0 ee hs | tote — |kl) + 6(kKo + lems Gop 


(k, x) is only a three-dimensional scalar product. kg now no longer 
occurs in the integrand except in the 6-functions, and therefore the 
integration with respect to kg is trivial. The result is 


gas 


ae 


= 4r 53(a). 
0 


ag= 


We can now understand the nature of the singularity in the A-function 
at the origin. The function itself vanishes at the origin and its derivative 
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with respect to ao is something which has the ordinary three-dimen- 
sional 6-function value when we consider it as spread over a), do, a3. 
I would like you to note that 


ae 


da,0a" 


A(a) = 9. (10.7) 


This follows at once from the Fourier resolution. 

If we look upon the a’s as coordinates of space, this A-function 
consists entirely of waves traveling with the velocity of light. These 
waves form a sort of pulse, an ingoing pulse to begin with, a pulse 
that goes to the origin, goes through the origin, and comes out again 
still as a pulse. That is how the A-function is to be pictured. Just at 
the origin, the function itself vanishes but its time derivative does not 
vanish. Those are the important properties of the A-function. 

You see why it comes into field theory in the equation for the com- 
mutator [V(x), V(x’)]. You obviously need to have for this commuta- 
tor something that can be resolved into Fourier components traveling 
entirely with the velocity of light, because (07/dx, dx“) V(x)'= 0, and 
if you apply the operator 07/dx" dx, to [V(x), V(x’)], the result is 


Q° 9° 
[V(x), Vix')] = i COAee))| = 0. 
OX, ax" OXp xe 
This checks with 
2 
g A(x — x’) = 0. 
@nenche 


That completes the proof that this scalar field theory is a relativistic 
theory even in the quantum theory. All our commutation relations are 
Lorentz-invariant. One could do a corresponding calculation for a 
nonzero rest mass. The results are essentially the same. We no longer 
have the A-function occurring; we have another expression involving 
a Bessel function, but it can likewise be proved to be Lorentz-invariant. 


$11. The Electron Field 


With the electron field we have four complex field functions 
wa (a = 1,2, 3,4). They may conveniently be considered as the four 
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elements of a column matrix y. They are assumed to satisfy the field 
equation 


in? + a") — man = 0, (11.1) 


where the a’s are 4 X 4 Hermitian matrices multiplied into the column 
matrix y by matrix multiplication. Any two of them anticommute 
and the square of each is unity. The field equations (11.1) are rela- 
tivistic when one makes the four components of y transform cor- 
rectly under Lorentz transformations. They have to transform ac- 
cording to the laws of a spinor. 

The conjugate complex field quantities Y. are considered as a row 
matrix y. Thus the matrix product Yay, for any of the a matrices, is 
a single element. The conjugate complex equation to (11.1) is 


in? + Par) + mhan, = 0. (UT) 


Let us first consider the field classically. Because the field equations 
are first degree in 0/dx9, it is very easy to put them in Hamiltonian 
form. One finds that one gets the correct field equations (11.1) and 
(11.2) if one takes 


H= | irda” +. mY 2emWz) d°x, (11.3) 
with the Poisson brackets 


[Wars Yo’) = 0, ware Woz] = dab 53(x a x): (11.4) 


Here yz is to be understood as y, at the point x,,x2,X3, for a 
particular time x9 which is understood, like the notation V, in Section 
9. The # occurring in (11.3) and (11.4) is just a parameter in the 
classical theory. 

It should be noted that the Hamiltonian (11.3) is real, because if 
one takes the conjugate complex integrand in (11.3) and makes an 
integration by parts, one gets back to the original expression. 

We have here the classical theory, and to apply quantization to it 
it is convenient to introduce Fourier components. In the case of the 
scalar field we made the Fourier resolution first four-dimensionally, 
and then we passed over from the four-dimensional to the three- 
dimensional Fourier resolution. Here we can make the three-dimen- 
sional Fourier resolution directly. We put 


Vn = jaie | ginny, dp, 
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with the conjugate equation- 


pe prst2 fete, ee 


This applies for each of the four components. , has four components 
like y, and is also understood as a column matrix. 

Now let us express the Hamiltonian in terms of the Fourier com- 
ponents. We get 


a i Vx(OrPr + ammye?", dx d°p 


, | tolaup, = mM Wp d° p. (ey 


We see that a,p, + a,m is the energy operator associated with each 
Fourier component p. 

Now there are four Fourier components ~p. a,;p; + amm has 
positive eigenvalues for two of them and negative eigenvalues for the 
other two. There’s a mixing up of the positive- and negative-energy 
Fourier components, and to carry on we have to separate them out 
and treat them differently. To do that we must transform the matrix 
Q;Pr + a&m,m into the diagonal form. This can be done by a suitable 
unitary transformation. 

We choose U, so that U,U, = 1 and such that U,(a,p, + anm)U, 
is a diagonal matrix. The diagonal elements are plus or minus the 
square root of p? + m? because (a;p, + amm)” = p® + m®, so 
(a,p; + amm) has the eigenvalues +(p? + m”)1/?. So the diagonal 
form of the matrix a,p, + amin 1S 


Gan) 0 0 4 
0 (p? + m?)}!2 0 0 
0 OSG ae 0 
0 0 0 —(p* He m?yil? 


Let us suppose we have found this unitary matrix U. There are many 
possible such matrices; it’s not uniquely determined at all, but any 
one of them will do for our purposes. Now we can write 


Uy (apr + amm)U, = (p? ae eee (11.6) 
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where V is the matrix 


lo @ 0 0 
al 0 0 
OO 0 
0 0 OF 


We now express H in the form 


H = i, Op U,(arPr + amm) Uy Upp dp. 


I’ve merely introduced U,U,, which is equal to 1. We get 


Te [ote SE mn). (VO ae 


Now let us put Upp = ¢p- $p also has four components. We get 


y= [oe = m’)'!?6 Vey ay 


There we’ve done the transformation which separates out the positive- 
and negative-energy parts of the field function y. If we write H out 
in full it reads 


H= fo oe m?)"!?(b1pO1p ae Poph2p coat $3pP3p oe P4ph4p) d° p. 
able) 


We might quantize this H according to Bose quantization; the result 
would be a theory in which some of the particles have negative energy. 
The Hamiltonian here would contain some negative-energy terms, 
which would result in a nonphysical theory with negative-energy 
particles, which we don’t want. Alternatively we might use Bose 
quantization and change the Hamiltonian, putting plus signs in the 
third and fourth terms of (11.7). That would give us a theory with 
only positive-energy particles, but then we should have destroyed the 
localized expression for the Hamiltonian. 

We cannot apply Bose quantization to this theory without either 
bringing in negative-energy particles or spoiling localization. We can, 


11. The Electron Field 57 


however, apply Fermi quantization while keeping the localization 


and keeping only positive-energy particles. That is why we use Fermi 
quantization for this kind of field. 


’ With Fermi quantization we may take 


H= [oe =e im)" (Bisbip ale P2ph2p ae $3pO3p Gi b4pP4p) dp. 
(ules) 


This differs from (11.7) by an infinite constant. We have interchanged 
the order of ¢3, and ¢3, and of ¢4, and ¢4,, which involves bringing 
in a minus sign with Fermi quantization and also involves a constant 
to be added on. So this new Hamiltonian is equal to the preceding 
one apart from an infinite c-number. An infinite c-number, of course, 
doesn’t invalidate the use of this Hamiltonian in the Heisenberg 
equations of motion. It will only come into play when we try to 
interpret the Hamiltonian as an energy. 

Now if we’re going to use (11.8) as our Hamiltonian and are going 
to have only positive-energy particles, each of the operators that occurs 
to the right of a product must be an absorption operator in order that 
it give zero when applied to the vacuum ket. We then have a vacuum 
ket |0) with H|0) = 0 by taking ¢1pl0) = $2p|0) = $3p|0) = 
$4p|0) = 0. That’s how we have to define our vacuum ket in this case. 
It’s the only way to do it if we are to have all our particles with a 
positive energy. It involves counting ¢1p, d2p, d3p, G4p aS absorption 
operators and $1y, dap, o3p, dap aS Emission operators. As far as 
the ¢, and ¢2 components of ¢ are concerned, we just have the usual 
kind of second quantization referring to the emission and absorption 
of the positive-energy electrons, which appear as ordinary electrons. 
As far as the 3 and the 4 components are concerned, we have inter- 
changed the operations of emission and absorption that a direct 
quantization procedure would give us. This means we are working 
with a picture where the real physical things are just holes in what we 
would get with a direct method of second quantization. These holes 
are the positrons. 

We have to interpret $3, and ¢4, as emission of a positron with 
momentum —p. In the primitive method of quantizing, where ¢3p 
and $4, are absorption operators of electrons in negative-energy 
states, these operators imply the absorption of a quantity of momentum 
p, and if we are going to reinterpret them as emission operators they 
have to be the emission operators of a particle with momentum —p. 


58 Lectures on Quantum Field Theory 


Of course, $1, and 2p are the emission operators of an electron with 


momentum: +p. 
The classical Poisson bracket relations (11.4) get replaced by the 
following with Fermi quantization: 


[pan > bop + a dapPbp' ala Pop Pap = 0 
and 


[bans Pop']+ = Sap 53(p — P’)- 


These are the basic anticommutation relations coming into the Fermi 
quantization of the electron field. 

How does it go if we transform back to the y’s? Let’s work that 
out. We have evidently [Wap, Wop']4 = 0. Further, 


aes Wop ly = [i patiee ’ bap Uprar]s - 


The matrix elements of U, and U, are just numbers. They can be taken 
out of the anticommutator, so that we get 


ead Wop'l+ — OU eUsray bed 53(p 2 p) 
Sab 53(p a DP); 


where we have used 
Coee Ces = bab: 


So the anticommutation relations referred to the y’s are just the same 
as referred to the ¢’s. This kind of unitary transformation doesn’t 
disturb the anticommutation relations. 

We can now transform a stage further and go back to the original 
y,’s, and we find 


[Woz » Vox] = 0 and [Wars Wor] = bab 53(x ae x), 


The Fourier transformation back from the p’s to the x’s doesn’t affect 
the anticommutation relations. 

The Hamiltonian is the same as the original Hamiltonian (11.3) 
which we had in the classical theory, except for an infinite constant. 
Of course, this infinite constant just represents the energy of the sea 
of negative-energy electrons. 

I want now to get a formula for the total charge. The charge is made 
up from a contribution of —e for each electron and +e for each 
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positron: 


Q <e|liiey, + dopbo» — b3pdap — b4pb4p] dp 


aes + dopdop + Fanb3ap + Papp] d?p + inf. const. 


This infinite constant corresponds to the total charge of all the negative- 
energy electrons. Now the above expression, neglecting the infinite 
constant, becomes 


—e| Yate dp a =e | tut ax 


So we can say that we have a charge density jp) = —ey,,, neglecting 
an infinite constant charge density. We thus have both a localized 
energy density and a localized charge density, neglecting infinite 
constants. 

We know how the y’s transform under Lorentz transformations, 
from the one-electron theory. In the present theory, with second 
quantization, they transform in the same way. If we now make Lorentz 
transformations, we find that the charge density is associated with a 
current density 


Jr = meer oe 


The charge density and current density can be combined into one 
formula 


We — —eprays, CUE) 


where ag is understood to be the unit matrix. We have dropped the 
infinite constant added to jo, which corresponds to the charge density 
of all the negative-energy electrons. 


§12. The Electron Field: Relativistic Invariance 


Now, starting from classical equations and using the rules of quan- 
tizing according to the Fermi statistics, we obtained the quantum 
theory of Section 11. The question remains whether this theory is 
really relativistic. The procedure of quantization involves working 
with one particular Lorentz frame of reference; it is not a procedure 
that is manifestly Lorentz-invariant, and it becomes necessary to 
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check the Lorentz invariance of the quantum theory by a direct cal- 
culation in the quantum theory. 

We can do that by a method analogous to what we did for the 
scalar field quantized by the method of Bose statistics. We take two 
field quantities at two general points of space-time and see what their 
anticommutator is. We evidently have 


[Wa(x), Pox’) J+ = 0 


(x means here a point in space-time x = X90, X1, X2, X3, in contrast 
with the suffix x, which refers to just a point in three-dimensionsl space 
at a time that is understood). Let us put [Wa(x), Yo(x’)]4 = Kav(x, x’). 
Now, we could evaluate this commutator by expressing the two field 
functions in terms of their Fourier components. We know the anti- 
commutation relations between the Fourier coefficients, so we have 
enough material to do the calculation. We can get the result more 
directly by noting that this quantity K has certain properties which 
are sufficient to define it and that there is a pretty obvious answer to 
what K must be to have these properties. 
The properties of K are: 


(1) Kao(x, x’) involves the x’s and x’’s only through (x, — x,-). It 
is obvious that it is only the differences that can come in, because the 
whole thing is invariant under translations in space. 

(2) It satisfies the wave equation 


[o,ih(d/AXy) — amM]anKy.(x — x’) = 0. 


The operator [a,ih(d/dX,) — o&mM]qp acts only on the ¥;(x), whether it 
occurs on the left or on the right of Y,(x), and applied to y,(x) it gives 0, 
because ¥(x) satisfies the wave equation (11.1), so K(x — x’) satisfies 
the wave equation, too. 

(3) For x9 = x0’, Kae(x, x’) = 63(x — x’) 5q4. This third equation 
is just the anticommutation relation for the field functions at two 
points at the same time, which we had previously. 


K satisfies those three conditions and conversely those three con- 
ditions are sufficient to fix K. We can see that in this way (3) fixes K 
when x9 = Xo’. Now if you keep xo’ fixed and allow xq to vary, the 
wave equation tells you how K varies. So now we have K determined 
for all xo for the given x9’. Then condition (1) shows you how K 
varies when x9’ is allowed to vary. These three conditions thus fix K 
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completely, and if we can write down a K that satisfies those three 
conditions we know that it is the correct one. 


Well, there is a K which pretty obviously satisfies those three 
conditions: 


Kao(x, x‘) = h-* / DS ; (1 ee ~ curt) CaN * dp. 
(12.1) 


Here po means -k(p? + m?)!/? and the summation refers to the sum 
over the two values for pg. Now let us check to see if this K does satisfy 
the three conditions. Well, (1) is obvious; (2) is also fairly obvious, 
because if we apply the differential operator to K it will bring into the 
integrand on the left a factor (po — a,;p, — am,m), which, multiplied 
into [1 + (ap; + amm)/po], gives zero. Let us see whether (3) is 
satisfied. Put x9 equal to x9’. That means that the term involving 
Po in the exponential disappears, and we get 


Sash f eite—s'ni ia = Opler —= 2). 


Therefore all three conditions are fulfilled and this K (12.1) is the 
value of the anticommutator. 
Now we want to show that K is Lorentz-invariant. We can write it 


IRC) = da 42 GEIL a= Ga (12.2) 
where 


Jo = prs f ete an 
W(ar—ax’)- d® 
J, = aaD? pen pin oe. COLS) 


3 
= —1(x—2x’)- 4d p 
= [> me ‘ 2 Wir 


Let us see how these three integrals transform under Lorentz trans- 
formations. 
Well, let’s take J,,. We may write it 


Jn = n-8 f mente 98 Ny = = ae 


where po is now an independent variable. To see that both expressions 
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for J, are equivalent we transform the 6-function, 


3(po” — p? — m”) 

= 8lpy ep? em) po eee a 
[pee @ 2) |" Ope el 
+ |po — (p? + m?)¥?|-! apo + (p? + m?)"?] 
[2(p? + m)"?)-! {8[po — (P? + oy] 
+ apo + (p? + m?)"7J}. 


So the 6-function gives rise to two contributions according to whether 
po is +(p? + m)"? or —(p? + m’”)"!?; and that gives you just the 
two contributions occurring in (12.3). Well, with J,, in this form it is 
obvious that it is Lorentz-invariant; all four p’s are now treated on the 
same footing, so J», equals an invariant; we do not need to work it out. 

Now, the other J’s can be treated in a corresponding way. We can 
combine them into a single formula for J, : 


J, = wef pye 2 ?* 5(p2 — p? — m?) 2% dp. 

|Pol 
This we see equals a four-vector. Well, that tells us how these co- 
efficients J are affected by a Lorentz transformation. 

You might, when you look at expression (12.2) for K, think that it 
is wrong, because, with ag = 1, it is K = aopJo + an, + OnJm, 
and you might be inclined to think that there should be —a,J, to 
give a Lorentz-invariant quantity. That plus sign is quite correct, 
though. We may write 


K = OQyOmdI” + otmdin - (12.4) 


You see that this gives the correct signs, because, putting u = 1, 2, 3, 
we have in a, a quantity that anticommutes with a, and gives us the 
minus to compensate for the minus that occurs when we raise the 
suffix. 

Let us apply a Lorentz transformation to K given by (12.4) and see 
what happens toit. We can make a Lorentz transformation that changes 
vectors according to A,* = a,’A,, where the coefficients satisfy the 
equation a,"a"? = g’?, The inverse transformation reads A, = a’,A,*. 
One easily checks that these relationships are consistent with each 
other and preserve the squared length of the vector. 

Now suppose that under this Lorentz transformation y transforms 
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into y* = Ly, where L is some 4 X 4 matrix. We then have Y* = JT. 
Now Kz,» transforms like ya¥,, so 


K* = LKL = Lend ym LJ" + LoinlJin - (12.5) 


From one electron theory we know that the quantities Ja, transform 
like a four-vector. This in fact is the charge-current density four- 
vector. That means that ~*a,y* = a,’Pa,y. The left-hand side here 
is equal to ~La,Ly. This result holds for any wave function y and 
we get, therefore, 


bated, = Ay a,. (12.6) 
Now we also know that a,,W is an invariant and is equal to 
Vr amy* = PLomLy, 
Se 
Dak = Ohya (12.7) 


Equations (12.6) and (12.7) have Z on the left and L on the right, 
and we want to have some equations in which it is the other way 
around, L on the left and Z on the right, to evaluate (12.5)... Multiply 
Gea abyela som the leit; we gel La,,la,t — £. We can cancel J 
on the right on both sides, because it is a matrix without vanishing 
determinant. Then we can transfer a,, to the other side and get 


[etnil, = Che. (12.8) 
Multiply (12.6) by La, on the left and am on the right: 


Lotm Lo, LotmL = Oty’ LotmOtyOmL. 


Using (12.8) we get 


Cm Oy = Ay? Lotyy Oty Om L- (12.9) 
Now we substitute for J” in terms of J#* in the expression (12.5): 
K*® = LamayomLay’I* + Lom ESn*. 
Using equations (12.9) and (12.8) we get 
K*® = ctmOy0tml**® 4+ oimdm®. 


There you see we have K* of the same form as K, and that proves 
the Lorentz invariance of the anticommutation relations. The Lorentz 
invariance of the whole theory is thus established. 
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§13. Fields with Interaction: Relativistic Invariance 


As a result of this work we have a satisfactory quantum theory for a 
scalar field, a field of bosons, or the electron field with Fermi quan- 
tization. But we have a satisfactory quantum theory for these fields 
only by themselves. To get something interesting, we have to put in an 
interaction. 

The method of putting in an interaction is to add up the various 
Hamiltonians for the various fields alone and then add a further term 
to the Hamiltonian to correspond to the interaction energy. The 
question is: What shall we take for the interaction energy? Well, we 
can work from the classical theory, which will suggest an interaction 
energy for us to take. We can take over this interaction energy from 
the classical theory, and then we have a quantum theory. We then 
have to check, from calculations performed entirely in quantum theory, 
whether it is all right. In particular, we have to check whether our 
theory is really a Lorentz-invariant theory with the interaction. To 
do that we must work out the field equations, given by the Heisenberg 
equations of motion, and check whether they are Lorentz-invariant. 
Then we have to check that our commutation or anticommutation 
relations are Lorentz-invariant. 

I want now to discuss the question of how we deal with the Lorentz 
invariance of the commutation or the anticommutation relations when 
there is interaction. The work in Sections 10 and 12 shows that it is 
all right when there is no interaction; we have to generalize this work 
to take interaction into account. 

Well, you might say, let us use the same method when there is inter- 
action. That is going to be extremely difficult, because we have in the 
commutator or anticommutator two different times; and if we want 
to compare quantities referring to two different times we have to 
integrate the Heisenberg equations of motion. To integrate the Heisen- 
berg equations of motion when there is interaction is a very com- 
plicated process. We need some other method. 

Instead of taking two arbitrary points in space-time, which is really 
more than necessary, let us just start from the commutation relations 
connected with field quantities at one particular time, that is, field 
quantities on a three-dimensional surface x9 = constant in space- 
time, and let us apply an infinitesimal Lorentz transformation, so as 
to go over to another three-dimensional surface inclined with respect 
to the first one. Let us then see what are the commutation or anti- 
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commutation relations for the field quantities on the new  three- 
dimensional surface. If we can check that they are all right, we should 
know that our commutation or anticommutation relations remain in- 
variant under infinitesimal Lorentz transformations; and if they are 
invariant under infinitesimal Lorentz transformations, they will be 
invariant under finite ones. 

We start off with our various field quantities on a surface xy) = con- 
stant in space-time and we make an infinitesimal Lorentz transforma- 
tion to the surface x9* = constant as in this space-time diagram: 


X,* = Constant 
Infinitesimal Lorentz transformation 


Xo = Constant 


We have known commutation or anticommutation relations between 
the quantities for a given x), which are assumed not to be affected by 
the interaction. We shall calculate from them the commutation or 
anticommutation relations at x)* = constant. 

The Lorentz transformation will involve a change in the time co- 
ordinate x9* = Xo + €V,x,, where e€ is infinitesimal and eV, gives 
the relative velocity of the two frames. The V,’s are c-numbers. We 
shall neglect e”. Now we take two field quantities K and \; they may 
be U or V or or Wag. We have 


ae ate caee Gs" = x0) Fe = iC GT aeawall 


where H is the total Hamiltonian including interaction. With \ we 
do the same thing: 


Az (Xo*) = Az’(Xo) + EV 3x5 [Nears Ai}. 


Now let us work out the Poisson bracket for those two quantities. 


Recenter 1K Ca), ha eo) eV GIL Ky 721, X27] 
AP eV X5'[Kz, [Ax 3 veal 
= [K,(xo), Az’(Xo)] += EV Og! == Bee) Gay [Az’s veal 


Bee Koa (13.1) 


where the last line was obtained with the help of Jacobi’s identity. 
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Now there are some occasions when we are interested not in the 
Poisson bracket of two field quantities but in their anticommutator; 
if they are both y’s or one of them is y and the other . So we have 
to do the calculation also for the anticommutator. We have 


[Kz(x0"), Nx(X0)" 1+. = [Ke(X0), Ae (Xo) 4 

= eV xllKe, A], Xe a ee [ee ee 
In this formula you note we have sometimes the anticommutator and 
sometimes the Poisson bracket; where we have the Hamiltonian coming 
in, we have to keep to the Poisson bracket. The formula involves a 
mixture of anticommutators and Poisson brackets, but we can handle 
it as before. 


EGER MACRO = NEE ere. 
+ V(x, — x;,Ke, Dor Ale 
+ &V,x-[Kz, (or, Aly + ([Ke, A214}. 
(iB) 


Using the notation of the footnote on page 37, we have the formula 
lie: bl ele lao e lesa lene) oeralie (13.3) 


You can check it just by multiplying it out with straightforward non- 
commutative algebra; you find that everything fits. We have here a 
formula for anticommutators analogous to the Jacobi identity. With 
the help of it we can write 


{[Kz, Dwr, HI, + ([Kz, A], X21} = [[Kz, dv], A 


The expression (13.2) then becomes similar to (13.1). It becomes 
possible to combine them into a single formula: 


[K2(x0"), \x(X0 My. = [Kx(x0), Ax"(Xo) 4 
ar eV,(x,’ iar x,)[Kz, [Az ’ AY],,, 
See Albscen Molle, Jeb (13.4) 


This one formula sums up the two results for the Poisson bracket or 
commutator and for the anticommutator. 


K and } so far have been quite arbitrary field quantities. I am now 
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going to take them to be basic field variables. Take K, to be one of 
i ey eva andi Aeato be one of V.-, U,V... vase. Winder those 
conditions [Kz, dz], = 0 or +63(x — x’). In fact, all the basic field 
quantities we have been concerned with satisfy commutation or anti- 
commutation relations of this type. We either have naught on the 
right side or we have the 6-function or minus the 6-function. Thus 
[Az, Az]; = c-number and [[Kz, \x]i, H] = 0. If we took more 
complicated local things, e.g., if we took K and } to be functions of U 
and V, that wouldn’t hold, of course. But we are only interested in 
the commutation relations of the basic field quantities. 

The last term in (13.4) has thus been cut out. Now let’s see what 
we are left with: 


eG ee Gy eG eal 
Si EIAs) = SI | Nas Je celle 
Sie eV, ee IK, Dee, Hells, 


where I have taken into account that H is equal to the proper energy 
Hp plus an interaction energy Hg. Proper energy means the energy 
of the fields by themselves, without any interaction; and Hg is the 
interaction energy which is added to Hp. 

We now see the condition for the commutation or anticommutation 
relations to remain Lorentz-invariant. We know that they are Lorentz- 
invariant when there is no interaction; that we proved previously. 
So we want the interaction not to disturb the commutation or anti- 
commutation relations, and that means that 


(x,’ co Xr )Kz , [Az , Heal], = 0, (eke)) 
[eo oll multiple of 6;(@¢ — x’). 


Take Hg = {Woz d®x, where Wo, is the interaction energy density. 
I am assuming a local theory in which the interaction is given by a 
localized energy density, integrated over three-dimensional space. We 
have now to choose this function Wg so as to satisfy the condition 
(13.5). We see that the condition is satisfied if Wg, equals a function 
of U,, Vz, Wax, Vax With no spatial derivatives. In general, of course, 
you might have an interaction energy density involving spatial deriva- 
tives, derivatives with respect to x,, X2, x3; but suppose we take the 
case when this is not so. In that case you see that the Poisson bracket 
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of Hg with \,, will involve the field quantities at x’ without any deriva- 
tives, and the commutator or anticommutator of this Poisson bracket 
with K, will give some multiple of the 6-function without derivatives 
of the 6-function. 

So here is quite a substantial class of interaction-energy densities 
we can take, which we can be confident will not destroy the Lorentz- 
invariance of the commutation or anticommutation relations. This 
class is sufficient for practical purposes. Theoretically one can 
generalize a little bit but that generalization hasn’t turned out to be 
useful. 

That is a formal proof of Lorentz invariance. It could be criticized 
on the grounds that we are not making proper use of the 6-function 
in this kind of work, because what we are doing is taking two field 
points close together and expanding the commutator or anticommutator 
in a Taylor series in (x — x’). If these two field points are very close 
together, [K,, \z’],. cannot in general be expanded as a Taylor series 
in (x — x’). It is obvious that it cannot, because something very 
singular happens when one of the two points x and x’ is on the light 
cone of the other. 

The above use of the 6-function is not a correct one, but we can put 
the calculation in order in a rather simple way. We form fa.K» ak 
and {bz/\,- d°x,, where a and 6 are continuous c-number functions of 
the three x’s, and then we form [fa,K, d?x, [b,x d°x’]. That is 
something that does not involve any 6-functions. It just involves inte- 
gration over a and 5, and it is quite permissible to calculate how that 
quantity transforms when we make a Lorentz transformation, using 
a Taylor expansion to connect the things on x,* with the things on 
X9. We can then go through work that would be exactly the same as 
the work above with all our equations multiplied by a,, b,- and inte- 
grated over x or x’. Then, of course, there wouldn’t be any 6-functions 
coming in at all. 

You can put the calculation right in that way. It would be rather 
tedious to write out all these equations with these extra factors coming 
in—the equations are quite long as it is; but you see that it will be quite 
possible, and then we should not be making any unjustifiable Taylor 
expansions. So although the calculation as I gave it is not really 
justifiable as it stands, it can be put in order in a fairly easy way. One 
gets the same conclusion that, provided the interaction energy satisfies 
the above-mentioned condition, the commutation or anticommutation 
relations are Lorentz-invariant. 
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§14. Constraints 


I now want to go on to another subject, the question of constraints 
in the quantum theory. When we work from a classical theory and 
put it in the Hamiltonian form, we might very well get constraints 
appearing—equations involving the q’s and p’s only. These constraints 
are of two kinds: first class and second class. 

The second-class constraints can be eliminated. With the help of a 
suitable modification in the definition of Poisson brackets, the second- 
class constraint equations can be taken to be strong equations or 
definitions. We are left with first-class constraints. We have now to 
consider the general problem of a Hamiltonian with first-class con- 
straints. How are we to take that over into the quantum theory? 

A full account of classical constraints is given in my book, referred 
to on page 42. The quantum theory has to be reformulated in 
terms of the Heisenberg picture, without the help of the Schrédinger 
picture. 

The constraint equations are weak equations, written ¢,, ~ 0, which 
means they must be used in accordance with certain rules. I'll first 
give the classical rules. 

Constraints can be added. They can be multiplied by factors: the 
factors can be functions of the q’s and p’s. We must not use them in 
Poisson brackets, 1.e., [¢m, /] ~ 0 in general. 

Now we have to set up suitable quantum rules. We maintain the 
addition rule. Now what about multiplying? If we can multiply them 
by factors on both sides, we can form a quantum Poisson bracket; 
and that we must not do, as we are not allowed to form Poisson 
brackets in the classical theory. How can we avoid the possibility of 
forming Poisson brackets? We do that by adopting the rule that we 
can multiply the weak equations by factors only on one side, for 
example, the left side. So we take as our second rule that we can 
multiply by factors, functions of the q’s and p’s, on the left and that 
we cannot multiply by factors on the right in general. So we cannot 
form Poisson brackets. 

These are the rules we take to be the rules for working with weak 
equations or constraints in the quantum theory. These rules destroy 
the symmetry between left and right multiplication. This symmetry 
between left and right, which we have had all the time in our quantum 
theory up to the present, gets destroyed when we have a quantum 
theory with constraints in it. 
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When we have weak equations in the classical theory we need certain 
consistency conditions between them. We have the classical con- 
sistency conditions 


These consistency conditions we assume to be the same in the quantum 
theory. 

In the classical theory the consistency conditions may be written 
as strong equations, 


[Om > on] a Amnj9j - (14.2) 


Anything that is weakly equal to zero has to be strongly equal to some 
linear combination of the ¢’s, because the ¢’s are the only independent 
quantities weakly equal to zero. In the quantum theory the consistency 
conditions may also be written as strong equations like (14.2), but now 
the coefficients have to be on the left. In the classical theory it doesn’t 
matter whether these coefficients a,,,,; are on the left or right, but when 
we go over to the quantum theory we must insist that the coefficients 
are on the left. 

In the classical theory we have the result that the Hamiltonian is 
first class. This means that 


[om, H] = 0. (14.3) 


In the quantum theory we maintain this condition. These weak 
equations can be written as strong equations: 


and in the quantum theory we must again insist that the coefficients 
shall be on the left. 

Those are the rules by which one handles constraints in the quantum 
theory. Now, you may be rather perturbed by our giving up the sym- 
metry between left and right multiplication; you will say, does it 
matter that we have to give up this symmetry? It matters in one 
respect, with regard to our definition of real g-numbers. We shall have 
to go into this question of reality again. 

Before we do that, I should like to develop the theory further. Let 
iE use the constraints ¢, ~ 0 to impose conditions on the ket vector 
les 

on|2) — (14.5) 
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The question arises: Can we find a ket to satisfy all those conditions? 
Let us work on those conditions a little and see whether we get an 
inconsistency. From ¢,,|P) = 0 and ¢,\/P) = 0 we can deduce 
bndm|P) = 0 and dnon|P) = 0; therefore [6,, ¢m]|/P) = 0. This 
condition is not independent of the previous ones. From (14.2) it is 
just Gnmj;|P) = 0. So working with the restrictions (14.5) on the 
ket in this way we cannot get new restrictions, and under those con- 
ditions one can assume that the ket |P) exists. 

With our rules for working with the constraints we can thus use the 
constraints to define kets by (14.5), but we cannot use them to define 
bras; it just wouldn’t work out, because the factors would be on the 
wrong side. A ket that satisfies these conditions is called a physical 
ket. That’s a definition. 

We now define a g-number to be physical if, when it is multiplied 
into any physical ket, it gives another physical ket. The physical kets 
form a subspace in all the ket vector space and the physical g-numbers 
are operators that operate in this subspace. If a is a physical g-number 
we have for any ¢,, dna|P) = 0, and hence (dna — adn)|P) = 0 for 
all physical kets |P); so 


[dns a] = 0. (14.6) 


Conversely, the conditions (14.6) ensure that a is a physical g-number. 

Perhaps I should say at this point, to give you some physical insight 
into all this, that we shall later be applying this theory to electro- 
dynamics and we shall then find that the quality of being physical is 
just the quality of being gauge-invariant. The physical g-numbers cor- 
respond to quantities that are gauge-invariant and the physical kets 
are kets that are gauge-invariant. 

You see that if a and 8 are physical so is a + 6, and of course af 
and Ba are also physical. The Poisson bracket of a, 8 is also physical. 
So the physical operators form a ring. You can add and multiply 
them freely and get other physical operators from them. 

Equations (14.3) shows that the Hamiltonian H is physical. We now 
see that, if a is physical, sois &@ = [a, H]. Thus if ais initially physical, 
it always remains physical. 

Now let us consider the question of conjugate complex q-numbers. 
At this stage we see that we can define the conjugate complex of a 
g-number if it is a physical g-number. We can define it in the way we 
did in Section 3 for the case of no constraints. We assume that every 
physical a has a complex conjugate @ that is also physical, and we 
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assume that 


atB=-a+f and of = Ba. (14.7) 


Each ¢, itself is physical, according to the definition. It is to be 
counted as real, because it has the value zero for all physical kets. 

These assumptions for conjugate complex g-numbers cannot be 
extended in general to nonphysical g-numbers. One can see this as 
follows. We saw that ¢, is real. The equation expressing this should 
be written as a weak equation, 


Similarly, for any g-number a, am is weakly real, 


abm ~ abm- 
Take a to be real and apply the second of the rules (14.7). We get 


adm ~~ Pme. 


For a general g-number a@ this conclusion would be false, but for a 
physical g-number a it is all right, according to (14.6). So we have a 
consistent scheme only as long as we restrict our ideas of conjugate 
complex g-numbers to the physical g-numbers. That is really sufficient 
for applications of the theory. 

The importance of a definition for reality is that it will enable us to 
deduce that real physical g-numbers have real eigenvalues, and that is 
most important when we connect our theory with measurements. The 
fact that we don’t have a definition of reality for the nonphysical 
dynamical variables doesn’t matter, because they are not quantities 
which can be measured. Well, this is the general situation with regard 
to constraints in the quantum theory and also with regard to conjugate 
complex g-numbers. 

There is a special case where we can introduce the notion of the 
conjugate complex for nonphysical g-numbers. That is the special 
case which occurs when the coefficients @mnjs Uma, in (14.2) and 
(14.4) commute with all the ¢’s, for example, if they are c-numbers. 
Then 

[om, Gn] = AnnjPj = OjAmnj, 
lm, 1] = bmnodn = bnbinn - 


If we have this special case, then in our basic equations there is no 
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dissymmetry between left and right multiplication and we could con- 
tinue to develop the theory symmetrically between left and right 
multiplication. We can then maintain the definition for conjugate 
complex q-numbers even when the q-numbers are not physical. So 
we have this special case in which it is possible to preserve the notions 
of conjugate complex and reality for nonphysical g-numbers. 

It so happens that this special case is realized in the case of electro- 
dynamics. Electrodynamics is one of the important applications of 
constraints. Another important application is to the gravitational 
field or to any theory involving curvilinear coordinates. We have then 
constraints to express the conditions that our equations must remain 
invariant under transformations of coordinates. In case of these con- 
straints associated with curvilinear coordinates, we do not have the 
special case, and then we do not have the possibility of defining real 
g-numbers among the nonphysical q-numbers. We can define real q- 
numbers then only for the physical ones, which would be the ones 
independent of the system of coordinates. 
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I should like now to proceed to develop electrodynamics. Let us 
first take the electromagnetic field by itself. We can use the classical 
theory as a basis for getting a quantum theory. There are different 
forms in which the classical theory can be put. I propose to work now 
with Fermi’s form, in which we use the constraint 


Ody 


oe (15.1) 


This is the form that is mostly used in practice; it has the advantage 
that there’s rather more symmetry between the four coordinates of 
space-time than with other forms. 
The classical action density & for this form of electrodynamics is 
given by 
A eee AG, (15.2) 


where F,, = Ay,» — Ay». Well let us work out what we get with it. 
First, we see that it can be transformed into 


—1F,,FY” + $A,A, + perf. diff. 


74 Lectures on Quantum Field Theory 
Now we can make a further change to 
4A, (A? — A) — BAP AS. 
plus another perfect differential, and we are left with 
4xr& = —$A,,A*” + perf. diff. (s2)) 


The per-ect differential does not affect the field equations and can be 
neglected. 

Our action density is now such that for each value of yu it is just the 
action density (9.1) for a scalar field with no mass term, except that 
there is a different sign for » = 0. We have here the action density 
for four scalar fields. So far as the action density is concerned, and 
ignoring the constraints (15.1), we have four independent scalar fields ; 
they just happen to be labeled uw = 0, 1, 2, 3, but that labeling of the 
four fields as the four components of a four-vector is not relevant to 
the derivation of equations of motion from the action principle. You 
see at once that the equations of motion are just [_]A, = 0, like (9.2). 
It is obvious that these equations of motion are consistent with the 
supplementary condition (15.1). So what we do is work with the 
action density (15.2) or (15.3) and bring in the Lorentz equation 
(15.1) as an extra condition which we impose on our dynamical 
variables; and we see that this extra condition is consistent with the 
equations of motion. 

Let us now proceed to apply the usual method for going over to 
the Hamiltonian. We introduce the momentum variables —47B* 
conjugate to A,, as with the scalar field in Section 9. We put the 
minus sign here so as to have 47B, conjugate to A, for the case of 
uw = 1,2,3. The B, then have the values B, = A,'° and satisfy the 
Poisson-bracket relations 


[Aye, Byer] = —4rgy,53(x — x’). (15.4) 


We now get for the Hamiltonian (call it Hp), 


Ae = ne | (B,BY + A,"AB) dx, (15.5) 


which is just the Hamiltonian (9.3) for four scalar fields with a minus 
sign for the one with u» = 0. The degree of freedom associated with 
Ao, Bo give a negative contribution to the energy. If we take p = 1, 2, 
or 3, we get the usual positive contribution. 
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So far this is just a treatment of four scalar fields. We now have to 
bring in the restriction on the potentials (15.1). Let us see what con- 
ditions they give us on the Hamiltonian variables at a certain time. 
These conditions will be the Hamiltonian constraints. At one par- 
ticular time we must take 


AG =) (15.6) 
and 


d HN 
a (15.7) 


These will be independent if we think of them just as conditions holding 
at one particular time. So there we have two independent conditions 
which we are imposing on the A’s. They are all the independent ones, 
because the further condition (07/dx90x9)A,"* = 0 is a consequence 
of the previous ones and the field equations [ JA, = 0. This further 
condition and all the further ones we get by successive differentiations 
with respect to time are not independent but are consequences of the 
equations of motion. 
The condition (15.6), expressed in terms of Hamiltonian variables, 
gives 
Bot A,’ = 0. (15.8) 


The condition (15.7) may be written with the help of the field equation 
[ ]Ag = 0 as Ag’ + A,’ = 0. This is just the Maxwell equation 
div § = 0. Expressed in terms of the Hamiltonian variables it reads 


Bee 0, (15.9) 


We have in (15.8) and (15.9) all the Hamiltonian constraints. Each 
of these equations represents a threefold infinity of constraints, one for 
each point x;, X2, X3 in three-dimensional space. We may write them 


1 a 0, 2 a CF 
where ¢, and @2 are defined by 
¢, = Bo + A,”,7 os = By AG’. (15.10) 


To check the consistency conditions (14.1), we have to check that 


[Piz Piz’) * 0, [Pics ox] = YO, [por 22'] = 0. dail) 
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These come out very directly. In fact, they hold as strong equations, 
with all the coefficients a,,,; of (14.2) vanishing. 

We quantize that theory by treating the four A’s like four scalar 
fields and using the Bose statistics. The four fields are quite detached 
as long as one does not use the constraints. The constraints are the 
only equations which give any relationship between them. If we forget 
about the constraints for the time being, we can quantize the theory 
just like four scalar fields. 

We introduce Fourier coefficients A,;, just as we did before in the 
case of the scalar field. We take the classical expression (15.5) for the 
Hamiltonian into the quantum theory, and then if we put it in terms 
of the Fourier coefficients we get 


Hp = = 40? | kA A d®,, + infinite c-number, 


corresponding to (9.16). The Fourier transformation gives directly 
the symmetrical product $(44 + AA), and to replace it by the un- 
symmetrical one with the 4’s all on the right involves the introduction 
of an infinite c-number. 

The point to be noted about this Hamiltonian is that the contribution 
for the zero degree of freedom is negative. I’d like to write that out 
explicitly: 


Hp = an? [ eAnA + Aor Aon + AgrAzsx — AorAox) ak 
+ infinite constant. (loa12) 


We have this rather peculiar thing: a negative energy associated with 
the Ao waves. 

The commutation relations between the Fourier coefficients are, 
corresponding to (9.12) and (9.13), 


Annee ui = (0). Ape doe = —!8w es y 
[Ants Aree] [Aue Avil 42 \k| 53(k — k’) (15.13) 


They lead to the commutation relation between the 4’s at two general 
points in space-time: 


[Au(%), 4,(X)] = Bu Ae — x’), (15.14) 


corresponding to (10.6). 
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We get the above results simply from quantizing the four scalar 
fields, and now we have to bring in the constraint equations. We 
must first check that they are consistent with each other. We already 
did that checking in our classical work, and it goes the same in the 
quantum theory. Alternatively, we may check the consistency directly 
from (15.14). We can do that by forming 


feasco za) = rey ey 


Ox, Ox,’ Ox, OX, 
a2 
= Ax, ax,’ Eu A(x — x) 
a2 
es Ox Ox a) 


from (10.7). 

The various Poisson brackets (15.10) all vanish strongly. That’s 
really more than is necessary. We need them to vanish only weakly 
and find them to vanish strongly, with all the coefficients a,j; zero. 
So we have the special case where it is possible to preserve the ideas of 
conjugate complex g-numbers also for the nonphysical ones. 

Well let us find out now which are the physical dynamical variables. 
A physical dynamical variable has to commute with all these 0A,/dx, 
quantities; at any rate it must commute weakly with them. Suppose 
we form 


aA” AAP ue, 
On Ox, Oboe 


Well we can work that out as 


ae ; ; 92 ‘ 
9x, Ox, [A'), Aol — 55 9,7 IA (x), A(x’) 
a” 2 ; 
3 pi es v el, , i> Poe Se = 
= ee AG = x’) dx, dx,’ °° A(x — x’) 


= 0, 


on account of the antisymmetry of A(x — x’) between x and x’. 

Now (0A4’/dx,) — (0A"/dx,) give the electromagnetic fields 6 and 
5¢, We see that the electric and magnetic field quantities are physical 
in the sense in which I have defined physical quantities. The potentials 
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are not physical. That corresponds to the physical quantities being 
gauge-invariant quantities. 

To express the constraints in terms of the Fourier coefficients, we 
have 


Ale) = i (Age? + Aine") a®k. 
We now form 


2 = —j / Ge Al em ee 

The constraints (15.1) give KXAS, ~ O and k*Ag, ~ 0. We can pass 
from the constant Fourier coefficients to the time-dependent ones 
because the same time factor is involved for all four of them. So we 
may write the constraint equations between the Fourier components 
in the form 


kK Auk = 0, ys ze (15.15) 


One can separate the waves A,, Ao, Ag into longitudinal and 
transverse parts. If we consider waves moving in direction 3, then 
Az, is a longitudinal Fourier component and A,;,, Ao, are the 
transverse Fourier components. It is convenient to count Ag also as 
a longitudinal wave. That separation is not relativistic. If we make a 
Lorentz transformation the separation would not be the same, but in 
spite of its not being Lorentz-invariant it is a useful separation for 
the discussion of various features of the theory. 

If we take waves moving in the direction of the axis 3, so that 
k, = ko = 0,k3 = ko, the equations (15.15) give 


Agr — Aar = 9; Aor — Age ~ 9. 


We shall make some deductions from these weak equations, remem- 
bering that we can multiply them only by factors on the left. We have 


(Ag + Az3)(Ao — Az) + (Ao + As)(Ao — Az) = 0. 


This gives us Apdo + Apdo — A3Ad3 — AzAd3 ~ 0. We have 
AgAo — Apdo equal to some c-number; this happens to be infinite 
if we take exactly the same k values in both factors, but that doesn’t 
bother us, because in any case it cancels with 4343; — A343. The 
result is that 274949 — A3A3) ~ 0. 
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We may now write (15.12) 


Hy = an? [UA + AgrAor + Asudan — AozAox) a’k, 
(15.16) 


ignoring an infinite constant. This Hy has no longitudinal energy 
when applied to a physical ket. I have proved it only for waves moving 
in the direction of axis 3. But the same thing is true for waves moving 
in any direction, because with this way of writing H, there’s nothing 
to distinguish axis 3 from the general axis in three-dimensional space. 
So we adopt this form for Hy, and then for any physical ket the 
longitudinal part of the energy vanishes. 

This mysterious negative energy now gets accounted for. The 
physical kets are the only ones we shall want to make use of in the 
interpretation of the theory; and for these kets, the negative energy 
associated with the Ap waves is just compensated for by the longitu- 
dinal part of the energy associated with the A,, Ag, A3 waves. 

We can express A,, = Au, + AX, the transverse and longitu- 
dinal parts, respectively. There are just two independent components 
among the transverse part, because k"A%, = 0. We have, with Hp 
given by (15.16), 


Hp = 4r’ i AL Aled 


This Hp applied to |0) gives 0, where |0) is defined to satisfy A7,|0) = 0 
and also to be a physical ket, which means 


k#A|0) = 0,  R*A,,|0) = 0. 


816. Quantum Electrodynamics 
Let us consider how the theory of the electromagnetic field gets 
altered when there are charges present. Beginning with the classical 


theory, we have to work with an action / equal to the action of the 
field plus the action of the matter, 


i= ic + £4) d*x, 


where £p is the action density for the field alone, given by (15.2) or 
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(15.3), and Ly is the action density for the matter, including the 
action density for the interaction between the matter and the electro- 
magnetic field. 

In all the theories of practical interest £,, involves the electromagnetic 
field variables only through involving the potentials and not through 
involving any derivatives of the potentials. I shall assume that now. 
There’s no reason why one should not consider a moze general case in 
which derivatives of the potentials appear, but this mcre general case 
has not yet turned out to be of any practical importance. 

Let us continue with our general methods and leave it undetermined 
how the matter variables occur in £,, or just what the matter variables 
are. We shall require, of course, only that £,, shall be invariant under 
Lorentz transformations, and then we have a relativistic theory. 

If we vary A” we get, with the form (15.3) for £p, 


[Canam + O8m 5 4) dt 
aA" 


= [ (Gora aE 22) OAC dee 


At 


él 


so 6f = 0 gives 
Axyv — 4rj, = 0, 


where j, = —0£,,/0A*. With the help of the Lorentz condition (15.1), 
it may be written 


Flag 4a je: (16.1) 


We get Maxwell’s equations if we interpret this j, as a current density, 
and we see how the current density is determined from the action 
principle: It is just the derivative of £,, with respect to A#, with a 
minus sign. 

We can proceed to pass to the Hamiltonian formulation by the usual 
method. With £y, not involving derivatives of the A, , the introduction 
of momenta conjugate to A, goes the same as in the case of no charges. 
We have again B, = A,°° satisfying the Poisson-bracket relations 
(15.4). Denoting by g, the dynamical coordinates needed to describe 
the matter, whatever they may be, and by p,, their conjugate momenta, 
we get for the Hamiltonian 


B 
eS I(- pee ac Pies i Lp oat i) doe 
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We may conveniently write it as a sum of two parts: H = Hp + Hy,. 
Hy comes from the first and third terms and is given by just the same 


equation as in the case of no charges, (15.5). Hj, consists of the 
remaining terms 


Bin = | conte — Lm) d>x. 


That's as far as one can go in completely general terms. One may 
specialize by assuming that £,, is linear in the potentials; that is, 
those terms that do contain potentials are linear in them. In the 
important examples it is linear in the potentials. Then j, is independent 
of the A", and we have the current density involving only the matter 
variables. Under these conditions, of course, 


Lm = —A,j* + terms independent of the A, . 


If j, does not contain velocities, Hy, consists of an interaction part, 


Hg = if Af dx (16.2) 


plus a part involving only the matter variables. 

In quantum electrodynamics we consider the electromagnetic field 
interaction with electrons and positions. We have already considered 
the electromagnetic field by itself in Section 15 and electrons and 
positrons by themselves in Section 11, and our problem now is to put 
them together and to have them interacting. For that purpose we 
need the total Hamiltonian 


ip Sse ites hy (16.3) 


where Hg is the Hamiltonian for the electrons and positrons alone, 
Hy the Hamiltonian for the electromagnetic field alone, and Hg the 
interaction energy. Hg is a function of y and y alone, and Hy is a 
function of a A,, B,, the Hamiltonian variables describing the electro- 
magnetic field. The interaction energy Hg is a function of all the 
variables y, ¥, Ay, Bu. 

Our problem now is to decide what Hg shall be. In the classical 
theory we considered the case when the interaction is described by a 
term in the action integral involving the potentials only and no deriva- 
tives of the potentials, and we found for Hg the expression (16.2) if 
ju is independent of the A’s and of the velocities Gn'°. These condi- 
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tions for j, are satisfied by (11.9), which gives 


Hg = —e i A“Da, d°x. (16.4) 


Let us take over this expression as it stands into the quantum theory. 

The classical theory is only to be used as a hint for getting a quantum 
theory, and we have to check and see that we really have a relativistic 
quantum theory when we work with the Hamiltonian (16.3) with 
Hg given by (16.4). We have to check whether the field equations to 
which it leads are relativistic, and we must check whether the com- 
mutation and anticommutation relations to which it leads are relativistic. 
Now the question of whether the commutation and anticommutation 
relations are relativistic is something we discussed in a general way in 
Section 13, and we found that when the interaction Hamiltonian does 
not involve derivatives of the dynamical variables, it does not disturb 
the commutation and anticommutation relations. Well, there are no 
derivatives occurring in expression (16.4) for Hg, so this is the simple 
kind of intersection which does not disturb commutation and anti- 
commutation relations and it is quite all right from that point of view. 

We still have to check whether our field equations are satisfactory. 
We work them out from Heisenberg’s equations of motion. Let’s 
take all the different variables we have in our problem, and see what 
equations of motion we get for them. Start with A,. A, commutes 
with Hz + Hg, so 


aA 
ae = ae Hr] = ie Hy] = Be 
which is the same formula that we have for the electromagnetic field 
alone. The B,’s retain the meaning of being the time derivatives of the 
A,,’s. The interaction does not disturb that relationship. 
Now let us form 


OB 

ae = (Bie Hy] a [Bi 5 Hy] a (ay HQ); 

the electron part of the energy does not contribute anything. The term 
[B,, Hr] is the same as with no interaction and is thus 4,’. To 
calculate [B,, Hg] we see that B, commutes with j, and that 


AR B,x'] aad =a ea, 53(x = 2% 
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sO ([B, Ho] = 4nj, and 0B,/dx9 = A," + 4nj,. Now when you 
combine this equation with our previous B, = 0A,/dxXo you get 

0°A, 
0X02 


eA eA (16.5) 


the Maxwell equation for the electromagnetic field produced by the 
current density j, . 

Now let’s take one of the components of y and write the Heisenberg 
equation of motion: 


oe ; 

it 5, Vez — inlWaz, Hy) = [Woes Ag + Hg)_. 
We have to compute the commutator now. To do that we multiply 
Hy + Hg by az first on the left, then on the right, and subtract. 
Now 


Hy + Hg = | Patast—ing” + eAW) + man — eAop)s ax, 


SO we get 


in a Yo = if (Ya, Vo lylar(— ity” + eA,'p’) 
+ mam’ — eAo’W’]y d?x’ 
= / Sav 63(x — x’o,(—ih'” + eA,'p’) 
+ mom’ — eAo'b’y dx’, 
and the result is 


ih ses = [o,(—iy” + eA) + mop — eAoble- 

This equation is precisely the ordinary one-electron equation involving 
the g-numbers A, for the field. We know from the one-electron theory 
that this equation remains invariant under Lorentz transformations. 
Therefore, the Heisenberg equation for az is all right. The equation 
for ¥ must be just the conjugate complex equation, and that must also 
be all right. 

So all the field equations that follow from the total Hamiltonian 
(16.3) are relativistic, as one would expect from classical analogy. 
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The interaction energy given by the classical theory is apparently the 
correct one to use. 

The Hamiltonian by itself does not give a complete theory; to get a 
complete theory we must bring in the constraints. For the case of no 
charges we had the constraints 0A,/dx, ~ 0. We must determine if 
they are still all right or if they get disturbed by the interaction. The 
question of whether they are all right or not is the question of whether 
they are consistent or not. . 

In the case of no interaction we checked the consistency in two 
ways, one of them consisting of taking them at two general points in 
space-time and working out the Poisson bracket, from (15.14); we found 
that it vanishes. Well, we can’t do that now; we can’t take the 
dA,,/dx, at two different times and compare them, because that would 
involve integrating the Heisenberg equations of motion, which is a 
very complicated problem. We have to check the consistency con- 
ditions without integrating the Heisenberg equations of motion and 
must use the other method, which consists of taking all the constraints 
that hold at one time and seeing whether they are consistent with the 
equations of motion and with one another. 

At a particular time we have 


dA, ddA. a AA, 


ma 


Cra or Oxgnon gst Oxo" OX, 
and so on. Now if you look at these equations you see that, as in the 
case of no charges, only the first two are independent when you take 
into account the field equations (16.5). The third one is a consequence 
of the first two, because it gives 


Aytt” + Anjt = 0. 


Now j,," vanishes because of conservation of charge. To prove it in 
the present theory we refer to the one-electron theory. We know that 
the charge is conserved in the one-electron theory. From that we can 
deduce the conservation of charge in the present theory, because the 
field equations for y in the present theory are of exactly the same form 
as in the one-electron theory. Just taking over the one-electron theory 
and making the potentials into g-numbers instead of c-numbers does 
not in any way affect the proof of the conservation of charge. 

Well we are left with 4,47" ~ 0, which is valid as a consequence of 
A," = 0 holding at all points in three-dimensional space. The result 
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is that when the field equations are used, the constraint A,"°° ~ 0 
doesn’t give us anything new. Similarly, the higher derivatives with 
respect to xo don’t give us anything new. The fact that these higher 
derivatives don’t give us anything new is the check that the constraints 
are consistent with the field equations, and thus that they commute 
with the Hamiltonian. 

We are left with two constraint equations: each of them is a three- 
fold infinity of conditions, because each of them holds for all points 
X1,%X2,%X3. We let us write them out in terms of our Hamiltonian 
variables and check the Poisson-bracket relationships. The first gives 
us 


Bo —+- A," = oi = 0, (16.6) 


as in (15.10). The second gives 07A9/dx0” + B,” =~ 0. Now, using 
the field equations, we replace the second time derivative by Ao’ + 
4mjo and write the constraint 


(Ao” + B,)” + Anjo = do = 0, G7) 


like (15.10) with the extra term 4zj9. You will recognize that equation 
as —div & + 4njo ~ 0, one of the Maxwell’s equation for &, the 
electric field. 

To check whether these constraints are consistent, we have to check 
whether the equations (15.11) still hold with the extra term 42/9 in 
goo. The first two are obvious. For the third we have 


[25 62] = [(Ao” + By)” + 4ajo, (Ao + Be’) + Anjo’). 


We get one new term 16r7[j9, jo’]. So our work for the case of no 
charges tells us that everything is all right, except for this new term 
involving the charge density at two different points. Now what is this 
Poisson bracket? There is no need to bother to work it out in detail, 
because we can see on general grounds that it has to vanish. The anti- 
commutator of y and y’, in which they are not differentiated, involves 
only the delta function: so [jo, /o’] = some coefficient times 63(x — x’). 
Now [jo,/o’] has to be antisymmetrical between x and x’. Any 
multiple of the 6-function is symmetrical unless it 1s zero; hence 


(Vosvo l= 0: 


So we have checked that the constraints commute weakly with the 
Hamiltonian and with each other; actually they commute strongly 
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with each other, as in the case of no charges. So we can still give a 
meaning to the reality of all g-numbers, not merely the physical ones. 

There are two important applications of this theory of quantum 
electrodynamics: the anomalous magnetic moment of the electron and 
the Lamb shift. These are applications that one cannot make with a 
more primitive theory. One needs the full quantum electrodynamics to 
get the results. I shall now show how the calculations of the anomalous 
magnetic moment and the Lamb shift are to be made by working from 
the Hamiltonian (16.3) and (16.4). 

We have to solve the equations of motion, and as I have discussed 
before, it is no use trying to solve the Schrédinger equations of motion. 
We might set up formally a ket vector satisfying the Schrédinger 
equation, but that ket vector would not stay in Hilbert space. The 
only way to carry on with this work is to confine our attention to the 
Heisenberg picture and obtain solutions of the Heisenberg equations 
of motion with the Hamiltonian (16.3) and (16.4). 


§17. Solution of the Heisenberg Equations of Motion: General 


Now let us first discuss in general terms what solutions of the Heisen- 
berg equations will look like. If we have any Hamiltonian H, we can 
take the Heisenberg equation for a dynamical variable & i#(dt/dt) = 
fH — H, and try to get a solution of it. That solution would be in 
the form & = f(t, constants of the motion). These constants of the 
motion would be q-numbers, noncommuting, of course. What should 
we take them to be? 

Well, we may take them to be the dynamical variables at one par- 
ticular time, so as to get 


E = f(t, q(to), P(to)). 


We might be given all the field quantities at time fg and calculate the 
field quantities at a general time in terms of them. | call this method 1A. 

Of course, you might say that from the relativistic point of view this 
is not very good. It means taking £, a field quantity at a general point 
in space-time, and expressing it in terms of field quantities over all 
positions of space at one particular time. That is a very unrelativistic 
way of proceeding. 

Alternatively, we might express £ as a function of the four x’s and 
of the ingoing fields, to get Exorjzezg = J (Xz, ingoing fields). This I call 
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method 1B. How are ingoing fields defined? We define them as in the 
classical theory. Consider the electromagnetic field, for instance. We 
have the equation : 

LJAu = Agi 


holding in the quantum theory as in the classical theory. We can solve 
this equation in the quantum theory as in the classical theory and 
express Ay = Ayin + Apret, Where A, ret is defined in terms of j, by 
the classical formula. The noncommutation of the j,’s at different 
places doesn’t matter here because everything is linear in the j,’s, and 
A, ret Appears as the usual linear function of the j,’s. We then define 
ae a ee We then have |A,3, = 0. 

That is how one defines the ingoing field in the electromagnetic case 
and one can do a corresponding thing for the y fields. These definitions 
go just the same as in the classical theory, because they involve only 
linear relationships. 

You might try to solve the Heisenberg equations of motion by 
getting each field quantity at a particular point in space-time expressed 
in terms of the coordinates at that point in space-time and the ingoing 
fields or the Fourier coefficients of the ingoing fields. You might say 
that this is a nice method to work with because it is relativistic, and it 
would enable one to set up the whole theory in a manifestly covariant 
form. However, the ingoing field is a rather complicated concept when 
one starts from the Heisenberg equations of motion. Its existence 
depends on the convergence of the integral for the retarded field—an 
integral that goes back to the infinite past. There are examples for 
which the concept fails. A further discussion of this question will be 
given in Section 32. 

It is not desirable to have the theory depending at this early stage in 
its development on such a shaky concept as the ingoing field. For this 
reason one should reject method IB and fall back on IA. We can im- 
prove on method IA by adopting a different method, which I call 
method 2. 

Method 2: Find a quantity K involving ¢ explicitly such that it is a 
constant of the motion. Express K as a function of ¢ explicitly and the 
q’s and the p’s, the dynamical variables, at time t. The advantage of 
method 2 compared with IA is that there is only one time ¢ involved. 

The quantity K involves the time explicitly and so it doesn’t satisfy 
the same equation of motion that ¢ satisfies. The equation for § is a 
valid equation only for dynamical variables and functions of dynamical 
variables that do not involve the time explicitly. K does involve the 
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time explicitly, so it must satisfy the more general Heisenberg equation 


dK Oe 
Feo KH — HK + in 


The extra term comes from the explicit dependence of K on t. 
The requirement that K is constant means 0K/dt = 0 or 


ih 


in ox = A= KH. (aa) 


We get this new equation of motion to work with. You see that there 
is the opposite sign coming in, and also the total differentiation with 
respect to time is replaced by partial differentiation. You must 
remember that in this partial differentiation gq, and p, count as con- 
stants; i.e., 0g,/dt = 0, Op,/dt = 0. gq, and p, have total derivatives 
with respect to time but no partial derivatives. 

Method 2 is the method with which we are going to work. 

We now have to consider the problem of getting solutions of (17.1) 
with H equal to the Hamiltonian for quantum electrodynamics. We 
shall use a perturbation method which assumes that Ag is small. 
Hg involves the parameter e (the charge of the electron), and that is 
small compared with the quantity of the same dimensions constructed 
with the Planck constant. On the basis of e being small, we can use 
this perturbation method and express everything as a power series 
in e, and we just work with the early terms in that power series. 

The first thing we do is pass over to the interaction representation. 
I suppose you are quite familiar with the interaction representation 
as used in connection with the Schrodinger picture. One can use a 
corresponding interaction representation in connection with the 
Heisenberg picture, and it goes like this. 

Let us write H = Hy + Hg, where Ho is the proper energy. For 
any q-number a, define a* = e'#0t/%ge—*#ot/*® This is just an ordinary 
unitary transformation, because Hy can be considered as an oper- 
ator in Hilbert space. 

Let us apply that transformation to the quantity K and make a 
partial differentiation with respect to f: 

2K 
in on sg Oe (—Hok + ine + KH) me 
Se Cele Wile = Ge tye 
= Oe Hinks _ KHjes = Hoke Kip 


(17.2) 
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This is the equation of motion for K in the interaction representation. 
It achieves the usual result which one achieves by going over to the 
interaction representation: We get equations of motion in which 
things vary only on account of the interaction; the proper part of the 
energy doesn’t contribute anything to the variation. 

This passage to the interaction representation is exact. We are now 
going to make use of the assumption that the interaction energy is 
small and express K* = Ko* + K,* + Ko* +.---, the various 
terms being of different orders of smallness in e. Now we substitute 
this expression for K* in the equation of motion (17.2) and equate 
terms of corresponding orders of smallness. We get 


* 

e oto ay (17.3) 
* 

ih ont SS aH oe (17.4) 
* 

if os Se Ke ek Ge. (17.5) 


and so on. 

Now what do these equations mean? First of all, (17.3) shows that 
Ko* must not involve the time explicitly. So Ko* equals a function of 
the q;’s and p;’s, the dynamical variables at time f. We choose a suitable 
Ko* to start off the calculation. Then we proceed to get K,* by 
evaluating the [Ko*, 4 *]_ commutator and making a time integration. 
We proceed to get Ky* by evaluating the [K,*, Hg*]_. commutator 
and making another time integration, and so on. In that way we get 
a solution of the Heisenberg equations of motion by this step-by-step 
method. 

When we have got these K* quantities we can go back to the K 
quantities by just inverting the unitary transformation, and then we 
have our solution. One point I want to emphasize is that Ko* involves 
only dynamical variables and does not involve the time explicitly. 
When you go back to Ko you will get something that does involve the 
time explicitly. 


§18. Solution of the Heisenberg Equations of Motion: Quantum Electro- 
dynamics 


We shall be concerned with the one-electron problem in the calcu- 
lation of the anomalous magnetic moment and the Lamb shift, and 
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for that purpose we shall take Ko* to be the operator of the emission 
of an electron into a certain state. 

It will be necessary in these calculations to suppose the existence 
of a static magnetic or electric field. To calculate the anomalous 
magnetic moment, we suppose that there is a static magnetic field, 
and then we examine how this static magnetic field influences the 
energy of the electron that is being emitted. Similarly, for the Lamb 
shift we suppose that there is a static electric field and we determine 
how that static electric field influences the energy of the electron that 
is being emitted. 

In our derivation of the Hamiltonian we didn’t consider static 
fields, but we can bring them in in a trivial way. A static field is 
described by potentials @, which don’t vary with the time; they are not 
dynamical variables; they are just c-numbers, functions of x1, x2, x3. 
So to bring in the static fields we add on these @’s to the other 4’s. 
Thereby we bring in some extra terms to the Hamiltonian. 

Let us simplify future work by taking # = 1 as well as c = 1. 
When one gets final results, one can fill them in to make the dimensions 
right. 

We shall consider how to introduce the absorption and emission 
operators which are going to be our working variables. They must be 
defined as dynamical variables that do not involve the time explicitly. 

Let us take first the electromagnetic field variables. For them the 
work is the same as in Section 15 for the case of no charges. We start 
with A,, and B,, = 0A,2/dX9 and make a Fourier resolution like 
(9.4) and (9.5). It reads 


Ae i [Apney 9 eee slid ie 


(18.1) 
Bas = if ile = Ae | d°k. 


These equations define A,;, and A,,. We just have to invert the 
Fourier resolutions and we get two equations which determine A,;, and 
Ane 

What I'd like to emphasize about these equations is that there’s no 
time t occurring in them anywhere. The equations that define A, 
and A, don’t involve the time and so we have (0/df)A,, = 0, 
dA,1/0t = 0. This way of introducing Fourier coefficients gives 
Fourier coefficients that are not explicitly dependent on the time. Of 
course, the total (d/dt)A;, is something which is not zero. 
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We know what that total dependence on the time is in the case of no 
interaction. For no interaction A,, is proportional to etlklt, so we 
then have (d/d1)Aux = i|k|A,.. This means 


[Aus, Ho) = i|klA,,. (18.2) 


This equation must also hold when there is interaction. It just tells 
us something about the commutation relationship of A, with Ho, 
which must hold independent of whether there is interaction or not. 

Alternatively, one could obtain (18.2) by working out the Poisson 
bracket directly. We have 


[Ane 0 
and 


(Aye, Axel = — gS aa(k — k’) (18.3) 
from our work on the electromagnetic field by itself, and the com- 
mutation relations are not disturbed by the interaction. Equation (18.2) 
is a consequence of the commutation relations (18.3) and therefore 
holds also when there is interaction. 

Equation (18.2) can be written 


ArHo => (Ho ee |k|).A;. 
We can generalize that to 


Axf (Ho) am S(Ao = |k|)A;, 


for any function f; it holds certainly for a power series, and we don’t 
need to apply it to other functions. We can use this now to calculate 


* LH yt —itll gt Tot —i(fy—ikl)t i\k| ¢ 
A; ae Ape os = 2 WE HED) Ae ae 


There we have calculated the starred field variables that occur in the 
interaction representation. 

To introduce the absorption and emission operators for the electrons 
and positrons, we must take into account the static electric and magnetic 
fields. They give a further term in the Hamiltonian: fax, d°x. This 
further term is going to be put in Hy and not Hg. The perturbation is 
going to remain the same as before when we had no @, field. 

With the term {@*j, d?x included in Ho, we shall be wanting to use 
the wave functions for the one-electron problem in the static field. 
Let’s suppose that these functions have been worked out. Finding 
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the different eigenfunctions of the energy for one electron in the static 
field is just a trivial problem, without any field quantization involved. 

In this one-electron problem let £,, be an eigenvalue of the energy 
and (xa|n) the corresponding eigenfunction. The eigenfunction is a 
function of the three x’s and a function of a; a takes on the four values, 
corresponding to the four components of the wave function. This 
n Is just a symbol to remind us of which eigenfunction we are talking 
about. I’m going to write down the equations as though the energy 
spectrum is discrete; if it is not discrete, it just means the usual change 
in the way of writing, with the sums over the discrete n values replaced 
by integrals. 

I would like to write down now the wave equation which the eigen- 
function satisfies. We have as the energy operator for one electron 
in the static field, 


x (= 


We take the ab matrix element of (18.4) and multiply it into (xb|7). 
We get 


) + amin — eQo. (18.4) 


E (-i 2. + eG) + a,m — e00| le) = 8 Cale. 
Xy jab 


This equation is the usual equation for determining eigenfunctions of 
the operator (18.4). I am going to suppose that we have worked out 
all these eigenfunctions and we’ve got them all normalized, so that 
they satisfy the usual orthonormal relations: 


J caxay¢xaln' os San ee (18.5) 
a summation being implied over the four values of a; and 
D> (xalny(n|x’b) = Bap 63(x — x") 


with (n|xa) the conjugate complex of oa 

Now let us define y,, = f(alxa)var d°x, a summation being implied 
Over a; Or we can invert that equation and write Von = cal 
The work with eigenfunctions previously has just involved c-numbers, 
and now I’m introducing g-numbers. These Wax S are g-numbers, and 
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I'm defining these new q-numbers y,,. We shall be using also the 
conjugate complexes 


Vn = | Faxtxcley d°x, Wan = S Yn {n|xa). 


Here we introduce the g-numbers y,,, Wn; what are their anticom- 
mutation relationships? One very easily finds that 


[vn ; ve = CL vn ; Vridy ona 


Now we started off with ¥,,, and that is something which does not 
involve the time explicitly: (0/f)Y.z = 0, (0/8t)Par = 0. The equa- 
tions that introduce y,,, ¥,, do not involve the time explicitly; so we 
can infer that dy,,/dt = 0, dy,,/dr = 0. 

Now Ho, the proper energy, is equal to Hy + Ep, the electron part 
and the electromagnetic field part. Hy, the electron part, is 


Ay = | doe (-: 


When we go over to these new y,,, J, variables, we get 


ts =) [tatetsay [ar (rh + 0) + nm — 0] 
X (xb'|n' nr Bx. 


re) 
5 se eG) + a,m — ec| Wor Ox. 
X> jab 


Now we make use of the equation the eigenfunction satisfies, so that 
we get 


iy = YS | Batnbxa) Eyal Wb de 


an! 


Now the only things involving x here are the eigenfunctions, and we 
can make use of the relation (18.5) and get 


Ay = > Vee Onn nr! = SS Erbin : (18.6) 


nn’ n 


This is what you would expect the energy of the electrons to be, referred 
to the y,,, YW, variables. 

The variables y,, and y,, are the operators of emission and of absorp- 
tion into the different states n. For £,, greater than 0, y, is the operator 
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of absorption of an electron and y, is the emission of an electron; 
and for E,, less than 0, ¥» is emission of a positron and Y, is the absorp- 
tion of a positron. That is the physical meaning of these yn, Vn 
quantities. 

Now we must do a piece of work corresponding to what we did in 
the case of the electromagnetic field variables, to find the starred 
quantities in the interaction representation. Again we can get the 
results simply by considering the case of no interaction, when 


Wy 
I dt oe 
So [¥n, Hz]. = Enn- This equation holds also when there is inter- 


action. We can deduce from it ¥,H9 = (Ho + £En)¥n; and more 
generally ¥nf(Ho) = f(Ho + En)\¥n, 80 


Wnt = ef nem= 0! = e tEnty and Vr* = goo 


I’m going to use an abbreviated notation for the Fourier resolution 
of the electromagnetic field, 


= | A emg dale (18.7) 


I’ve introduced this new index v, which takes on the values +1. For 
y = 1, we define A’) = A,,3 for vy = —1 we define 44,” = Ay_y- 
Furthermore, an average over both v-values is implied, going with the 
integration over k. It is more convenient to work with an average 
over both values of y and not a sum, because if vy drops out of a term, 
you can forget all about »v. The point of this notation is that we can 
write the Fourier expansion of A,, as a single term instead of two 
terms. You might say that it is not a very great saving in writing 1f we 
put down just one term instead of two—but we shall have to deal with 
quadratic functions of the potentials and then this notation will enable 
us to write down one term instead of four, and that’s a considerable 
saving. 

In terms of the absorption and emission operators, Hg appears as 


Ho = —2e ii Alc 'em Vd dike ete Sl i AQ Bkdn 


Mn 


X (alate In a. 


Here I’ve introduced a new kind of quantity, (njate~*™” |n’), which is 
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just a matrix element for the one-electron problem, defined by 


| clax)ecne"(oxln' are 


We now need Hg*, to appear in our perturbation theory: 


Ho" = —2e >) i (alate In!) BRAG Dn Wn 


nn’ 


—2e Xf (alate ">In" PAY Dae CT Bn Bt 


nn’ 


I have used there the formula 4%* = A®e*!*l’ You see, when v 
is equal to +1, we need the factor e''*!*; when v is equal to —1, we 
have the conjugate complex thing and we need the factor e~"'*!": so if 
we put eK! it’s correct for both values of v. 

Now we are interested in a one-electron problem. We therefore take 
our Ky* to be the operator of emission of one electron into a positive 
energy state. Let us take Ko* = y,, the operator of emission of the 
electron into the initial state 7. We assume the energy of state 7 is 
positive; otherwise Y; would be the absorption operator of a positron 
instead of the emission operator of an electron. This choice of Ko* 
makes 0K,*/dt = 0, so it is suitable for starting off the perturbation 
calculation. 

We now proceed to find Ky*: 


.OKy" E 
Came ee he I — [Ho ye 
= —2e Xf tolate 0" BPKAY 
x [v. y : HE) gO 
Now 


nw’, Vile = Vrbn Wi os Winn’ = YnWnWi at Witn') 
= Wn On'i- 


We have 6,; coming in, so we substitute 7 for n’ and drop out the 
summation with respect to n’. We get 


dKy. i poi, x)) \ 73, 4%) 7 ptVlkKI+E,—E)t 
i = —2e nla’ e 1) AKA LE ne aa 
_ Df (later Phi) ALT 
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Now we can carry out integration with respect to 7; everything is con- 
stant except the exponential factor, so the integration is trivial: 


i 7 ge Wega 
Ke == Wp | (late ale PKA Vn v|k| Seep ee— Eo 


Now we must go on to Ky*. We have 


I 


[Ho", Ke te 
he > flare ™ ty ahote ak dk! 


nnin'’? 


i(0K»* /dt) 


oiolkl4+B,—E pt 
vik| + E, — E; 

X [AS burda, Ave nd. (18.8) 
All we have to do now is a time integration to get K2*. That would 


give us the solution of the Heisenberg equations of motion to the 
second order. 


« ST | Ey Egret 


§19. Normal Ordering 


Now this solution involves, of course, g-numbers, and would take 
different forms according to the order in which we arrange those q- 
numbers. We need to write this solution in a certain standard way 
and for that purpose we introduce a normal ordering of the g-numbers. 
We take it to be the ordering in which all the emission operators are 
to the left, all the absorption operators to the right. I take over this 
normal ordering from the usual work in the Schrédinger picture and 
assume it is physically significant. Some justification for it will be given 
in Section 31. 

We must now, therefore, arrange the g-numbers in (18.8) so as to 
put all the emission operators to the left, all the absorption operators 
to the right. This rearranging of the factors, to get them into the 
normal order, will bring in some extra terms, and we must see just 
what they are. It does not matter about the order in which different 
absorption operators appear, or the order in which different emission 
operators appear, changing that order would at most involve a change 
of sign when we are dealing with fermion operators, and such a change 
does not bring in any extra terms; that change is to be looked upon 
as trivial. 
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The absorption operators are y, for E, greater than 0, y, for E, 
less than 0. Now what about the electromagnetic field variables? In 
the case of the transverse variables, we have operators of emission and 
absorption of photons, and there it is quite clear what is to be an 
emission operator and what an absorption operator. We get the A; as 
absorption operators. How about the longitudinal variables? Well, 
there it is not clear what are emission operators and what are absorption 
operators. We shall again take A, as an absorption operator. This 
should be looked upon as a definition of absorption operators in the 
case of the longitudinal variables. It is a convenient definition to work 
with, and makes it clear what the normal ordering is in all cases. 

Question: Could you expand your explanation for the need of the 
normal ordering? 

Answer: We want to get some numerical coefficients to which we can 
give a physical meaning. Now, you do not get any definite numerical 
coefficients at all if you leave the ordering unspecified, because changing 
the order changes the numerical coefficients. A normal ordering is for 
the purpose of getting well-defined numerical coefficients, which you 
have to do when you want to get something that is to be compared 
with experiment. g-numbers are not something you can compare 
directly with the experiment; you must take out numerical coefficients, 
and the only way to get numerical coefficients is to have some definite 
system of normal ordering. It is, of course, an assumption that this 
normal ordering is physically significant. I shall discuss this physical 
assumption later, in connection with the interpretation of the theory. 

We adopt this definition of normal ordering and we assume that it is 
physically significant. We must now see what extra terms arise when 
we do the normal ordering. 

We have the commutation relations (18.3) for the electromagnetic 
field quantities. Consider the product Aiea We see that if vy and 
v’ are both equal to | or —1, this is already in the normal order. Also 
A DAR? is in normal order. The only case in which we do not have 
normal ordering is where we have y = —]1 on the left and »’ = | on 
the right. And 4%, AX}. differs from the normal ordering of the factors 
by the commutator in the second of equations (18.3), so we get 
AG Axe, = normal order — g,,/4m7|k| 53(K + k’). We have to 
change k — k’ tok + k’ in the 6-function because of our definitions 
of Aj? with various v-values. We can write one formula that includes 
all these results, 


@) gy — Bur ai le ey laa 
A,tAxx? = normal order — 4n2|k] 63(k + K( D )( 5} 
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There are different ways of writing the v-factors, because 
d= yd $2) =] (=e) = a ae 


These are just algebraic formulas that hold because v and »’ are able 
to take on only the values | and —1. 

Now consider the product ¥,W,-. This is in the normal order for all 
nif E, is negative. In general we have 


Wn¥n1 = normal order + dyn OCEn’). 


The @-function is the step function, which is used quite a bit in field 
theory. It is defined by 


OC) — a ele 8 NG), Ox) = 0 fon 620; 
Alternatively, @(x) = 4(1 + x/|x|). Note that 
OQ + O(-—x»=1, O() — O(-x = x/lx1. (19.1) 


We can also write ¥,W,. = normal order + 6n,/O(E,,), because nn’ 
have to be the same for 6,,- to contribute anything. Then we have the 
corresponding formula 


YrWn = normal order + 6n,/O(—E£,.). 


We must now take the commutator in (18.8) and put it in normal 
order. If we write it out in full it is like this: 


(¥!) giv) 7 7, DMO TT 
AX Anita Vn Wn —s AQRA Wnt Writ 


Putting it in normal order using the formulas for normal order given 
above, we get 


_ _&r 
472\k| 
x x@ ate Vn’ bnnv O(E,) =e Vn bain O(— E,)] 
=F @ — Vn Buin O(— En’) a Vn’ bnn O(— E, )]} 
-+ extra terms, 


ORCS rd) = 21) 


the extra terms involving niore than one emission or absorption operator. 
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We must substitute this commutator in (18.8), and that gives us 


0K" e° 
g a Ye V(E,—E;)t $ (UE, — las 
1 yi a Wn’ ae n 1) ale > Vale E;)t 
T ; v 
nm n 


+ extra terms 


2 
@ = : - 
ee Me Ynt Me ats Zale om! ap Sur) terms, (19.2) 


where the coefficients Y,; and Z,,; have the following values: 


2 : ey: 
oe = (ay Ff erage" (alate aed 
n 


v|\k| + E, — E; 
Pak 
ala OC ae = OF, TA (19.3) 
and 
Line = (ee @(—Ey) | Glage"™™In'yinlote "Up 


yp — yp! d°k 


x Ik E, — Ei Tk 


(19.4) 


There is one immediate simplification we can make. We have to 
average over all values of v’; that means that the terms containing v’ 
linearly drop out. Thus the vy — vy’ factor in the numerator, for both 
Y and Z, can be replaced by v. 


§20. The Energy Change 


The extra terms in (19.2) all refer to more complicated physical 
processes with emission and absorption of more than one particle, 
whereas the Y and Z terms just mean a modification in our original 
emission operator produced by the perturbation. Suppose we pick out 
the one term in (19.2) for which is equal toi. Itis (e?/m)(Yu + Zi yj. 
This term represents a change in the energy of the emitted particle. 

Let us see just what energy change is represented by the coefficient 
of y; here. The most definite way of getting the desired result is to 
consider the effect of another perturbation in the Hamiltonian. Con- 
sider a perturbation V of the form V = >onWnbr, where b, is some 
small c-number for each value of n. Let us suppose a perturbation like 
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this to be added to our Hamiltonian. Such a perturbation will obviously 
cause a change in the emission energy of amount 5b,, because, 
from (18.6), the part of our Hamiltonian that refers to the particles 
only is just > Wv,nEaw,. Since &, gets changed 104 sa 0,mel 
means that 5, is the small change in the energy of the emitted particle. 

If the problem is handled by a perturbation method, this V gives 
i(0K ,*/ot) = V*¥Ko* — Ko*V*. Talang ko, — ys bel one, 
i(0K ,*/dt) works out to be (0K ,*/dt) = yb; to first order in the b’s. 
This 0K ,*/dt can be compared with (19.2), and it enables us to infer 
that the change in the energy of the emitted particle is precisely the 
term forn = 7: (e?/r)(Yix;, + Zi). 

This is the expression that we need to evaluate to get the anomalous 
magnetic moment and the Lamb shift. We have the general formulas 
(19.3) and (19.4) for the matrix elements of Y and Z, and it is the 
diagonal elements that will be the important ones. 

We can write Z in the form 


—i(k,x)| - 3 
ta all Bb (nloye— |i) dk f 
Znai= 7 ie KES gene on (2019) 


where Jy? = $>0,/O(— En )(n’ |ate'™” |n’). Now J, can be written 
as a diagonal sum, 


Je! = ${(ate™ @(—E))), (20.2) 


where ((_ )) means diagonal sum and E is the energy operator for one 
particle in the external field, 


E = a,(p; + e@,) + amm — eQo. 


E is not a q-number here, but it is a matrix for the one-electron 
problem. It is a matrix whose elements are all c-numbers. O(—£), of 
course, is (1 — E/|E|). How is |E| defined? Well, this is the modulus 
of a one-electron operator, a fairly simple operator. We have 
|E| = +(E)"’*. So |E| is quite a well-defined operator, perhaps a 
bit awkward to work out; in fact, the only practical way to work it 
out in general is to use a representation in which E is diagonal. 

The diagonal sum of a matrix is the same in all representations. 
The diagonal sum (20.2) which we find to give the correct expression 
for J," when referred to the representation with E diagonal, can be 
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evaluated more conveniently in that representation in which the 
momentum variables are diagonal. Let us proceed to work in the 
representation in which the p’s are diagonal. We have 


(pile |p") = 3(p! — p” — bk). (20.3) 


This follows from p,e°” = e®(p, + k,). So we ell 


1 a(K,x)) 7 , " , 
Jt = Lif ort *|p')(p'etO(—B) |p") dp! d°p" 
1 
= LE Gp at@(—E)Ip' + k)d'p’. (20.4) 


Now let us take first of all the case of no external field. For @, = 0, 
E and ©(— E) commute with the three p’s, so that J,” has only diagonal 
elements when referred to the momentum representation: thus i“ = 0 
fork ~ 0. Fork = 0 weare not interested in the value of J;“, because 
our whole perturbation theory is concerned only with the values of k 
that are not zero. The zero value of k would correspond to static 
fields, which are not included in the perturbation. 

Let us now consider the case of a uniform magnetic field. For a 
uniform magnetic field, the @,’s are linear in the x’s. Then (p’|E|p’’) = 0 
for p’’ differing from p’. Under these conditions |E| also involves no 
matrix elements differing appreciably from the diagonal, but only ele- 
ments infinitely close to the diagonal, and we again have J," = Ofor k 
differing from 0. Thus the Z term contributes nothing to the anomalous 
magnetic moment. 

For the Lamb-shift calculation we have an electrostatic field that is 
not uniform, and in that case we do have off-diagonal elements con- 
tributing to J;”. 

The Z terms are called the vacuum-polarization terms by the people 
who work in terms of the Schrédinger picture and Feynman diagrams. 
I do not know whether that name has much significance when one works 
in the Heisenberg picture, but still it is good notation. The vacuum- 
polarization term contributes nothing to the energy of the electron in 
the absence of a static field and nothing in a uniform magnetic field, but 
it does contribute something in a nonuniform electric field and thus 
gives a contribution to the Lamb shift. We shall have to work it out 
later when we come to work out the Lamb shift in detail. 
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Well, let us now study Y. We can simplify Y a little using the proper- 
ties (19.1) of the © function, and then it becomes 


Hey a = 6a D | rte ln) (n|o" mo py v|k| + TA =a E; 


E,, \ a®k 
; 20.5 
x (> TE, +) Tk ed) 


We can write this in a more concise form, 


_ a(k,x ate ENE —i(k,x ak 
5 oe | oe sed iki + oe Wes alt) TE 


We can make a further simplification if we replace vE by |E| in the 
denominator, which is permissible because the numerator vanishes 


when vE = —|E|. Y,,;can be looked upon as the n, i element of some 

operator 
= =! ide) YE (E/|El) —i(k,x) ee 
Y = (47) foe 1 == (l= vp” a ik] 
(20.6) 
For any given electron operator », we define fi = e7'” ye’, 
Then 
vt (E/E) ak 
Y = (40 tf . va 20.7 
a ae” 


There is a change in sign in k in (20.7) compared with (20.6), but that 
does not matter because k is just a dummy variable of integration. 
The sign in (20.7) is the more convenient one to take. 


§21. Cutoffs 


The expression for Y involves an integration over k, and when you 
proceed to carry out that integration you get an infinity coming in. 
This working with the Heisenberg picture has eliminated many of the 
infinities but not all of them. There are still some which survive, and, 
of course, if we want to have a logical theory, we cannot tolerate any 
infinities at all. What we have to do now is to make a cutoff in the 
interaction energy so as to get rid of this infinity. 
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We are changing our Hamiltonian, of course, when we put this 
cutoff into the interaction energy. Now you may object that changing 
the Hamiltonian in this way spoils the relativistic invariance of the 
theory. That is quite true; however, if we do the cutoff properly the 
relativistic invariance will be spoiled only for high-energy phenomena. 

Quantum electrodynamics in any case cannot be a completely accu- 
rate physical theory because it neglects all sorts of particles that turn 
up in high-energy processes, mesons, neutrinos, and so on. Quantum 
electrodynamics cannot be a complete theory in itself. We find that 
we cannot make quantum electrodynamics consistently relativistically 
invariant for high-energy processes. Well, that does not really matter. 
We shall find that we can make the cutoff at a fairly high energy value, 
of the order of several hundred million volts. It is reasonable to make 
the cutoff at an energy value somewhere there and we shall then have 
a theory that fails for processes involving energies of more than some 
hundreds of millions of volts, but which we can hope will be accurate 
for lower energies. We may thus hope that it will be accurate for 
calculating things like the anomalous magnetic moment and the Lamb 
shift, where there are no high-energy processes involved. 

I look upon it in this way: We have a quantum electrodynamics 
that cannot be accurate for high-energy processes, because we do not 
know how to take into account high-energy interactions, and our 
cutting off the interaction for high energies is just an expression of our 
ignorance of how high-energy processes should be taken into account. 

We have to be careful in how we cut off the interaction energy. The 
usual treatment of quantum field theory is quite sloppy over this point. 
I want to do it in an accurate logical way, and for that purpose we 
have to specify accurately just how the cutoff is to be made. When one 
looks into it, one sees that there is more than one cutoff process 
worthy of consideration. 

Our interaction energy is Hg = —efPa,wA"d*x. We express it in 
terms of the Fourier coefficients and we see that each term involves 
three particles. It involves Y,, Wn, AY , which are emission or absorp- 
tion operators for three particles. 

Now one way to cut off is just to put an upper limit to the magnitude 
of |k| in the interaction energy. | call it cutoff (a). 

Cutoff (a): Limit |k|, say |k| < g, where g is some large number. 
That is the most natural and obvious way to apply a cutoff, and of 
course if you do that for an integral like Y you immediately get a 
finite quantity, and are no longer bothered by infinities. However, it 
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turns out that this one is not so satisfactory, and one should consider 
other cutoff processes. 

Let us take cutoff (b). Here we limit the sum of the moduli of the 
energies of the three particles involved in each term of Hg; that means 
|E,| + |E:| + |k| < 2g. We shall be interested in the final state 
being equal to or rather close to 7; and the three energies involved are 
thus E,, the energy of the intermediate state; E;; and |k|. £; is a 
constant as far as the integration process is concerned, and so let us 
take it over to the other side and get |E,| + |k| < 2g — E;. That 
is a reasonable method of applying the cutoff. 

One should be careful, of course, when one thinks of various possible 
ways of applying the cutoff, that one preserves the reality of Hg. You 
must not cut off Hg in some way which would turn it into something 
that is not physically real. 

There is a third cutoff process to be considered. This is where we 
limit the sum of the kinetic energies of the three particles. Previously 
we limited the total energies. I am now going to limit only the kinetic 
energies to get cutoff (c). E contained a part —e@ which was potential 
energy and which we are now going to exclude in the cutoff process. 
That has the effect |E + e@o| + |k| + E; + e@o < 2g. 

There are thus various cutoff processes we might consider; one cannot 
say offhand which is the right one to use. We have to investigate each 
of them. We would like to have our cutoff process lead to physical 
results that are independent of just where we chose g. So we must 
carry on with the calculations and see what luck we are going to have 
in that respect, and we also have to decide later which of those three 
cutoff processes we are going to adopt. But already at this stage we 
can work out the differences these cutoff processes will lead to. I 
shall do that now and work out the differences in the expression Y 
according to these three processes. 

Let us go back to the integral (20.5). Cutoff (b) involves the restricting 
of variables by the inequality |E,| + E; + |k| < 2g. If we want to 
bring in such an inequality, the easiest way is to introduce into the 
integrand in (20.5) the factor O(2g — |E,| — E; — |k|). Equation 
(20.6) then becomes 


ee (4 =| idx) Y + (E/E) 
te AAT 3) = 1k 


—?(k,x) ae . 
[k) 


x OQg — |E| — E; — |k|)e%e 
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and (20.7) becomes 
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+ (E/|E)) 
Ya = 4 if ao 
(a) (43) oe TEE OE 
xX Og — |E| — E; — [khea" or 
Similarly for the cutoff (a) we get 
= all d°k 
Va = (4 if is a C k al : 


We have E = ay + (a, 7) — 
LS 


EP = 


So for large |k|, of the order g, 


e@o, where r, = p,; + e@,. Thus 


Amin + (a, 7) + (a, k) — e@o, 
k? + 2(k, 7) — 2(a, k)e@p + terms independent of k. 


|E| = |k| + G7) — (@, Ne@o + O(g) 
where 7, = k,/|k|. Thus 
y+ (E/JE|) _ »+ (a1) = 
Eee, % ° @ ? 
= > =e O(g a) 


when one averages over v. Also 


Cy — 2) = 195 


O(2g — 2\k| — (7) 
+ (a, De@o — E; + O(g~')) 
O(g — |k| + Foy) 


say, where 
Fo, = —40,7) + 3a, De@o — 4£; + O(g™’). 
Hence 
{O(g — |k| + Foy) - — {k\)} oe 


Let us put d°k = |k|? d|k| dQ, where dQ = 


element of solid angle. 


Then we carry through the integration first with respect to |k| and then 
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over all solid angles. We get 
1a = 1 = (4) tay | Ha, DF) AQa*. 


This quantity F,,) is, of course, an operator, but that doesn’t in- 
validate the argument. For each value of the /’s, that is to say, for 
each direction in k-space, we may take a representation for which 
Fw is diagonal; then the integration over |k| becomes trivial. If we 
use a different representation, the result must be the same. 

When we carry out the integration over all solid angles, many of the 
terms will vanish. The only terms of Fy) that contribute anything are 
those which are linear in the /’s. We get 


ay = ey = oe E ue Al 3ecto| aa 


Then carrying out the summation over the four values of » gives us 
the final result 


Yo) — Ya) = — — = tego: (21.2) 


With cutoff (c) we follow the same method; we just have to replace 
E by the kinetic part of E. Thus we replace E by E + e@o, and 
E; by |E; + e@o| in the argument of © in (21.1). Then we have to 
replace Fy) by Fce), which is Fi, = —4$(I,7) — $|E; + e@o|. That 
gives us the result 


Vics = Vay = ae eT) (ie) 


Well that is the method of calculating the differences caused by the 
different methods of cutoff. You note that the differences don’t depend 
ong. So long as g is large, these results are valid, and different kinds 
of cutoff lead to finite differences in our expression for Y. For the 
present we’re going to leave open the question of which of these cut- 
offs is the correct one, and we shall decide that later. 


§22. Case of No Static Field 


Let us go on now to take the case of no static field and work out Y 
explicitly. With @, = 0, we have 


E = (@,p)+amm, E? = p? + m’, 
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and |E| = (p*? + m*)"!? = pg, po being defined as the positive square 
root of p? + m?. So (20.6) becomes 


Is) 
l 


iS ~ i} 
4 =i) YP9 ae (a, p) + a&mm ue d k P 
Oe a rke= a SE) kl 
Let 


2 ~ a 
ea bs “1f PPo Ge (ex, i) ats AmIN 4 d k 
ee ee ee y=.) Ik 


It differs from Y in having po instead of E; in the denominator. Now we 
need to know Y only to the extent that we need Y|), which suffices to 
fix all the matrix elements Y,,;. We have pol’) = E,|7), so Yo|’) = Yi), 
and thus Yo is as good as Y for giving us the results we want. We 
shall proceed to evaluate the integral Yo instead of Y. This way of 
setting out the work is convenient for the later applications. 

We can carry out this summation over yu fairly easily. We have 


one = ==\5 
0,0; = d,, ps i,4, (22.1) 
omenyee —= ey. 


This gives us 
Y= any | =i se (GAO se Bee es 


Bo(|k| a Po — Ypa) Ik| 


And this we can write in a more convenient form 


ox i 
mae =i) / BLN aOR eee 
at) Eudes: is on 


which equals the previous one because the difference between them 
vanishes when one averages over v. 

To evaluate this integral we can proceed as follows. We choose a 
system of coordinates for which p; and pg are zero; i.e., we choose our 
k-axis to be the direction p3. Now you might say: Is such a thing 
justifiable when one takes into account that the p’s are not just numbers 
but are operators? In general that is not justifiable; but in this piece 
of work the p’s commute with everything, so we could very well work 
in the representation in which the p’s are diagonal and then take p,; and 
P2 to be zero. Later on, when we have to do corresponding calculations 
where we have field functions coming in which don’t commute with the 
p’s, the situation is more complicated and we cannot use this trick. 
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We then have 


Yo = ory | vik| — Po + as(ps + kg) + 2amm Bk 
Bo(lk| + Bo — ¥Po) k| 


We neglect terms linear in k,, kg, because upon integration over all 
k-values, terms linear in k,, kz will contribute nothing, because of the 
symmetry between positive and negative k-values. Of course, we 
cannot neglect terms that are linear in k3, because fo involves kz: 


Bo = [m? + (pa + ka)? + ka? + ke?) 
= (m? + ps? + 2psks + |kl?)"?. 
Now one can carry out the integration by making a suitable change 


of variables. The trick is to introduce a variable w = fo + |k|. First 
we pass from the three k’s to cylindrical variables. We put 


ky = p cos ¢, ky = p sin ¢, pS (ky? a7 koe 


and we pass from k;, kg, k3 tok3, p,¢. We have d*k = p dp do dk3. 
Then we pass to w, ¢, kz. Now for the second passage here we need a 
Jacobian, 


(wks) _ (*2) ees hee oe 
O(pk3) OP /k,—const Po [k| Polk 
And we then get d*k/|k| = (fo/w) dw do dk3. 


Well, carrying out that change of variables, the expression for Yo 
becomes 


YS ff ss = 20 ill = Porde v|k| + az3ks de 


@ — VPo (63) 


ie (ODD 


Let us now work with the cutoff (b), because it is the convenient 
one to use here. Then we have to integrate first with respect to kg 
and finally with respect to w. You see, w is just the total energy variable 
apart from the constant term E;, so the cutoff with respect to total 
energy means a cutoff with respect to the w variable. 

Now it is fairly straightforward integration. First, to integrate with 
respect to k3, we must know what the limits of kg are for a given w. 
For a given w, the limits of k3 occur when k3 = +|k|. We have this 
connection: (w — |k|)? = po? + 2psk3 + |k|?; that gives 2w|k| = 
w” — po” — 2p3k3. Then you see that the limits of integration are 
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kz = (@” — po”)/2(w + ps) and —(w? — po”)/2(w — ps3). Now if 
you look at the expression (22.2) for Yo, you see that there are some 
terms independent of k3, some linear in kg, and some linear in |k|; 
and so we need three integrals, 


/ pos Pp 2 
d So ee 0 
k3 AS p32 W, 


| ee me (w* — po’) 
Jradts = 3 GB hare, fists = $e Oe 


This gives us 


ly ae [ests ate 2am = (26 aye y (w” = Po dw _— os P3(w me Po’) 
a? — ps? 2a? — pa? 2w? — pa®? 


X (@ + vpo) dw. 


At this stage it is convenient to do the v averaging. We could have 
done this v-averaging at any time during the calculation, but it would 
be inconvenient to do it earlier; it would mean our terms would get 
longer. But here we can do it conveniently. We find 


-. a3p3 + 20mm — po , (Po — asp3s)(w” — el 
tore /| Oe jie = P(E — TURAN ee 


Now the w? — po inside the [ ]canbe replaced by w? — p3* — m? 
and we get 


Fe 1 files Be Get — 70 ee Poe eb | Pie 


2 OO ike (Go — jee 


3 
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This is something we can quite easily integrate with respect to w and 
the result is, when we put in the limits of integration, 


< 2g9—E 
WV, = | apatite i= polyeGe = eee ee es m| 
4 w P3 Po 


= $(a3p3 + 4amim — po) log (2g/m) — a(Po — a3Ps). 
We did this calculation for the case of p; and p2 vanishing, but we 
can go back to general p’s. For general p;, p2, p3 Obviously we just 


have to replace a3)3 by the scalar product of (e, p). We get the general 
expression 


Yo = 4[(@, p) + 4anm — po] log (2g/m) — 4lpo — fe, pl (22.3) 
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Now we have [(e, p) + am]|i) = E;|i) = poli). We then find that 


Yi) = anm[Z log (2g/m) — 4]li) = gmam|i), 


where q = 2 log (2g/m) — ¢. 
That is the result of this calculation. Let’s see what effect it has. 
From (19.2), 


Oke. 
aor 


2 . 
Nes qm ys Wale |e rama a 
WT 
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This shows us the way in which the perturbation affects the emission 
operator we’re interested in. 


§23. Mass Renormalization 


This effect can be eliminated by a procedure known as mass re- 
normalization. This procedure consists of adding a further term to the 
Hamiltonian corresponding to a variation in the mass parameter. Let 
us keep m to be the observed mass. In H replace m by m + 6m; that 
is, the mass parameter as it occurs in the Hamiltonian is now going to 
be different from the observed mass. This variation in the Hamiltonian 
is going to be treated as a second-order perturbation. We get an extra 
term in H of amount 


Hy = im | Toa dom Se Pn (nlom|n' Wn 


nn! 


This is to be treated as an extra second-order perturbation, giving an 
extra variation in K*, 


oe") ee 
l = [H. »>WVil—> 
( or extra [ 2 ie ] 


Ha" = 6m D1 Tn (nlam|e! ret Pn ¥n 
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where 


Using this H2* in the commutator, we get 


* 
,OKe = Ye ( | \ (EH, —E)t 
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This extra i(0K»*/dt) is to be added to the previous one. 

Now we can make the total i(@K2*/dt) = 0, provided we choose 
dm suitably, that is, provided we take 6m = (—e?/x)qm. This pro- 
cedure is mass renormalization. The result is that the perturbation 
does not effect the electron emission operator to the second order of 
accuracy. 

At this stage we might discuss what value we should take for g; 
g must not be too large because we have to be able to treat H»* as a 
perturbation; that stops us from making g go to infinity. But we can 
take g to be fairly large and still keep 6m/m rather small, because of 
the logarithmic factor. In our units, we have e? = 1/137. I worked 
out a few figures. With 2g¢/m = 1000 we get log (2¢/m) = 6.9 and 
6m/m = 1/43, which is small enough to be reasonably counted as a 
perturbationeffect. We might take 2g/m = 10.000; then log (2g/m) = 9.2 
and 6m/m = 1/32, which is still reasonably small: for a perturbation 
theory. g is of the order of 250 X 10° eV in the first case and 10 times 
as big in the second. One must expect electrodynamics to break down 
with energies somewhere between these two orders because there will 
be further effects coming in at these high energies due to the possibility 
of creating mesons, neutrinos, etc. So we should imagine 2g/m to 
have a value somewhere in this range, which means that our electro- 
dynamics is limited from the point of view of high-energy processes. 
But we may expect it to be a pretty good theory for processes with 
energies small compared with these energies. 

With this cutoff, 67/m can reasonably be counted as a perturbation 
effect. That I think is how the mass renormalization ought to be 
handled. Handled in this way we do get a logical theory where we’re 
solving equations, not merely using working rules, and we’re neglecting 
only terms we can feel are small and for which the neglect is justifiable. 

Now this calculation is done with a cutoff (b). How will things go 
if we use a different cutoff? You see that we did use cutoff.(b), because 
our last variable of integration was w, equal to |k| + fo, and the 
cutoff of this variable means a limitation of |k| + fo. That just cor- 
responds to what we have with the cutoff (b). 

With cutoff (a) we get, according to (21.2), Y = qamm + §(q, Pp). 
Now we can get rid of the first term by mass renormalization but we 
can’t do anything about (a, p). You might say: Can we not introduce 
a further term in the Hamiltonian, a sort of counterterm to (a, p) 
which will cut it out? We could introduce such a further term, but that 
would have to be of the form of some multiple of (@, p), and such a 
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further term is not of a relativistic nature. If we put in such a counter- 
term in the Hamiltonian to get rid of (a, p) our theory would no 
longer be relativistic even for small-energy processes. Well, we infer 
that cutoff (a) is a bad cutoff. 

At present we cannot distinguish between cutoff (b) and cutoff (), 
because they are equivalent for the case of no static field. But we can 
distinguish between them by an argument depending on gauge in- 
variance. Let us suppose the static magnetic potential is zero, @, = 0, 
and the static electric potential is constant in space as well as in time, 
@ = a c-number. Now such a potential, of course, means no field 
at all. It ought not to affect things at all. Let us bring such a potential 
into our equations and see what effect it has on them. This means 
working with an energy operator 


E= (a, p) + amm — hoy = ie eho, 


where U is the energy operator in the case of no potentials. 
Now in evaluating (20.6) we are working with the quantity 


De CeVa) 
Ik] + (Z| — vEy (23.1) 


multiplying it by suitable factors on the left and the right, and then 
integrating; and we must see how the result is affected by bringing 
in the external potential. Let us make the restriction |e@o| < m for 
simplicity. This means that the potential must not be a terribly 
big thing; but it still can be fairly big. Any eigenfunction of U is an 
eigenfunction of E, because the difference between U and E is just a 
c-number. If U> 0, then E> 0, in virtue of the inequality 
|e@o| < m; therefore, E/|E| = U/|U| = U/po, where po, you remem- 
ber, is just (p? + m?)!/?_ - 

First, let us replace |E| by vE in the denominator of (23.1); we get 
[fy + (E/|EDVUk| + v(E — E)]. The replacement is permissible, 
because |E| = +vE and if |E| = —vE we have the numerator 
vanishing, so that the equation is valid in any case. Let us now write 
it [vy + (U/po))/[|k| + v(U — e@o — E;)]. Now we can make the 
same sort of change again, replacing vU by |U| = po, using the same 
argument again. We thus get [vy + (U/po)]/[|k| + po — v(e@o + Ej). 

The result we have here is the same as in the case where there is no 
static potential, except that the number E; is replaced by E; + e@o. 
In the evaluation of the integral you see that everything would go as 


24. The Anomalous Magnetic Moment Wi 


before, with a number E£; + e@p replacing the previous number Jee 
So our previous calculations for the case of no potentials can be taken 
over with this substitution. 

When we previously carried out the calculation, we made use of the 
weak equation (@, p) + amin ~ E;, valid with both sides applied to 
|). Now we have instead 


(a, P) So Amn =~ ic; - eo, 


because we still have E ~ E;. So here again we have to substitute 
the number E; + e@o for E;. The whole calculation goes through as 
before with this substitution, and we get the final result that Y =~ gan, 
with the same q as before. With mass renormalization we get 
dK»*/dt = 0, as before. 

You see everything goes through all right with the substitution. We 
ought to see what kind of limiting process we are now using. Previously 
we were using a variable of integration defined as w = fy + |kl, 
and we did the integration with respect to w last. That means we had 
a cutoff with respect to w. Now this variable w is unaffected by the 
substitution of E; + e@o for E;, so we are still working with a cutoff 
with respect to w. The variable w just refers to the total kinetic energy, 
so we are working with cutoff (c). 

Cutoff (c) is the physically correct cutoff which gives a gauge- 
invariant result. Cutoff (b) is not gauge-invariant. 

We did this calculation using an electrostatic potential constant in 
space and time. You might say: What would happen if we used a 
magnetic potential that is constant in space and time? You easily 
see that that has quite a trivial effect on the calculation; it just means 
replacing p, by z,. After the integration p, disappeared from the final 
result, so the z, will disappear equally well. So a constant magnetic 
potential has no influence on the calculation. Well, that completes 
the theory for the case of no static field or the case of constant potentials. 


§24. The Anomalous Magnetic Moment 


We now go on to consider cases in which these potentials are not 
constant. Let us first consider the problem of the anomalous magnetic 
moment. For that purpose we have to take the potentials @, to represent 
a uniform magnetic field, uniform and constant over the whole of 
space, and take @y = 0. 
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Let us suppose the magnetic field is of magnitude % in the direction 
of x3. Now we have 


E = (a,p + e@) + anm = (a, 7) + amm, 


where z, as before, is p + e@. From the one-particle theory we get 
E? = x? + m? + e3Co3, where ajag = i¢3. 
We shall assume 3¢ is small and neglect 3¢?. Thus 


|E] = (mn? + 2?)? 4 desta gm? + 37)”. 


We have 
oO 
Y = (40) | oye! i’ 
with 
¥= (r+ a) (E| + [k| — vB". (24.1) 


We shall suppose the momentum of the electron is small, so the 
as are small, and neglect terms of the third degree in the 7’s. We 
shall also neglect 73C. 

The denominator in Y sometimes vanishes ify = 1. It cannot vanish 
for y = —1. We must work out the order of magnitude of the con- 
tribution to Y of the domain of integration for which the denominator 
is small. We shall find that it is negligible. 

We have 


(Ee aak| 2 
(\E| = |k| poe (Be ee et 2k ee 
(PSRs 129) 2 Ai) <> fe Bm 


Write 
Ds 4(E? = am) + (kr) a Ik|E;. 


We can work with this denominator instead of the one in (24.1) with 
y=1. 
We need to know YF only to the extent that we need Y|i). We shall 
write 
AzB if Ali) = Bii). 


Let a, denote a longitudinal a with respect to the direction of the 
magnetic field, a, = 1 or a3. Let a, denote a transverse a with respect’ 
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to the direction of the field,a, = a; or ag. Then a, commutes with E2, 


oE*a, = E?, 
but 
el = on? — eRo3 = E? — 2e8a3. 


In these equations there is no summation for } or 7. 
We have a,Day = (kr) + |k|E;. Thus a,Da, cannot vanish, 
because of the smallness of (kr). We have 


a,Da, ~ (kr) + |k|E; — estos, (24.2) 


and this can vanish for small values of |k| of order 3¢. Its vanishing 
corresponds physically to the electron emitting or absorbing a photon 
and changing its spin direction. 

For the transverse a, , the numerator in a,Ya, (with y = 1) is small. 
It is (apart from a factor of order unity) 


ar {1 + [(ar) + (ak) + amm)|E|}0,, 


which is of order (1 — am) + terms of order 7 or k, and (1 — am) 
applied to |) is also of order 7. 

The domain of integration for which the denominator (24.2) is small 
can be pictured as a thick eggshell enclosing the origin, with a radius 
of order e3C/m. We get the order of magnitude of the integral over 
this domain if we replace the thick eggshell by a thin one, putting 


(kr) + |k|E; — e3o3] for [(kr) + |k|E; — e3Ca3]}—'. 


We then get for the integral 


3 
/ [O(k) + O(n] (km) + |k|E; — estos] WT = OK) + Ofkr). 


Here we must substitute for k its value in the shell, which is of order 3c. 
The result is negligible. 
Define 


Yo = + [(at) + an] 
x (m? et my V Cm? = myil2 4 Ik| i: v(m? a a) le", 
ile, = ae ae a, ela 


= a d°k 
Yo = (41) | aS%00 1K : 
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Suppose the 7’s are small and expand GY and Yo as power Series in 
the x’s as far as terms of the second degree. The only second-degree 
terms in Yo involve (ax)(km) or (kx)? or w*. On integrating over all 
k-directions (ar)(km) gives zero and (km)? gives 3k?m7. So the terms of 
second degree in Yo all involve 77. Hence to this accuracy Yo is the 
same function of the z’s whether the 7’s commute or not. 

If we look upon the z’s as commuting, UY is the same as the 
corresponding expression of Section 22 with 7, for p,, and so Yo is 
given by the Yo of (22.3) with z, for p, which gives 


Yo = Alan) + amm — m — (x*/2m)][log (2g/m) + 3] 
+ amm[z log 2g/m) — 4]. (24.3) 


This is obtained with the cutoff (b) or (c), because the last variable of 
integration is w = |k| + (m? + #?)1/?, which differs by a negligible 
amount from |k| + |E|. 

Let us write e3Co3 = h for brevity. We have from (24.1), 


= {vy + [(am) + anm\(m? + #? + Ay?) 
X (Gm? + 2? + AY? + [kl — vE I’. 


Yo is the same expression with h = 0 and (m? + x”)"? for E;. 
We shall expand the difference Y — Yo in terms of the small quan- 
tities A and E; — (m? + 71”)?. The latter is weakly small, since 


E; fue (m? at ny il? pa |E| ooh (m? ali mele = 4h(m? = ie 
(24.4) 


In the expansion we do not have E; — (m? + x7)'/? occurring on the 
extreme right, so we cannot immediately use the expression (24.4) for 
it. We get, ignoring some factors on the left, 


[Es — Gm? + 97)" ? Gn? + #7)? + [kl — vOn? + 2?)?777 
This differs from the expression 
[@m? + ?)N? + [kl] — v@n? + 9?) ?71LE; — (mn? + 2?)!7), 


to which (24.4) is applicable, by quantities involving [(rk), x], and 
such quantities are of order 73¢ and are negligible. Thus the expansion 
is justified. 

In the expansion of Y— Yo with the use of the weak Sgqutetste 
(24.4), every term involves / or 3. Thus we can neglect the z’s alto- 
gether, because we are neglecting 73¢. 


Put 
Then 
Y-— Yo = 
We have 
Thus 
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(m? al yur = i, 


—amin(h/2p*)\(u a [k| aad vn)" — V0 mIN/[L) 
x (ue + [k| — vin)~*[(h/2u) — v(E; — m)). 


mn 


i = 
a,a" = —2, Ayame’ = dap, 


a,eza* = 203, AyAmT 3a" = 0. 


an(Y — Yojat = —(u + |k| — vm)? [v(h/u) + 2(E; — m) 


Now 


Ce | — 


Squaring, 


—4Vam(m/p)(E; — m)] 
= —(u + |k| — vin)~*[v(A/u) + (h/m) — 2v(h/p)] 
= h(u + |k| — »m)~?@m — p)/mu. 


vm)(u — |k| + vm) = wp? — k? — m? 4+ 2v|k\m 
2v|k|m. 


(uw + |k| — vm)?[u? + k? + m? + 2vm(u — |k|) — 2u|k|] = 442m? 


or 


(u + |k| — vm)?(u + vm)\(u — |kl) = 2k?m?. 


Thus 


a(Y — Yola* ~ h(u + vn)\(u — |k\)m — p)/2mPk*y 


= —h(u — |k\)/2mp, 


Y — Yo = —(h/2m*)\(40)7! i (1 = |k\/u) a°k/[k| 


—(h/2m*) / (1 — [k|/u)Ik| dk. 


It does not matter which cutoff we use for evaluating this small 
quantity. Use cutoff (a). Then we get for the integral, 


HK? — km? + KP)? + m? log (e+ Qn? + 4?)"?f} 09 


Thus 


= $m’*[log (2g/m) — $). 


Y— Yo = —(h/4m)[log (2g/m) — 4]. 
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Taking the expression for Yo given by (24.3) and omitting the last 
term, which is canceled by mass renormalization, we get 


~ Alam) + amm — m — (r* + h)/2m] 
X [log (2g/m) + 3] + h/4m. 
Now 
rl =p Ga — he (x? + A)/2m = E — |E| = 0. 
Thus 
Y = h/4m = eXo3/4m. (24.5) 


This corresponds to an extra energy 
(e?/n)¥;; = e?o3/4am. (24.6) 


The elementary theory gives for the energy in the magnetic field 
eKeo3/2m. So the magnetic moment is increased by the factor 
(1 + e?/2n), or, bringing in the constants, (1 + e?/2rhc). This value 
for the anomalous magnetic moment of the electron is in good agree- 
ment with experiment. 

Formula (24.5) contains a rather stronger result than we used in 
(24.6)—that Y doesn’t have any nondiagonal elements. That means 
that the eigenfunctions are undisturbed by the interaction to the 
accuracy to which we have worked things out. The interaction in the 
case of a static field consisting of a uniform magnetic field is thus a 
very smooth thing which doesn’t perturb the eigenfunctions of the 
energy. Its only effect is to give an extra energy term corresponding to 
the anomalous magnetic moment. 

If you want to compare the theory I have given here with the usual 
theory working in terms of the Schrédinger picture, I think the best 
book to use is W. Heitler, The Quantum Theory of Radiation, 3rd ed., 
Oxford University Press, London and New York, 1954. If you make the 
comparison, you will find that the old work is not really a logical de- 
duction such as we’ve been giving here, but still you have the same 
sort of integrals turning up, and the same substitutions have been 
made for evaluating those integrals. There’s a great deal of similarity 
in the details of the calculations, but the argument of course is very 
different. I think the argument is a poor one when one tries to 
work in terms of the Schrédinger picture. 

Calculations of the anomalous magnetic moment have been made 
to the fourth order, that is, to e*/h?c”, by R. Karpus and N. M. Kroll, 
Phys. Rev. T7, 536 (1950). These calculations are extremely heavy. 
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The Karpus and Kroll calculations contain some numerical mistakes. 
They have been amended by A. Sommerfeld, Phys. Rev. 107, 328 


(1957), and Ann. Phys. 5, 26 (1958), and also by Petermann, Helv. 
Phys. Acta, 30, 407. (1957). 


§25. The Lamb Shift: Preliminary 


I would like now to go over to the other main application of quantum 
electrodynamics to the Lamb shift. We now have a static electric 
field but no static magnetic field; that means we have a potential @p 
which can be any function of x;, x2, x3, @p = @o(x1, X2, Xa), but 
is not a function of time, and the static magnetic field is now zero, 
Ga— Uv. 

That gives us for the one-electron energy problem E = (ae, p) + 
a,m — eQ@, where @ is written for @y, because we’ve only got one 
potential coming in this work and there’s no need to write the suffix 
all the way through. I write, for short, E = U — e@, where U stands 
for the operator (a, p) + amin, so U? = po’. 

We get E? = po” — e(U@ + GU), neglecting @7. @ is a general 
function of x,;, x2, X3, so of course it does not commute with U; 
that is why we have to write these two terms U@ and @U separately. 
We’ve now got to find the square root of E” to get |E|. How do we 
find the square root of a general operator like E?? We can do it when 
we take U@ + QU small, and the result is 


|E| = po — 3(UB + BU), (25) 


with a suitable quantity B of order @. One can find out just what 
this quantity B is by squaring |E| and comparing with E”. If one does 
this, one sees that one needs 


PoB + Bpo = 2eQ. (25.2) 

Let us resolve @ into Fourier components: 
Co — ‘| Grem ee dia (253) 
@,, of course, is just a number. I ought to point out that this 


Fourier resolution is essentially different from the Fourier resolution 
that we made for the electromagnetic field variables. Here there is no 
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variation with time; it is a straightforward Fourier resolution in three 
dimensions. From (20.3) we find that 


(p'|@|p”) = ic b3( pp” kya hae aes 


Now let us take (25.2) and express it in this matrix representation; 

we get 
(pol + po'’)(p'|Blp") = 2e(p'|@|p") = 2e@p—p’, 

where po’ = (m2 + p’?)"2, po” = (m + p’*)"2. Well, that fixes B: 
DXA ON eer ; 
Po + Po" 

Now from (25.1) |E|~* = (po)~' + (2p0) (UB + BU)(po)—’. 
We have used the general formula (x + n)~! = x7! — x7!nx7! for 
expanding a reciprocal with n small, when we have noncommuting 
quantities coming in. We now form E/|E|; E and |E| commute with 


each other, so we don’t need to specify an order for the product. If 
we evaluate the product with £ on the left, we get 


E/\E| = (U — e@)(po)~* + (2po) (po? B + UBU)(po)~", 


(p'|Blp”) = (25.4) 


with neglect of second-order quantities. Using (25.2), we get 


EO. UO Boo ee eee ee 
par Ve Ren is? 25.5 
[Z| po 9 2" 296 pops po po ee) 


where JN is defined to be 


N = 4(poBpo — UBU). (25.6) 


§26. The Lamb Shift: Vacuum-Polarization Term 


J should now like to proceed to calculate the vacuum-polarization 
term. This term gave no contribution to the anomalous magnetic 
moment, but it does give a contribution to the Lamb shift and we have 
to work it out. 

The formula for the vacuum-polarization term is (20.1). We only 
want the diagonal elements of Z for calculating an energy change, and 
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I shall confine myself to working out the diagonal elements because 
it’s too complicated otherwise. They are 


dike 
||? 


Zi = x} i Ue Glen |) (26.1) 

J," is given by (20.2). We can evaluate the diagonal sum in any 
representation, and of course we now choose the momentum repre- 
sentation, as we already have our quantities A and B expressed in 
terms of this representation. We get 


mt} tk, x)| 7 , E " , " 
t= aff prietornnrot (1 — EY yp dora 


Here we have taken the diagonal sum as far as it concerns momentum 
variables. We also have to take the diagonal sum with respect to spin 
variables; this is indicated by (e“{1 — (E/|E|)]). Using (20.3) and 
(25.4), we get 


U l | 
SA aff in Hi 7 me y f 
k 4 (p'|\@ Po ean he 
~ 5(p’ =, ae k) aa don 


The first two terms in the [ ] contribute nothing, because we take k 
to be different from zero. We are not interested in the case of k = 0; 
that would just correspond to a constant potential, which we’ve already 
dealt with in our discussion of cutoff (c); we need to take into account 
only those terms of [ ] which have a nonvanishing Fourier coefficient 
fork # 0. 

Now let us go back to the expression (25.6) for N. It involves only 
even powers Of a1, a2, @3, Am. It follows that (a,V) = O(r = 1, 2, 3), 
because a@,N will contain only odd powers of the a’s, and so has no 
diagonal sum. Hence J,” = 0. We are left only with the problem of 
evaluating J;,°: 


Ie = 3 / / (p'l(Pp0)(N)(Po)'1p") 8p! =p" + k) dp! dp" 


Now we need to evaluate (4). It is only those terms in N that don’t 
contain any @ matrices at all which contribute to the diagonal sum. 
So we have 


4(N) = 4(poBpo — prBpr — mB). 
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Substituting in /;,” and using (25.4), we have 
0 , ” CBp''_p! 
— 6 — ie 
J // 3(p P + ae, aie Don 
HOBO Er or eae 
Po Po 
It is convenient to go to new variables of integration. Put 


x 


oD a = ae ale. 
so that p’” — p’ = k. With this change we get 


Fe ff Pee ye) 2 
ie = can | a are eee "Bp = eQh (26.2) 
say. 

If you compare this with Heitler’s book, /;, is Heitler’s L44,f with 
ky = 0. Using quite different arguments in terms of the Schrodinger 
picture, one turns up with the same integral in the end. The integral 
T;, is the thing that determines the vacuum-polarization contribution 
to the Lamb shift. 

We have to evaluate the integral /, and we can do it by making 
suitable substitutions. First of all, k is a constant as far as concerns 
this integration, so we can choose k in a special direction. Take k 
in the direction of the axis x3, kK = 0,0, k3, to simplify things; and 
put p: = pcos ¢, po = psindg. Thus 


Po’ = (p? + ps? + m? + dk3? — paks)"?, 
Po’ = (p° 4 ps? + om? Fah = ska) 


Now we are going to introduce a new variable by putting 
Po + Po” = o. 


You will note some analogy in this work to what we did in evaluating 
the integral in the case of the free electron. There again we introduced 
a variable w which was the total energy. 

We are going to pass to p, y, w as variables of integration. We have 
d°p = pdp dd dp3, and we have to work out the Jacobian 


d(, P3) _ (2) _ pw 
a(p, Ps) op p3=const PoP 0” 
{ W. Heitler, The Quantum Theory of Radiation, 3rd ed., Oxford University 
Press, London and New York, 1954, p. 322, eq. (16). 
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The result is that d°p/po'po" = (dw/w)dédp3. With this change of 
variables, we have 


i = 2m | (po' po! =p eK mya? des dps 


/ ‘i tA vt “a d. 
a [po po” — 4(po”” + po’) + 3k3") dps 


d 
= rf Us" — (po" — po!) 1-5 aps. 


Now we note that 


m2 


(Po Poe = po = pa” = 2pskz, 


and this enables us to write 


Ae Ved 
= aka? |{ (1 - Ba) Be apy. 


Let us carry out the integration first with respect to p3. We need to 
know the limits of integration for each value of w. The limits of ps3 
occur when p = 0, and then 


w = [(p3 + 3ks)? + m7]? + [(p3 — $k3)?2 + m?]}4/? 
and 
gw? = p3* + 4k3? + m? + [(p3 + 4k3)? + m?] 
X (ps3 — 3k3)? + m?]"2, 
sO 


(—4w? + pa? + 4k3? +m)? = (p32 + $k3? + m?)? — p3?k3?. 


This reduces to the following: 


oan ae 2 Ane = ka? 
Ps Aw? = ke?) 


and plus and minus the square root of this quantity gives us the limita- 
tions on p3 for each value of w. 
We now find with those limits of integration, 


iE a... ane a 4m? = ks” yi 
P3 Tee 


2 2 Q\ 3/2 
QQ = 4m = kz 
|v! dps iso ( oe k32 ) : 


and 
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‘5 Be es 4m? = fee 1/2 1] Be — 4m? = ke ee dw 
Ty ates Ul = ee ee E 


At this stage in the work we can go back to the general & and re- 
place k3” by k?. To carry the integration further we put w* — k? = z?, 


and then /;, goes over to 


_ ank® f (2? + 2m® yz? — 4°)? 


3 z2(z2 + k?) see 


Ty, 


A further substitution is needed to get something one can integrate: 
1 — (4m?/z?) = y”. Then, with a little bit of reduction, we get 


serie cere y?dy | 
oo 3 4m? + k? — k2y2 1 — y? 


which can be integrated by standard methods. 

Now what are the limits of integration of our variables. Well, w? 
goes from 4m? + k? to 4g7, so 4m? < z” < 4g?, and y goes from 0 
to 1 — (m?/g”). We take the upper limit of w to be 2g, because we’re 
using cutoff (c). Our procedure of using the variable of integration w 
is really very suitable for working with cutoff (c), because w is just the 
total energy of the intermediate state, which goes up to 2g if we neglect 
|A| with respect to g. 

Carrying out the integration, we get the result 


2 
fy = 
cal 
¥ 


we 

x 27 log =8 — bak Cl — 2/( ) log [v7 + (y+ nt] 
where Y = k*/4m?. Now we have this rather awkward expression to 
deal with. Let us make the assumption that Y is small. You see what 
that means; it means that our electrostatic field does not have any 
components with very short wavelengths. That is quite reasonable 
when we apply it to the field in the hydrogen atom, for instance. Any 
wavelengths there will be long compared with m—', so the approxima- 
tion is quite good. 

We need to evaluate /;, to the accuracy of Y” and the second log to 
the accuracy y*/?. This log becomes y¥/? ~ dy3/2 4 3-5/2 and 
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the whole expression for J, becomes 


4 kA 
hk : y~ $7] 2 wk 
:) 5 as \5m?2 


a 2a 2g > 
a (ioe m =) 
Equation (26.2) becomes 
eae gi 
Lp” = (ak x) ey. 
Equation (26.1) becomes 
oe a nies ak? Qalile@™ |) ak 
ae oe 15m2) ~* 
(e F —1(k,x 2 (é . —tk,x a 
£ 44 | axe OO) OV ee gra Gif ane ed ile 


= qlil@s|i) + 7 UlV eli), (26.3) 


lI 
BB, 
3 
4 
(accom 
a 
nN 
lay 
i 
men 


with 


with the help of (25.3). Z;; multiplied by e?/z gives a contribution to 
the energy. 

Now you see that we have an infinity coming into the first term if 
we allow g to go to infinity, because we have log g coming into q. 
We can, however, get rid of this term by a procedure known as charge 
renormalization. It consists in assuming that the parameter e, which 
occurs in the Hamiltonian, is not the observed charge. We keep e to 
denote the physically observed charge and suppose that the charge 
parameter in the Hamiltonian, which we call e + 6e, differs from the 
physically observed charge. 

Now the Hamiltonian consists of two parts Hy + Hg. We are 
going to take 6e of order e®, so that adding the ée contribution to Hg 
would not affect things to the order of accuracy to which we are working. 
We are only working to the order of accuracy of e” as far as the pertur- 
bation is concerned. But e also occurs in Hp and 6e makes a significant 
change in Hy, 6H y = — bef apy d°k. Let us look upon this change 
as a further perturbation. We may write it 


5H) = —ée ye Yn (nl@ln Wr. 


nn! 
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Making the transformation to the starred representation, 


bHo” = —8e >) dn(nl ala’ ner. 
nn’ 


This gives 


(2) 
ot extra 


The term here for which n = j also corresponds to a change in energy, 
and we get cancelation with the energy change coming from the first 
term of (26.3) by suitably choosing de: 6e = (e?/1*)gq. 

By changing the Hamiltonian in this way, a way that does not really 
have any physical significance and means simply that the charge 
parameter in the Hamiltonian differs from the observed charge, we are 
left with just one term in Z which gives a contribution to the Lamb 
shift. 

Question: Must we get rid of that first term in Z,;; because it depends 
on g? Is that the reason? 

Answer: I think that the answer to your question should be no. 
Even if g did not involve anything that became infinite with an infinite 
cutoff, charge renormalization would still be the proper thing to do. 
Otherwise, the energy associated with a charge e in a static field @o 
would not be simply e@q; there would be an extra term and the energy 
would be (e + 6e)@5. Now one usually takes the definition of the 
observed charge to be such that, if multiplied by @o, it gives the energy 
of the electron in a field @y. If one is going to work with that definition 
of the observed charge, one is forced to the conclusion that the observed 
charge is different from the e in the Hamiltonian. 


[5Ho, vi] = —Se ~ P,(n|@iePn—# 0! 


§27. The Lamb Shift: Relativistic Term 


I now want to calculate the other contribution to the Lamb shift, 
which comes from the integral Y: 


3 Os py 
y= any | a moth, @ — v+ (E/E) 
ue cote yg 2s a) oe 


Let Yo be the corresponding operator with no external field @ and 
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with po written for E;, so 


Wg = v + (U/po) =e 


Pao [k| = PDO 


(so) 


Po 1s by definition positive. 

Now Yo gives an integral Yo which has already been worked out 
when we dealt with the case of no external field. The result was 
(22.3). When we bring in the usual mass renormalization we get an 
extra term 


— alls log (2g/m) = Z]. 
Adding this to (22.3), we get 
Yo = 31, p) + amin — polllog (2g/m) + 4}. 


Now we make use of weak equations that are appropriate to our 
present work. We have E/|E| ~ 1, which merely says we are starting 
off with a positive-energy state, and that leads to, from (25.5), 


U — po = N(po)’. (272) 
With U = (a, p) + a,m, we find that 
Yo ~ $(log (2g/m) + 5)N(po) |. 


I shall use the notation F ~ G to mean (i|F|i) = (i|G\i). We only 
want to know the average value of the operator Y for the state |/). That 
average value is going to give us the Lamb shift. So we use this notation 
for things with the same average. 

We have from (25.6), 


(i|2N(po) |i) = (i\(poB = UBU)(po) ' |i) = (Po = 0B 
= —(i\(po) 'NBii) 


from (27.2). Now, NB is of order @? and is negligible. So 


t), 


we (27.3) 


Po 


and Yo gives no contribution to the Lamb shift. 

We must now evaluate (47) 'fa(Y — Yoo"(d*k/|k|). We have 
got rid of the main part of the integral Yq, and are left with quantities 
of the same order of magnitude as @. 
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eet 
D = (denominator in Yo) = Bo + |k| — po. 
Then 
(denominator in Y) = D — 4(0B + BU) + v(po — £)). 


Now (den. in % — D is weakly of order @. If we make an expansion 

counting (den. in Y) — D as small we get 

G= (++ G) (Ot + DOB + BO) ~ v(p0 ~ BID. 
(27.4) 


I would like you to note that it is permissible to use the weak condition 
Po — E; = O(@) here, because (po — £;) commutes with D. 

The expression (27.4) for ¥ can be criticized because the denominator 
sometimes vanishes. Now we had corresponding trouble in the case of 
the calculation of the anomalous magnetic moment, but there the 
trouble was not serious, because it turned out that the region of k-space 
for which the denominator was small, and for which therefore the 
expansion was not valid, was such a small region that it did not affect 
the final result. However, for our present calculation that is no longer 
true. There is a region of k-space for which the expansion (27.4) is 
not valid. That region is too large for us to be able to use the expansion 
throughout the whole of k-space without getting into serious error. 
Under those conditions there is no method of approximation that is 
valid over the whole of k-space. What we have to do is split up our 
domain of k-space into two parts, the high-energy part and the low- 
energy part, and use the above expansion only for the high-energy 
part, in which it is legitimate, and use a different method of approxima- 
tion for the low-energy part. It is a complicated procedure but it is the 
only way. 

We separate the two domains of k by the condition that |k| > ym 
for the high-energy domain and |k| < Ym for the low-energy domain; 
and we have to choose Y suitably. Choose ¥ so that Y < 1, and at the 
same time Ym >> |e@|. In the case of the hydrogen atom, where 
|e@| is of the order of the binding energy, it is possible to choose Y 
in that way. Ym has to be some energy between a few volts, the binding 
energy of the electron, and a half-million volts, corresponding to the 
rest mass of the electron. The precise value we give to Y does not 
matter, provided these inequalities are fulfilled. In this way we get 
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two terms, which are called the relativistic ( YR) and the nonrelativistic 
(¥np) contributions to the Lamb shift. 

Let us proceed to calculate Yr, using the expansion (27.4) and using 
(25.5) for E/|E|. We have the equation 


3(UB + BU) = a,mB + 4a, pB + Bp), 
and from (25.2) and (25.6), 


LU 


i = = = 
5 oe + BU) = e@ — (po) N. 


These equations give us 


~ 


Y— Yo = Yi + Yo + Vs, 


where 
ag v mK Eee 5 v(@, p 
Vi = Zp (FB + BD — MP (p. — By, 075) 
Vo = —(o) (ON + No) Dy 0276) 


Amn 


Po 


Ys = 2 | an ~ (> + ) ~ Ej) + vee | De Or) 


I shall assume that the momentum of the electron is small and retain 
terms of order of p? but neglect terms of order p*. We have to go as 
far as the p? terms to get anything interesting. I shall let 


p= (kK? + m?)\N?, 


This is a different » from what we had in the magnetic-moment case. 
Then 


ik ,k)? 2 
Bo = [u’ + 2(p,k) + pe}? = a + He . Sr + FI 
(27.8) 
_ s Tie sep, ke z 
ect [fp 222 ce ~ 2 (27.9) 


We are now going to work with the cutoff (a), because the difference 
between the cutoffs (a) and (c) is negligible. It is important to use the 
right cutoff when calculating Y,, but the difference between Y and Yo 
is small, and for this difference it does not matter which cutoff we use. 
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Now let us calculate Y,. You should note that in the expansion of 
D~' and in B we can now neglect the p” terms. We need retain only 
the terms linear in p. The reason for that is that every term in Y, con- 
tains (a, p), i.e., (a, p) or (a, k). Now if a term already contains p, 
any factors in that term need be evaluated only to the order of p, 
because we want the result to the accuracy of p®. Also, if the term 
contains (a, k) and is multiplied by a p? factor, it will give zero when 
we average it over all k directions. 

We have from (25.4), 

ALP vO 2€Qp'—p! F 
(p'|Blp ) = Da +E Bo’ 
where 


Po = bro + (pk) = ha i 


Thus we have from (27.9), 


(p'| Bl p’’) = CQ me ime [ a oye) . 
2p? 
Also we have 


e.0 k) 


+ [kl — ym = + @0, 
mv 


D3) Kk) pg eet 
where \ = wp + |k| — vm. So 
fel (oh i = (P, w]. 
TON 


There are the two terms for UY, in (27.5), which I shall deal with 
separately, say Yy = U1, + Yyo. We have 


CR _pr 


(P'|Giilp”) = ar (a, p’ + p” + 2k) 


x f est 0]. 
Qu? 


which reduces after averaging over k directions to 


WO “= see " k? 2k? 
(PY iil’) aw + v)(1 - - aa - 34): 


We may write this result in symbolic form 


~ 


ve ke Dk? 
= Faas : a —- —-—- — ——]}. 
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Tie = Seep + )(1 22H) _ OW), — py 


ee eo lee 
~ ie (a, p) ¢ Rice *) (Bara). 


We add these two terms together and carry out the integration and the 
v-averaging. We get 


g 
= = = d*k 
og (47) : i on(Yi4 ae Yr2)a" Tk} 


m™" (log 3 ~ 3) Ue PIE: ~ po) + eles va + @p)} 


We can reduce this somewhat by considering expressions that have 
the same average. First, we see that (a#, p@ ~ (FE; — a,»m)a@ with 
neglect of second-order terms in @. This has the same average for 
the state |/) as @(E; — amm) ~ Q(a@, p). Thus we have 


(a, p(E; — po) + zel(a, pe + a, p)] 
~ (a, p(Ei — po + e@) 
~ (@, p(U — po) ~ (Ey — amm)(U — po) ~ (Ei — anm)N(po)~* 
~ —ammN(po)—', 

where we have used (27.2) and (27.3) in the last two steps. We get the 


final result 
Y, ~ —log [(2v)~! — 4JamN(po)'. (27.10) 


Let us consider Y next. We have from (25.6), 
= 4[p0Bpo — (a, p)Bla, p) — m°B — (@, p)Bamin — ommB(, p)). 
Cia: 
Now use the formula, valid for any two vectors P, Q, 


(a, P)(a, Q) = (P,Q) + ie, P,Q), 


where (oc, P, Q) is the triple scalar product. Taking the p’, p” matrix 
element of N, we get with the help of (25.4), 


e wr , 
, N a . DSP 
(p'|N |p”) Tae nT 


x [po'po” — (p’, p”) — ile, p’, p”) — m? — (a, p! — Pp’ )omm). 
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We are now going to use 
B= = —_ oe (0, 
yor 2, Oy Trot Ore by Ol OL gy OL 
With these results we get 


= DAE Ce Day 
Po + Po’ 
X [po' po’ — (p’, p”) + ils, p’, p”) — m7]. (27.12) 


(p'|a,Na*|p!") = 


Applying the ~ transformation to both sides we get, after some 
reduction, 


- / ap a k 2 
(p'|a,Nat!| p!’) = —eGpr_p EG SP oo 
Ho; D op ) ile. p sap 1)| ie 


(per op aoe 
‘s E 7 2p? 


I would like you to note that in the small terms of order p? we can 
carry through the averaging over k directions already at this stage. It 
is not permissible to do that for the larger terms, because we shall 
later have to multiply this expression into other factors containing k. 
So we get 


(p'|a,Na*|p’) = — 7 Oe [ae = 80) 46 a ea 
k? 
x ( ES i) ae BN = I Wf 
Now we form 
(Plo. Yoet|p”) = oe! {ae = 0 + ee 
k2 

x ¢ ay 5) eetlax p= Dp W) 

Xo Dn 1 eee mae 
Working out the factors in the last line to the first order, we get 


(P'la,Yoa'|p”) = © --)@AT AT? + wo) — @! + pp’, kK) 
X [G7A)0 Gee ean 
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We have to average over all k directions and that leads to 
ne [e@? =p)? + 7 ey) ] k* 
p2r p!'—p' 2 p p ) i(o, p 3 p )) = 3u2 
xX (T+ w") — Bile, p’, pK? 
* (Pay? + Fy? + wy} 


The integration over all A-values between the appropriate limits gives 


(p"| Yo|p’") = (4r)~ if (p’ lay, Yoo" [p’ ne if 
= CApi_p’ , 2 61 =e a ae 
5 e772 EP) ig top (ar) § log 2 + 355} 


+ i(o, p’, p’”’) log (2v)~"}. 
We may write that result symbolically as follows. We have 


(p'|V?a|p") = —(@’ — p”)?@pr_p. 
Also 


(p'|grad @|p”) = i(p’ — p’)Qpr_p’, 
SO 


(p'\(o, grad @, p)|p’’) = i(o, p’, p")@pr_p. 


Thus we get the result that 


Vo = = 


an? 4 log (27) + 4log2 + 3% 


ae sa log (2Y)~*(o, grad @, p). (27.13) 


Now let us take Y3 given by (27.7). We can write it 


Ys p*(1 + M28) (p @ = £7 reQ) 


ss = om n( BD! — (Df) ‘e@) ity “ (D"@ = GD) 


D 
Y3i + Ys2 + Yo3- 


Now let us take the first term here. We have to multiply it by a, on the 
left and the right and we get from (22.1), 


~ 2 2 m 
a,Y3 io" = rat — 2renl) (p Do — E; — eQ). 
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Now we can use the weak equation 
Do =< = Py ae 


We can use weak equations here because we have already multiplied 
by a, on left and right. With the help of the further weak equation 
(27.2) we get 


S a 2VanMN = 
asa rat ae ) Nive) : 
Dek Z(1 2a) N(p0)7?. (27.14) 


I have here used the zero-order approximations D = ) and fo = uy. 
This is permissible because N is a very small quantity, of order p@. 
If we kept higher-order terms, the next approximation would involve 
terms with (p, k) and these terms cancel when one averages over all k 
directions, and then the next terms with (p,k)? will be too small, 
because they will give a contribution to (27.14) of order p?. So (27.14) 
is valid after averaging over k directions and neglecting p® terms. 

From (27.3) the first term in the parentheses in (27.14) contributes 
zero to the average and we are left with 


et Ay = 
ay Ys3 ya" ~ ec eo y 


and 
G 3 


2, re dk 
Ys} = (47) u / au Yara ~ 2 log Qa) a Ning me (275) 
: Ym 
Let us go on to Y33. We have the equation 
(ilp*eali) = G|p?(U — Eli) = (i\(U — E)p? |i) = (ileap?|i). 
So Gp? ~ p’@. 


Now let us consider the quantity D~?: 


De ee 1 — v8 + p? terms | (27.16) 


By the same argument as before, we get D~?@ ~ @D~? after averaging 
over all k directions. Using this result we can write 


a UY33e" = 2e(D-?@ — D~'@D~!) 
~ e[D*@ — 2D—'@D + @D-?]. 
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We take matrix elements and get 


(p' lex gso|p"”) ~ e@yr—p[(D-? — 20D'D'")-! + (D')-?} 
= Gp» [(D)-! — (D")-Y?, 


and, retaining only the first-order terms in [(D’)~' — (D”)~!) given 


by (27.16), 
a , fon ak 2 
(Pla Yasa" |p") ~ ep —p! aa x 


If we average over k directions we get 


NIN 
(p' la, Ysga"|p’) ~ Ae Apt —p! a AA? : 
so that 


~ ek? 2 
ay Ysa" ~ — 342 VASE 


Now we can integrate 


g 
ae ak e 
Y33 = (47) a ay Yzgar" Tk ~ Tie 
X [log (2v)~* — 4log2 — iJ V’a. C717) 


Let us go on to Y32. As we did before, we can deduce 


(D?fo) '@ ~ Q(D*h)“', 
and so 


a,Yzca" ~ 2Vamm[2D~'BD~! — (D?ho)~'e@ — e@(D7fo)—')). 
Taking matrix elements we have 


(P' |e Y3o0"|p’”’) 
ane hat 


- 4 
x (D0) : 


(po! + po") D’D" 
2VOmMC Ap! —p’ 
Po = Po 
ee (Dyer ee  — po’) 
Se DS a Oe 


+ (D"*p0)| 


136 Lectures on Quantum Field Theory 
and to second order in (p, k) we have 


2VOimMNER yp’? p! 


Bo’ + Bo" 
yx 200’ — 1)” | @ — Pp”) 
pear Mh 
x (CH, 20K) _ (WK) _ 200") 
pPr2 ur? ue r2 p2r3 
VOmMEA pt 


a dia ay p” ky 
iM 


2 Z 
x (4 a utr? a7 As) : 


Averaging over k directions makes 


(p’ — p’, k)? > 4@p — p”)?k?, 
SO 


oe ge yhe o ] 2 
Y 


We integrate and get 


g 
= ae k pee 
1a) — (Gn) ae ayYar! Tk) ~ Hp Los Qn =) (ats) 

Well, now that we have done all the integrations, it just remains 
to collect the terms (27.15), (27.17) and (27.18) together, 


2@ 


+2 €,. Qye a, NCGioim 
Adding on Y, and Yo, given by (27.10) and (27.13), we get 


v2Q@ 
“—s vs log (2Y)~* + 4] on 72 08 (27) (6, grad @, p) 


te [log OY =e 4]JanN(po) 


In the last term in Yp the factor (po)! can be replaced by m7! 
without affecting the order of accuracy. This term can be connected 
with the first two terms. We have from (27.11), 


VE eo 


2N = PobPo a (a, p)Bla, Pp) al m?B ele a,n[ Bla, Pp) ai (a, p) 8]. 
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Now 


(a, p)B ~ (E; — anm)B ~ B(E; — ain) ~ Ble, p), 


SO 
2amiNN ~ [poBpo — (a, p)B(a, p) — mB] An IN 
As [PoBpo ay (a, p)B(a, Pp) pias m? BYE; ai (a, p)] 
ae [PoBPo _ (a, p)B(a@, Pp) aa m” Byn, 


as the (w,p) in [E; — (a, p)] gives terms of order @p® and is negligible. So 


2e wt ? 
2p'lamN|p”) ~ : Pp <a Jn ’ WE 8 pen 
(p | |p ) Po a8 Po [Po Po (a, p Yq, p ) m1 ] 
ea Dead 2 e 
ots = 00 = ICA 
So 
ee, 2 
AanN~ — res vad am (o, grad @, p). 


Using this result we have 


Yr~ — 


ey? nag, eee 

3m Y cilos (2s =| ee (o, grad @, p). (27.19) 

That is the final formula for the relativistic part of Y. It has these two 

terms, one of them involving V@, the other involving the electron spin. 
You note the appearance of Y in this formula; we can’t take Y 

equal to zero, so we have to make a cutoff for low values of k. 
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Here we use quite a different approximation, an approximation that 
is appropriate to nonrelativistic electron theory, because here we are 
using & values that are small compared with the rest energy of the 
electron. 

The nonrelativistic part of Y is defined as 


Ym 
os 2 Soe eg 
Yvr = ey f a(Y — Yo)a" : 
0 [k| 


with 
vy + (U/po) 


ee + (E/|E)) 
Poole (k| — VPo 


| (eae | 


Yo = 
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The first thing to note is that we get small denominators only for 
y = 1. With regard to the terms when vy = —1, there are no small 
denominators. We might have integrated the relativistic part for those 
terms right down to k = 0. The contribution of the terms with 
vy = —1 to the nonrelativistic part is negligible. So instead of having 
to average over both values of v, we just put vy = 1 and put ina 
factor 4. Thus we take 


1 + (U/fo) 


oJ tt 
2° Pe 4. |Ki = ea 


7 (f =e oe on 


Now, from (25.5) and (27.12) we know that 
ee 2 
ia G a eg 
There is not very much difference between the ~ quantities and non-~ 
quantities now, because k is small; the difference between them is of 
order k only. We can thus replace U/jo by E/|E| in Yo without 


disturbing the accuracy of our calculations, the error being of order 
@p?/|k| in a,Yoo", giving an error in Yyx of order 


Ym ak 
ra | \k|-'—— = 4ryp’a, 
0 Ik| 


which is negligible because we are taking y < 1. 
So now 


a(Y — Yoa* = ays ¢ a =| 
x (UK| 2) Se (ke eee) 


Expanding in |k|, the above expression becomes 


On E 
2 ( i al 
ii = [El = E: Be Sere Bo — Po F 
HL Welclel + 1B] — Bs) TKI TRICKY F Po = Pod} 
(28.1) 


For » = 1, 2,3 we have a, on the left and right, and the matrix 
elements of the velocity a, are of the same order of magnitude as p/m 
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in a nonrelativistic theory. To see this we write 


(ile,|n) = CE; + Bye Gl Boy + a,E|n) 
= m~\(ilp, — eGa,|n) = m—{ilp,|n), (28.2) 
because 
(ile@a,|n) << miila,|n). 


So we see that the velocity has approximately the same matrix elements 
as the momentum divided by the mass. The momentum is counted as 
small in this work and so for up = 1,2, 3 we have a small momentum 
occurring to the right and left. The last term in the [ ] in (28.1) is 
then negligible and we are left with 


aera 5 E |E| — E; 
a(Y — Yojo” = ta ¢ de =) = a, (28.3) 
LE AGS) 42 ie = ea 
For » = O we do not get this small momentum, so the term for 
#2 = 0 1s substantially larger than the other three terms, and we have 


to calculate it with higher accuracy. To do this we use 


EC Er . iy — .)2 
9-%= 304+ SZ) 0b - Ae y eee 


(E\/LkI |k|? Ik\3 
7 (|£| — E,)° _ ra Po — Po 
Iki@(k| + |£| — 2; — {ki (k|? 


_ (@o — Po)” au (Bo — Po) |. 
[k|? Ik|?(k| + fo — Po) 


The term in (fp — po)? is negligible because it involves p®. Also one 
can show that 


Cyr) ele) and 5 — po) — (2 — E)- 
make negligible contributions. To prove that, we see from (25.1) that 
(\E| — Es) — (bo — po) ~ — (UB + BU) + 3(UB + BU) 
= O(@k). 


When multiplied by [| + (£/|E))] and averaged over all k directions, 
this gives a quantity of order @k*. We must then multiply by 
\k|~2(d°k/|k|) and integrate, so we get a term of order 


Yn 
a | |k| d|k| = O(a), 
0 
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which is negligible. We can carry through a similar argument to show 
that (Bo — po)? — (\E| — £,)? gives a term of order Y@ and is also 


negligible. 
Now we are left with 


a hn (\£E| — E,° 
— 1+ 5) = 
ek al (é\/ (kPdk| + |E| — Ed 


= ~1(8| - By(1 +; a) 


2) = a ee 
\k[8(\k| + |Z] — £2) (|E| — Ei). (28.4) 


All that I have done here is to rearrange the factors, and we can do that 
because all the factors commute with each other. The only things we 
have here that are not numbers are E£ and |E|, and they commute with 
each other. 

We may replace |£| by £ in the factor on the right of (28.4) because 
we have got the factor 1 + (£/|E|), which secures that E = |E], 
because when £ = —|E| it gives zero. We thus get on the right of 
(28.4) the factor 


E~ E, = pth) + anm — e@ — E; ~ (@,k). 


When we take averages for the state i, we can treat the factor (|E| — E,) 
on the left of (28.4) in the same way, so that it also becomes (a, k). 
So the term (28.1) for u = O has the same average as 


E E— E; 
Seavey ae — a, k). 28.5 


In the expressions (28.3) and (28.5) we can drop the ~’s, because 
the a, are already small (their matrix elements are small) and the 
difference between the ~ and non-~ quantities involves k, which is 
also small. When we now average (28.5) over all k directions we get 


= qa Fp): 
ee el + alge ope yy 


Adding on (28.3), we get 


a i iE [El ges 
ou Gort ~ Sor (1+ te) ge Ter ED 
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Thus 


Ym Ym 
= :., re d°k 
Mon == (GS 4) au — Goa FE ~ sf e (: a 5) 


ae aaa 


TEE [ps] Sook 


We carry out the integration and get 


ee a: VG ae 8], = 1a 
Ynr bar (1 as i) (2 E;) Ea 


Now with sufficient accuracy we can neglect |E| — E; in the 
numerator of the log, because Ym is large compared with the energy 
differences that we are interested in, so 


(i| Ywrli) = 5 >> (ile,|n)(E, — E;) log =——4 a E; (nlar|i) 


cl are 


2 J, = 2 ym 
= —— 7 — |i pr|n)|~ log ———— > 
from (28.2). 

This expression must be evaluated by computation for the particular 
problem in which one is interested. The result can be expressed in the 
form 


Gl ¥weli) = = BS lpm? 


En>0 


where M is some average of the |E, — E;|, which must be obtained 
from the computation. The computation must, of course, take into 
account the values of the matrix elements that multiply the log term; 
for various n values. M occurs only in the logarithm, so variations in 
M are not too important. Heitler} gives the value log (m/M) = 7.6876 
for the Lamb shift in hydrogen with the main quantum number 2. 
The restriction £, > 0 is not important, as the terms for £, < 0 
would contribute only a negligible amount, because of the smallness 
of the matrix elements for E, < 0. Dropping this restriction, we can 


J W. Heitler, The Quantum Theory of Radiation, 3rd ed., Oxford University 
Press, London and New York, 1954, p. 346. 
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express the result in a more suitable form: 


P : 2 WH as ; 
G| Yyrli) = 3772 8 aye (ilp(LE — £;)p|i) 


2 ym : 
=> >7 == Ui = — ES re 
52 8 ag Milp(U — e@ )prli) 


= i ye ee ar Mil(U — Eidpy? — 2preap, + pe (U — Eili) 


1 WHE 7. 2 2 . 
= 5,2 log py Wileapr — 2pre@pr + Pr eQ|i) 


— slog = (ilv@li). (28.6) 
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So finally, adding the expressions (27.19) and (28.6) and the second 
term of (26.3), we get 


GVW) + (lz) = ~ 385 Wal (los Ht, + 5) 


zy 
é ‘ i 
= 4m2 (ile, grad Q, p)lZ). 
and 


2 
= (il Y + Z|i) = Lamb shift. 


The important thing to note about this expression is that Y has 
disappeared from it. Of course, Y has to disappear if we have not 
made a slip anywhere, because the precise value at which we make the 
distinction between the relativistic and nonrelativistic parts is un- 
important. Y must be of the right order of magnitude for the approxima- 
tions to be valid, but just where we chose Y doesn’t matter. 

Well, that gives the calculation of the Lamb shift up to a certain 
degree of accuracy. People have carried the calculations through with 
higher accuracy, working in the Schrédinger picture. Professor Lamb 
told me recently that he is now doing experiments again, with a dif- 
ferent apparatus, and hopes to improve the measurement of the Lamb 
shift by quite a substantial factor, and then it will be well worthwhile 
to have the theory of the Lamb shift worked out to a greater degree 
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of accuracy, to compare it with these expected experimental results. 
Of course, the calculation will be very complicated, because it involves 
working to order (e”/fc)? or maybe even to (e?/hc)? to get the calcula- 
tions with the high degree of accuracy with which the experiments will 
be performed. Also, we shall have to stop doing the other approxima- 
tions we have made, counting the momentum as small. In all our work 
we have counted the momentum of the particle as small and worked 
only as far as p”. It will be necessary work as far as p*, and it may be 
necessary to work as far as the square of the potential @, so that all 
together there is a lot more work to be done to get the calculation of 
the Lamb shift to the desired degree of accuracy. The treatment I 
have given shows the method and shows that one is not bothered by 
infinities if one carries through the calculation correctly. 

In our calculations of the anomalous magnetic moment and the 
Lamb shift we were dealing with one-electron problems, and the 
essential point in these calculations is to count the interaction energy 
as small. We counted { 4,j" dx as small, and it seems that this is a 
good approximation for one-electron problems. 

In a way this is rather surprising, because the quantity we are 
assuming to be small is not gauge-invariant. Our whole Hamiltonian 
is a gauge-invariant quantity, but we split it up into two parts, a proper- 
energy part and an interaction-energy part, and those two parts by 
themselves are not gauge-invariant. It is not altogether satisfactory to 
do that sort of thing, but it does seem to work for one-electron prob- 
lems. One would like to have a calculation that works entirely with 
gauge-invariant quantities; if one could do that, the theory would be on 
a more satisfactory basis. 

Another thing you should note about these calculations is that we 
have not made any use of the weak equations 


oA 


v0! (29.1) 
abe 


It seems that these weak equations do not play an important role in 
one-electron problems. 


830. The Coulomb Force 


If we proceed to problems with more than one electron, the situation 
is quite different. If we have more than one electron present, it is not 
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permissible to count Hg as small, because we have Coulomb inter- 
actions coming into play. Then we have to make use of the constraint 
equations (29.1). 

I should like to show you how the Coulomb interaction does come 
in. It involves a certain treatment of the longitudinal electromagnetic 
waves in the Hamiltonian. This work is given in all the books on the 
subject in connection with the Schrédinger picture, and it is called the 
transformation from the Lorentz gauge to the Coulomb gauge. I 
should like to show you that one can do corresponding work with the 
Heisenberg picture. 

The constraints (29.1) give us two conditions on the Hamiltonian 
variables at a particular time, (16.6) and (16.7). Of course, there is 
one of these constraint equations for each point in three-dimensional 
space. 

Let us pass to the Fourier components. We write 


Ae = / (Ait oe As en Se 
ee 
lie = fine ak, 


Equations (16.6) and (16.7) written in terms of the Fourier coefficients 
are 


i / KA =A em ee a (30.1) 


iIk|(A ox ~ Ao—z) — ik,(Arg + Ary) ~ 0, (30.2) 
Iklk-(Arz — Arn) ~ |k|?(Aox + Ao—x) ~ —4ajon- (30.3) 


Take—[equation (30.3)]—i|k| [equation (30.2)], and that gives us 
2\k|?A on — 2Ik|kpAre © 4mjox- (30.4) 
Let us put ky = |k|; then (30.4) becomes 
A ee (30.5) 


The conjugate complex equation is 


Auk ~ ko 2° (30.6) 
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Now we split up the electromagnetic field into longitudinal and 
transverse parts and get an expression for the longitudinal part of the 
field energy. We put 


ee An, A Ark 


where 

AL, = (kr/ko”)ksAsy by definition, and A, = A, — AX. We have 
k,A™, = k,As,, and therefore k,A?, = 0, which checks that the 
definition of the longitudinal part is correct. 

Let us go back to our expression for the electromagnetic field energy 
expressed in terms of the Fourier components, (15.16), and ignore the 
infinite constant, which means fixing it so that there is no zero-point 
energy. Wecan write it Hp = Hp, + Arr, with 


Hp, = an? | ko®(Ahial — AorAon) ak, 
Herr = 4r° ii fone AD dak. 


The cross terms containing 44,47, cancel out. This simply means 
that there is no interference energy between the longitudinal and 
transverse waves of the field. 

We can write Hry as 


Hry = 42° [ths AnnGvon — k*A,x) 


—koAonlk" Auk oe iewsla))| d’k 


An? / (Coa tee A, — kp Aone 4 uel cc. 
(30.7) 


Now we have to make use of the constraint equations. I discussed 
in Section 14 how we are to use constraints in the Heisenberg formula- 
tion; we have the rule that in general we can multiply constraints only 
by quantities on the left and not by quantities on the right. For the 
present work we find that we can get all the results we want if we. 
multiply the constraint equations only by quantities that commute 
with them. Of course, in that case it is permissible to multiply them on 
either side. 
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(Aox + Aox)k* (Aue + Aye) or + Aox)Vox + Jo—x) 
and 

(Aor — Aow)k"(Aue — Aun) ~ ~ do — Aox)Vox — jo—x)- 
Subtracting 


= = 2 , = « 
Aork* Aur + Aork"Aur ~ = (Aor jo—z + AoxJox)- 


Also, since k’A,, and k¥A,;, commute, we can infer from (30.5) and 


(30.6), 
es 


= Ar : 
k Ayk Aye = ees JOR Ome 
k uk jag J 0kJ0—k 


Using these results we find from (30.7), 


Ary ~ Bact | Cine ae + Aorjor) dk + 16a f tobias d’k. 


Then if we go back from the Fourier components to the field quan- 
tities, using (30.1), we get 


: 1 OatOnt 
ela mee d°x + 5 ||) bees ad 


The term — I Aozjox d°x cancels with one of the terms in the inter- 


action energy Hg, and 
ff Jado pax 2x! 
|x — 


is the Coulomb interaction energy between any particles that are 
present. In that way the Coulomb interaction energy between the 
charged particles appears from the longitudinal part of the field 
energy. 

The important point of this calculation is that you can carry it through 
in the Heisenberg picture using the constraint equations only to the 
extent that you multiply them by quantities that commute with them. 
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§31. The General Physical Interpretation 


Another problem I would like to consider now is, what is the general 
physical interpretation for this theory? In the usual quantum 
mechanics one takes the Schrédinger wave function, normalizes it, 
and forms the square of its modulus and interprets what one gets as a 
probability. We cannot do that anymore because we have not got a 
Schrédinger wave function in this quantum field theory. What can we 
do instead? 

In our calculation of the analomous magnetic moment and the 
Lamb shift, we had a normal ordering of our g-numbers. To make 
this normal ordering appear naturally in the calculations, one should 
introduce a certain ket vector at each time. At each time /, introduce 
a ket, which we may call the vacuum ket at time ¢ and denote by |O,), 
satisfying 


(any absorption operator at time 7) |O,) = 0. 


This |O;) would seem to play the role of the vacuum ket, but of 
course it differs very much from what one usually understands by the 
vacuum ket, because it is only an instantaneous thing; the ket that 
satisfies those conditions at time ¢ would be quite different from the 
ket that satisfies those conditions at a later time. People working with 
the Schrodinger picture assume that there is a ket that satisfies the 
vacuum conditions at all times, but such a ket does not exist. However, 
the vacuum ket at a particular time is a well-defined concept in the 
Heisenberg picture. 

I am now going to make two assumptions for the physical inter- 
pretation of the theory. 

Assumption 1. Each physical state corresponds to a q-number K 
which is a constant of the motion in the Heisenberg picture. K will 
involve the time explicitly, as well as the dynamical variables at time 7. 


dK = 


Hic Ee. 
we 0 Or i, + KH — HK = 0. 


In fact, we have been working with such a K all along, during our 
integration of the Heisenberg equations of motion. 

I would like you to note that this assumption differs very drastically 
from the usual assumption. It assumes that a physical state corresponds 
to a q-number and not to a ket vector. 
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Assumption 2. Anything that can be observed for the state K at 
time ¢ is described by K|O,). 

With this second assumption the normal ordering of K becomes 
important. This normal ordering involves putting all emission operators 
in any term to the left of all absorption operators. After the normal 
ordering has been carried out, those terms of K that contain any 
absorption operators will contribute nothing when K is multiplied 
into |O,). According to assumption 2, these terms wiil be latent at time f. 
The only things that show up are those terms in K involving only 
emission operators. 

Each of those terms will refer to certain particles in certain definite 
states. We can take the coefficient of each such term and form the 
square of its modulus. We then get a definite number and we can call 
that number an “intensity” for the appearance of these particles in 
the corresponding states at time ¢. I use the word intensity here and 
not the usual word probability, and there is a reason for that. These 
intensities cannot be normalized; or if you do normalize them at one 
time they will not remain normalized. Of course, they cannot remain 
normalized, because any latent terms that there are in K at one particular 
time may cease to be latent at a later time and contribute to observable 
things. 

This interpretation does not give such definite results as the usual 
interpretation with the Schrédinger theory because we do not have 
things we can definitely nail down as probabilities. We can only say 
that if the intensity is large, the probability of observing the particles 
is large; if the intensity is zero, the probability is zero. 

You will note that there is less determinism in the present theory 
than there is in the usual theory involving the Schrédinger wave 
function. The reason for this is that the quantity K|O,), which gives 
all that is observable at one particular time, does not fix what is observ- 
able at a later time, because of the latent terms in K—latent terms in K 
which you could provide with any coefficients you like and they will 
not contribute anything observable at time ¢, but they will still influence 
observable things at another time. 

In this theory there is no S-matrix. The S-matrix is something that 
is unitary and gives probabilities. I just do not know how to define an 
S-matrix working with the Heisenberg picture. Maybe future develop- 
ments will show some way of introducing an S-matrix, but the usual 
way of introducing it in field theory involves so much departure from 
logic that I do not see how one could take it over into a logical theory. 


32. The Present Field Theory and the Usual One 149 


The quantum field theory, which J shall leave in this state, is therefore 
far from being complete. There are various further ways in which one 
needs to develop it. One of the most important of these, of course, is 
to find some way of dealing with strong interactions. Electro- 
dynamics works fairly well, because there the coupling constant is 
small and we can set up a perturbation theory in terms of it, but for 
strong interactions the coupling constant is not small, and one cannot 
use the Hamiltonian theory with this kind of perturbation method. 
The main problem for the future development of the Hamiltonian 
theory is to find some method of approximation that will enable one 
to apply the Heisenberg picture to collision problems when the coupling 
constant is not small. There are many physicists who believe that the 
Hamiltonian theory breaks down for these problems; I do not think 
that there is any evidence for this breakdown; instead I think that the 
breakdown is only in the perturbation method. 

There is another unsolved problem in connection with electro- 
dynamics. I have done the one-electron calculations starting off with 
a zero-order approximation consisting of the emission operator of one 
electron. One might do a corresponding calculation for a one-photon 
problem, starting off with the zero-order approximation consisting 
of the emission operator of one photon. 

Well, if one does that one runs into trouble. One gets a self-energy 
for the photon which is quadratically infinite. If we bring in the cutoff 
it means we have the self-energy involving g”, an extremely large self- 
energy for the photon. This is incorrect, because, among other things, 
it violates gauge invariance. People usually just say that this self- 
energy should not be there and ignore it for that reason. 

People have worried a good deal about how, starting with a gauge- 
invariant theory, one can get a result that is not gauge-invariant. But 
from what I said before that is not at all surprising, because in these 
one-particle problems we use a perturbation method in which we assume 
that a quantity which is not gauge-invariant is small. There we have 
an assumption that is not gauge-invariant. It is surprising that this 
assumption works so well for one-electron problems, and it is not 
surprising that it fails for one-photon problems. 


§32. Relationship of the Present Field Theory to the Usual One 


The usual formulation of quantum field theory makes extensive use 
of the concepts of ingoing fields and outgoing fields. The connection 
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between these two sets of fields provides what is called the S-matrix. Let 
us examine how these fields are defined, to get an opinion of the value 
of this method. 

We define free fields to be fields that satisfy the field equations in 
the absence of interaction. Thus a free field quantity, B say, satisfies 


te Pa (32.1) 
dt 
where Hy is the sum of the energies of all the fields separately, without 
any interaction. 

We may consider the free fields which are equal to the actual fields 
of the Heisenberg picture at some standard time 7. This fixes them 
completely. Let us denote such free fields at time ¢ by Byy. Thus Bry 
equals the actual Heisenberg field at time 7. 

The free field variables B,7 for a given t and T must satisfy the same 
commutation (or anticommutation) relations as the variables Brr, 
since the B;r are connected with the Br,r by a unitary transformation 
built up from the infinitesimal unitary transformations described by 
(32.1). Thus the B,7 satisfy the same commutation relations as the 
Heisenberg variables at a given time. This result holds for any inter- 
action energy, even one with a cutoff. 

Let us suppose the variables B,7 go to definite limits as T > —o. 
These limiting values are the ingoing fields. They can be expressed as 
the actual fields minus retarded fields, the latter being similar to those 
that occur in classical electrodynamics. 

Of course it cannot be strictly true that the B,7 tend to definite 
limits as T —» — w, because the B;,7 will always be subject to substantial 
variations as T varies, but for T very large and negative it may well be 
that the variations are of so remote a character that they do not affect 
the calculations one is interested in. Under these conditions the ingoing 
fields become well-defined concepts. The outgoing fields are the cor- 
responding quantities with T— o. 

If the ingoing fields are well defined they are very useful. Their 
Fourier amplitudes are constants of the motion and may be taken as 
the basic variables in terms of which one integrates the Heisenberg 
equations. They have simple transformation laws under Lorentz trans- 
formations, so they lead to the possibility of the theory being set up in 
a manifestly covariant form. 

However, in the case of electrodynamics the ingoing fields are in 
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general not well defined. One can see this by taking the simple example 
of two electrons scattering one another, the force between them being 
mainly the Coulomb force. When the electrons are far apart the wave 
functions representing them do not tend to plane waves, because the 
Coulomb force falls off too slowly with increasing distance. The in- 
going waves for this scattering problem have definite frequencies and 
directions, but not definite phases, so they are not completely definable. 

For fields involving forces of interaction that fall off more rapidly 
than the inverse square law this difficulty does not arise. One can 
proceed to build up the theory in terms of the Fourier amplitudes of 
the ingoing fields. These Fourier amplitudes will consist of operators 
n of creation of ingoing particles and the conjugate operators 4 of 
annihilation of ingoing particles. One can keep to the usual physical 
interpretation of having the states correspond to kets, a state at all 
times corresponding to a ket in the Heisenberg picture. The vacuum 
state at all times then corresponds to the ket |0) satisfying 


for all %, and a general state at all times corresponds to the ket f(»)|0), 
the function f being some power series. 

One cannot set up a general quantum electrodynamics on these lines. 
One must avoid the use of ingoing fields. The present book shows how 
this can be done. One is forced to a new physical interpretation, but 
it is really not so different from the above one, because the K defining 
a state at all times takes the place of the f(n) above, both of them being 
constants of the motion in the Heisenberg picture. 

It is, of course, a severe handicap not to be able to use ingoing fields. 
In particular, one is prevented from setting up the theory in a manifestly 
covariant form. Maybe future progress will consist in making some 
generalization of the concept of ingoing fields so that it is applicable 
even to quantum electrodynamics. 
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